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Experimental Investigation of Influence 
of Edge Shape on the Aerodynamic 
Characteristics of Low Aspect 
Ratio Wings at Low Speeds’ 


G. E. BARTLETT? anp R. J. VIDAL? 


Cornell Aeronautical Laboratory, Inc. 


SUMMARY 


Some of the more significant results of an experimental in- 
vestigation of the influence of edge shape on the aerodynamic 
characteristics of a family of low aspect ratio wings (<R > 2), 
having straight trailing edges and taper ratios between zero and 
one, are presented. 

The results are shown to correlate with the approximate linear 
lifting-surface theory of Lawrence and with the semiempirical 
nonlinear theory of Flax and Lawrence. The discussion centers 
around (1) the effects of various leading and side-edge shapes, (2) 
the influence of leading edges, and (3) the effects of retreating 
side edges. Implications of the results as regards such design 
factors as lift-curve slope, stability and control, and leading 
edge suction and drag are discussed. In addition, the possibility 
of developing a new type of roll control for low aspect ratio wings, 
based on edge-flow-separation effects, is considered 


INTRODUCTION 


pr TO THE ADVENT of the transonic and super- 
sonic airborne vehicle, there was relatively little 
need for the low aspect ratio wing, and consequently 
only a minimum amount of research was forthcoming 
Requirements in the field 


on such configurations. 
of military aviation over the past few years have, how- 
ever, been responsible for the introduction of the low 
aspect ratio lifting surface as a desirable and necessary 
component of high-speed aircraft and missiles. 

Presented at the Aerodynamics Session, National Summer 
Meeting, IAS, Los Angeles, June 21-24, 1954. 

*This paper is based upon research carried out under Con 
tract Number AF33(038)-17397, E. O. Number 465-2 BR-1, 
administered by the Aeronautical Research Laboratory of the 
Wright Air Development Center. The project was under the 
technical direction of Fred L. Daum of that Laboratory 

t Head, Applied Aerodynamics Section, Aerodynamics 
Research Department. 

t Associate Aerodynamicist 


Considerable effort, both experimental and theoreti- 
cal, has been expended on low aspect ratio configura- 
tions at supersonic speeds, and methods have been 
developed that afford rapid and satisfactory predictions 
of the aerodynamic characteristics of these configura- 
tions. The aerodynamic behavior of these same con- 
figurations are not as well understood at subsonic 
and transonic speeds. Many of the design criteria for 
transonic and supersonic aircraft apply to the subsonic 
regime. For this reason, satisfactory methods for pre- 
dicting the aerodynamic characteristics of low aspect 
ratio lifting surfaces are of interest. 

The present paper includes experimental results to- 
gether with theoretical correlations for a series of low 
aspect ratio wing tests, run at low speeds. These tests 
were designed to attempt to fill some of the gaps in 
available data disclosed by previous investigations at 
Cornell Aeronautical Laboratory. 

The investigation included studies of the effects of 
variation of wing plan form, edge shape, Reynolds 
Number, and aspect ratio on the linear and nonlinear 
aerodynamic forces. In addition, pressure distribu- 
tions near the wing edges were obtained to aid in the 
analysis and interpretation of the data. Most of the 
wings tested had straight trailing edges; one series, 
however, had retreating trailing edges (trailing-edge 
span smaller than the maximum span). The questions 
investigated were mainly in connection with viscous 
effects on edge separation which may result in non- 
linear increments of lift. The bulk of the discussion in 
this paper pertains to the implications of sharp and 
round side edges, leading-edge influence, and retreating 
side edges on such design parameters as lift-curve slope, 
stability and control, and leading-edge suction and drag. 
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BACKGROUND INFORMATION 


A pronounced nonlinear relationship between aero- 
dynamic coefficients and angle of attack is typical of 
wings of nearly all plan forms when the aspect ratio is 
less than approximately 3. Since the usual lifting- 
surface theory predicts linear relationships, it is appar- 
ent that new expressions must be developed before low 
aspect ratio wing aerodynamic coefficients can be pre- 
dicted with suitable accuracy. 

Bollay,! Weinig,? Kriesis,* Scholz, and others have 
presented nonlinear lifting-surface theories for low as- 
pect ratio wings. Although the predictions from these 
theories are in accord with many experimental results, 
none are adequate for general usage. Recognizing the 
viscous origin of the nonlinear lift, Legendre,’ Adams,® 
Brown and Michael,’ Edwards,* Cheng,’ and others 
have recently attempted more rigorous solutions to this 
problem by representing the vortices due to edge sepa- 
ration as free vortices embedded in a potential flow at 
some position above the wing. By making suitable 
assumptions concerning the forces acting on the free 
vortices and the vortex sheets connecting the separa- 
tion point and the free vortices, these authors show that 
the viscous lift is of the order, (AR)'’*(a) 

These results presuppose a knowledge of the point at 
which separation occurs, and for practical purposes, 
are necessarily restricted to sharp edged wings. 

Flax and Lawrence,'’ in 1951, examined the available 
experimental data on low aspect ratio wings and ob- 
served that there were unmistakable correlations be- 
tween the characteristics of such wings and the geom- 
etry of the wing leading and side edges. Since none 
of the theories previously mentioned allow edge geom- 
etry as a parameter, they are limited in applicability to 
the class of airfoils for which they were developed. The 
theories of references 1—4 are further limited in that 
they offer no means of predicting moment characteristics. 

Since these early observations regarding edge in- 
fluences, Cornell Aeronautical Laboratory has been 
endeavoring to develop expressions that will have ade- 
quate versatility and be accurate enough for engineer- 
ing application. The C.A.L. approach to this problem 
has involved drawing on the linearized theory of a 
lifting surface in a steady flow for a basic term to which 
is added one or more nonlinear terms. Lawrence! 
has developed a theory suitable for predicting the linear 
aerodynamic characteristics for all plan forms of as- 
pect ratio below 4 and having straight trailing edges. 
Goodman'* has numerically evaluated the basic equa- 
tion of Lawrence’s theory and prepared charts of 
dC,/da and dCy/dC;,, for a large number of plan forms. 
Stone™ has carried out numerical calculations of the 
effects of flap deflection, aileron deflection, and dy- 
namic-stability derivatives. 

The development of general expressions suitable for 
predicting the nonlinear viscous effects that are inher 
ent in the aerodynamic characteristics of low aspect 
ratio wings is not in as satisfactory state as is the 


derivation of the linear expressions. The magnitudes 
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of the nonlinear effects are, however, large enough go 
that they can not be disregarded in calculations of low 
aspect ratio wing coefficients. 

The cross-flow drag concept suggested in reference 3 
and, earlier, by Betz’ for low aspect ratio wings, and 
suggested for use on slender bodies of revolution by 
Allen’® and Van Dyke,” appears to suffice providing 
adequate attention is given to the geometry of the 
wing edges. The cross-flow drag theory presumes that 
there is an additional normal force on each transverse 
cross section of a very low aspect ratio wing, stemming 
from the drag associated with the two-dimensional 
viscous flow around a long cylinder of similar section 
characteristics. The significant velocity is that com- 
ponent of the free-stream velocity normal to the wing. 

In applying viscous cross-flow theory, only terms of 
order a and a? need be retained, and the viscous con- 
tribution to normal force is taken as a constant times 
a’ rather than a constant times the more rigorous 
sin? a. The total lift coefficient on a low aspect ratio 


wing then is expressed by the equation 
Ce. (dC, da)jpa + Cp,.° l 


and the total moment coefficient (positive nose up 
about the apex is expressed by 


Cy = (dCy dC.)o (dC, da)o a — (& c)Cp,.a@ 2? 


where & is the distance from the wing apex or leading 
edge to the centroid of area, and c is the root chord. 
(dCy/dC_)o and (dC; da), are the slopes of the pitching 
moment and lift curves respectively at zero degrees 
angle of attack. 

The selection of a cross-flow drag coefficient, Cp,, 
which is appropriate for a given configuration, involves 
the consideration of many factors such as edge shape, 
plan form, thickness ratio, Reynolds Number, etc 
A thin sharp-edged wing of virtually any plan form will 
exhibit cross sections transverse to its fore and aft axis 
that closely approach a two-dimensional flat plate; 
consequently, one would be justified in selecting a 
Cp, of 2, the well-established cross-flow drag coeffi- 
cient for a two-dimensional flat plate. Furthermore, 
the sharp edges of the wing preclude the possibility of 
the flow behaving in any manner other than separating 
right at the edge, thereby removing the necessity for 
any consideration of the Reynolds Number influence. 
A thick wing having a varying cross section, on the 
other hand, should iiave its cross-flow drag force calcu- 
lated by integrating the drag force per unit length. 
Consider, for example, the delta wings tested by Lange 
and Wacke'’ which had NACA 0012 sections. The 
cross section near the nose of these wings more nearly 
approximates a body of revolution than a thin plate, 
while the cross section of the major portion of the 
wing may be likened to a plate of finite thickness with 
round edges. It can be seen that a good approxima- 
tion to the cross-flow effects could be obtained by ex- 
perimentally determining the cross-flow drag of each 
section falling between these two limits and by imnte- 
grating to obtain the average value for the wing. In- 
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INFLUENCE OF EDGE SHAPE 
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TABLE 1. Model geometric characteristics, 
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Aspect ratio 1.5 and 2.0 delta wings. 
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2. Typical wind tunnel installation aspect ratio 1.5 sharp- 
edged delta wing 
14 , 
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40 ASPECT RATIO |.5 DELTA WING | | | attack. It should be recognized that this assumption R 
WITH ELLIPTICAL EDGES ile is often erroneous, particularly for the case of wings — 
7 fe a with rounded edges. However, the cross-flow Mi post 
©RN2 3210 e , - Tag 
10 Fabs — coefficient for these wings is small enough to result jj SUBS 
Glen | small errors for small angles of attack before separatio —_ 
8 - a _— actually occurs. relat 
G Be oo THEORY | The fact that we are dealing with viscous flow whey plat 
6 ae examining the nonlinear effects on the aerodynamic = 
4 | coefficients of wings requires that the Reynolds Number war 
” ar | | of the cross flow be considered. It is mainly of im. be | 
P A a Ss | portance for wings with rounded edges since the nature ae 
"ad | | of the flow may well determine whether such an edge sai 
oO | | | | | acts as a round or sharp edge. Plates with sharp edges = 
-10 — a on do not exhibit any great variation in drag with Rey- drag 
| ia nolds Numbe tial | 
-8 7 1 = One might be inclined, and certainly entitled, to view alter 
| . . ‘ 3 - ae the application of the viscous cross-flow concept to _ 
* | pL ir wings with some reservations. Certainly the flow in pats 
“w | co al | | transverse planes cannot behave two-dimensionally - 
| KK | when the cross section is rapidly changing, as in a wide 
LINEAR THEORY ° a ; aes : 
-2 | delta wing. The leading and trailing edges must in 
Lo" | | fluence those portions of the wing near them. It seems 
oe = — sa a a as “5 rs even more incredible that there can be a coexistence o/ : 
a? the flow identified with the cross-flow drag concept arte 
Fic. 4. i.e., separated from the side edges, and the potential the | 
flow identified with lifting-surface theory. On the see 
other hand, one is forced to concede that even with some 
deed, this approximation would be useful if it were the introduction of a minimum of empirical constants, | : 
feasible to predict the value of Ua, the cross-flow ve- extremely good agreement between theory and experi 4) % 
locity, at which each section would first come under ment has been exhibited for both the lift and moment veins 
the influence of separation. If such a calculation could on wings and wing-body combinations. ' In addi- pln 
be made, a cross-flow drag coefficient as a function of tion, data from references 20-28 and others have been on 
angle of attack would be obtained, and as a conse- correlated with the cross-flow theory with similarly _ 
quence, it is expected that the use of such a coefficient good results. 
in predicting the nonlinear characteristics of a wing Edge 
would be laborious. Rather than attempt such a time- Fa een eee 7 ee le 
consuming calculation, it might be more desirable to | T | | ‘hone 
consider an ‘‘average’’ or effective cross-flow drag co- hPa a aa | 9) th 
efficient; that is, a coefficient which would not predict © RN=6 x10 | | | edges 
exactly the viscous effects at each angle of attack, but 10}— a alti 10% | stitut 
would approximate the effects at all angles of attack series 
with engineering accuracy. Recognizing the limita- lated 
tions of such an approach, in most cases it is not justi- of wi 
fiable to rigorously obtain an effective cross-flow drag aspec 
coefficient from chordwise integrations, but rather to hovel, 
obtain this average value by correlating experimental | These 
determinations of the nonlinear or viscous lift with physi 
Cp.a*. It is not surprising that the experimental lifts Fig. : 
of the Lange-Wacke wings show good correlation with solida 
theory for Cp, = 0.50 because this is of the order of an made 
average value between the turbulent drag coefficient of -8 ; can b 
a circular cylinder, namely 0.35, and the Cp, value of | from 1 
1.0 associated with a two-dimensional thick plate “6 } tical-c 
having round edges. For thin wings with rounded | ness-c 
edges, one would expect that this average cross-flow aes ae eee a a THEDRY + i cent, 
drag coefficient should lie closer to that of a two- ; | | a wings 
dimensional flat plate with rounded edges, and would | za | per ft 
therefore expect a degree of correlation with Cp, = 1. og 1 , L a “a the w: 
obtair 


The cross-flow theory as presented assumes that the 
cross flow separates immediately at small angles of Fic. 5 
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INFLUENCE OF EDGE SHAPE 


A plausible explanation of how there can be the 
seeming coexistence of the viscous cross flow and the 
potential flow requisite to lifting-surface theory is 
suggested in reference 10. The flow which has sepa- 
rated from the edges of the wing is assumed to have a 
relatively small component of velocity normal to the 
plate. This separated flow is carried in close proximity 
to the surface of the plate by the streamwise com 
ponent of velocity plus the downwash velocity. It can 
be likened, in most respects, to an addition to the 
upper-surface boundary layer and an effective addition 
to the airfoil thickness. The problem then becomes one 
of considering the wing plus the wake due to cross-flow 
drag as a new thicker rigid body embedded in a poten- 
tial how. Since this added thickness will result in no 
alteration of the formula for apparent mass, (pb*) 4, 
which is valid for a two-dimensional elliptic cylinder 
having any length minor axis, the slender-body poten 


tial lift should remain unchanged. 


DISCUSSION OF EXPERIMENTAL INVESTIGATION 


The present experimental program was initiated at 
C.A.L. in order to provide additional verification for 
the cross-flow concept, and in order to consider in 
greater detail some of the effects originally noted in 
reference 10. Among those effects investigated were 
1) the effects of various leading and side-edge shapes, 


2) the conditions under which a leading edge ceases 


to act as such and acts as a side edge, and (3) the effects 
of retreating side edges; that is, wings which have a 
maximum span at some point other than the trailing 


edge. 


Edge-Shape Effects 


In reference 10, an attempt was made to establish 
three general classes of edge shapes; sharp edges (Cp, = 
2), thin rounded edges (Cp, = 1), and thick rounded 
edges (Cp. = 1/2). In order to investigate what con 
stitute aerodynamically round and sharp edges, a 
series of six delta wings were tested and the data corre- 
lated with the existing cross-flow theory. 
of wings consisted of three aspect ratio 1.5 and three 


This series 


aspect ratio 2.0 constant thickness delta wings with 
beveled edges, elliptical edges, and thick round edges. 
These wings are illustrated in Fig. 1 and their pertinent 
physical characteristics are tabulated in Table 1. 
Fig. 2 shows a typical model installation in the Con 
solidated-Vultee wind tunnel where all these tests were 
made. Asa matter of interest, the separation vortices 
can be seen in this photograph as vapor trails emanating 
from the wing apex. The beveled-edge wings, the ellip 
tical-edge wings, and the round-edge wings had thick 
less-chord ratios (based on maximum chord) of 5 per 
cent, 7’ » per cent, and 10 per cent respectively. These 
Wings were tested at Reynolds Numbers of 0.9 & 10° 
per ft. and 1.S xX 10! per ft., and pressure data over 
the wing edges as well as six component force data were 
obtained. 
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The 1.5 aspect ratio wing lift and moment data have 
been correlated with the linear theory of reference 11 
and with the cross-flow theory previously outlined, 
and these correlations are shown in Figs. 38, 4, and 5. 
It should be noted that all Reynolds Numbers are based 
on the root chord. It can be seen that for the beveled 
edge wing, the cross-flow theory yields an excellent 
prediction for lift and moment for all angles of attack 
of practical interest, and that the results are essentially 
independent of Reynolds Number. The fact that 
Reynolds Number does not appreciably alter the data 
below the stall is expected since the flow must separate 
at the sharp edge of the wing regardless of the nature of 
the boundary layer. 

Examination of the lift and moment data for the wing 
with elliptical-edge shape indicates that for low angles 
of attack, these data correlate best with the linear 
theory and that the viscous effects are felt only at 
angles of attack greater than about 15 This indi 
cates that the cross flow does not begin to separate from 
the upper surface noticeably until relatively high angles 
of attack are attained. This conclusion is qualitatively 
borne out when the pressure distributions over the edges 
are examined (see Fig. 6). The cross-flow separation is 
manifested in a vortex motion of the fluid on the 
upper-wing surface, and consequently the center of the 
vortex or separated region is characterized as an unusu- 
ally high negative pressure region on the upper wing 
surface. It can be seen that this unusual negative 
pressure region is first apparent at angles of attack 
These pressure diagrams point 
It is observed 


between 15° and IS 
out two phenomena of general interest. 
that the separation first appears at the aft station 
: < 18°) and then becomes 


a ~ 


(x/Cmer. = 0.57, 15 : 
apparent at the forward station (¥ Car. = 0.25) at 
24° Sc aS 27 This indicates that the separation is a 
progressive phenomenon, and in this case its progression 
is relatively slow. For the beveled-edge wing, the pres 
sure distributions indicated that separation occurred 
almost simultaneously on all transverse wing sections. 

The pressure distributions also illustrate the fact that 
on the wing with elliptical edges, the separation does 
not occur at the transverse section extremity, but rather 
well inboard on the wing. The result is that outboard 
of this separation point the flow pattern is essentially 
that of a potential flow and that the peak negative 
pressures at the section extremity are preserved 
Since this leading-edge suction continues to exist after 
separation occurs, it is expected that the drag character 
istics will reflect this thrust component. 

The lift and moment data for the wing with the 
thick round edges are seen to correlate well with the 
cross-flow theory (Cp | 2) to high angles of attack, 
and unexpectedly well at low angles of attack. This 
low angle of attack correlation would indicate that 
separation occurred at very low angles of attack, and 
the pressure distributions (see Fig. 7) indicate that sepa 
ration first became apparent at the aft station for 9 
a < 12°, and at the forward station for 21° < a < 24 
At first it would appear that these pressure data are 
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not in agreement with the force data since the latter 
suggests separation first occurs at angles of attack well 
below 9°. However, recognizing the fact that the aft 
pressure station is approximately at the midchord 
position, it can be seen that apparently separation oc- 
curred at small angles of attack on the most aft por- 
tions of the wing and that this separation progressed 
forward such that it reached the midchord position at 
9° < a < 12°. Again it should be noted that for the 
round-edge shape, potential flow peak negative pres- 
sures were realized at the transverse section extremities. 
The aspect ratio 2.0 wing data for lift and moment 
have been correlated with the linear and the cross- 
flow theory and are shown in Figs. 8, 9, and 10. For 
the wing with beveled edges, it can be seen that use of 
the cross-flow theory results in a good prediction of 
both the lift and the moment for angles of attack up to 
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about 1S° where the lift and moment slopes begin 
lessen as in a gentle stall. The fact that these correla. 
tions are good indicates that the cross flow separated 
at very small angles of attack, as was expected and as 
is borne out by the pressure distributions. 

The lift and moment data for the round and ellip 
tical-edge shapes are seen to correlate best with linear 
theory, though a small nonlinear effect on the round. 
edge witig is noted at low angles of attack. On the 
basis of these force data alone, it would appear that 
there is little or no cross-flow separation on these wings: 
however, the pressure distributions (see Figs. 1] 


g and 
12) show that separation is first present at angles of 
attack roughly equal to, or in some cases less than the 
angles of attack at which separation became apparent 
on the aspect ratio 1.5 wings having comparable edge 
shapes. 


It appears that the aspect ratio 2.0 wings are 
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more strongly influenced by an additional separation 
phenomenon; that is, the separated flow permanently 
detaches from the upper surface earlier for the aspect 
ratio 2.0 wings than for the aspect ratio 1.5 wings. 
This conclusion is borne out by the beveled-edge wing 
data where it can be seen that the lift and moment data 
first deviate from the cross-flow theory at about a = 
24° for the aspect ratio 1.5 wing, and at about a = 
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IS° for the aspect ratio 2.0 wing. This effect has been 
previously noted!’ and it was proposed that the cross 
flow theory ceased to apply for delta wings at an aspect 
ratio between 2 and 3. For higher aspect ratios the 
experimental data showed that the lift and moment 
fell below the linear predictions, and these “‘critical 
aspect ratios’ were proposed as engineering criteria for 
the area where linear theory gives a better average pre 
diction than the cross-flow theory. The present data 
qualitatively substantiate these criteria but also indi- 
cate that these criteria must also include edge shape as 
well as aspect ratio. These data indicate that for a 
round or elliptical-edge shape, linear theory more 
accurately predicts the lift and moment than does the 
cross-flow theory if the aspect ratio is greater than 
2.0. For a beveled-edge wing, the data show that the 
critical aspect ratio is somewhat greater than 2.0, and 
that for aspect ratio 2.0 wings, the cross-flow theory 
satisfactorily predicts the aerodynamic characteristics 
to reasonably high lifts. 

It is interesting to examine the effects of edge shape 
on the drag of these low aspect ratio delta wings. It 
can be shown that for a thin flat plate model, the re- 
sultant force acts normal to the surface and that the 
induced drag of the wing should be C; tan a. For 
wings with elliptical span loadings at moderate angles 
of attack, classical wing theory indicates that the in- 
duced drag should be C,?/7-R. For low aspect ratio 
wings with elliptical loadings, it can be shown (see, for 
example, reference 28) that the resultant force lies 
halfway between the normal to the air stream and the 
normal to the surface, and the induced drag is C; tan 
(a/2). The drags of the aspect ratio 1.5 and 2.0 delta 
wings are compared with these drag formulations in 
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Figs. 15 and 14. It can be seen that the drag of both 
beveled-edge wings is slightly less than the flat plate 
variation, and it has been concluded that this is a re 
sult of the airfoil geometry. As is shown in Fig. 1, 
the sharp edges on these wings were realized by bevel 
ing the extremities. Due to the separation phenomena, 
there exists a relatively high negative pressure region 
over these edges, and as a consequence, there is a small 
chordwise thrust force. 

The drag of the wings with round and elliptical-edge 
shapes are seen to correlate well with C;, tan(a@ 2), the 
low aspect ratio induced drag for wings with elliptical 
loadings. It can be seen that the aspect ratio 1.5 data 
correlate well with this variation to the lift co- 
efficient where the lift, moment, and pressure data indi- 


up 


cate separation begins, and the drag data then falls 
roughly parallel to the C; tan a curve. The aspect 
ratio 2.0 data are seen to correlate well with the low 
‘spect ratio variation up to lifts somewhat higher than 
would be expected from the force and pressure data. 
Chis would indicate that the low aspect ratio variation 
loes not hold as well for the higher aspect ratio wing 
id that the correlation to high lifts is only fortuitous. 
his conclusion is based on the fact that separation is 
known to exist at the lower lifts and that this separation 
must be manifested in drags higher than those predicted 
by linear theory. 

A rather interesting and quite possibly practical idea 
stemmed from this study of the influence of edge shapes. 
Lawrence has suggested the utilization of edge effects 
to achieve roll control of low aspect ratio wings at high 
gles of attack, since the edge shape influences the non- 
linear part of the lift at those angles. The suggested 


means of achieving rolling moments is to use a side- 


edge “‘improver,”’ so called because it increases the lift 
in the vicinity of the wing side edge where it is used. 
\ simple “‘improver’’ could possibly operate by project- 
ing a sharp plate through a slot on a round side edge, 
thus causing this normally rounded side edge to become 
The added nonlinear 


lift resulting from the sharp edge would induce a rolling 


aerodynamically a sharp edge. 


moment which could be sizable at high angles of attack 
lor wings of low aspect ratio. 
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In order to check this idea, the aspect ratio 2.0 delta 
wing having round edges was fitted with a thin 1 /8-in 
steel plate that extended the entire length of the left 
The 
experimental rolling coefficients C,, referred to body axis 
The theory, in 


side and projected out about 2 in. (see Fig. 15 
are compared with theory in Fig. 16. 
this case, was based on the assumption that the center 
of viscous lift identified with a sharp edge acts at the 
centroid of area of one wing panel, and that identified 
with a round edge acts at the centroid of area of the 
other wing panel. Using this assumption, it can be 
readily shown that the rolling moment arising from the 
difference in viscous lift on the two panels is 


C, = (1/6)a7{(Cp Cp 2] (3 
where the subscripts | and 2 refer to the sharp and 
In the light of the experi- 
without the 


round edges respectively. 


mental results achieved on this wing 


sharp edge, it was assumed that Cp, 2 and Cp 
| 2, and the rolling moment for the wing with one sharp 
edge was 


C, = (1/8)a’ (4 


The experimental verification of the predicted roll 
ing moment coefficient is good, as can be seen from 
Fig. 16. The fact that the experimental rolling moment 
coefficient is higher than the predicted value at low angles 
can be attributed to the fact that, while the viscous lift 
develops on the sharp edge almost immediately upon de 
parture from zero degrees attack angle, the viscous 
lift does not begin to develop on the round edge until 
some finite angle has been reached. The lower rolling 
moment, relative to the predicted curve, at high angles 
of attack can be due to stall effects and/or changing 
cross-flow drag, on the round edge particularly. 

In order to illustrate more clearly the extent of the 
rolling influence realized on this test model there arc 
shown, spotted on Fig. 16, the theoretical deflections 
required of each of the two ailerons on a similar delta 
moment. The 


ailerons are assumed to comprise the area along the 


wing to give an equivalent rolling 


outboard 25 per cent span of the trailing edge having 
The 


aileron deflections were determined by the method re 


a chord equal to 10 per cent of the root chord. 


ported in reference 14 and are conservative in the sense 
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hibited on this curve are required ailavator deflections, 
as determined from experimental results presented by 
Hawes and May,” for a delta wing of aspect ratio 2.51. 
The ailavators in this case comprise the outboard 25 
per cent span of the trailing edge, having a chord equal 
to 12'/5 per cent of the root chord. 

The “‘improver’ might be used further as a substitute 
for a maneuvering flap on low aspect ratio wings at 
high angles of attack. By simultaneously projecting 
sharp plates from both wing side edges, added lift could 
be realized with very little pitching moment change. 


Leading-Edge Influence 

In examining the case of rectangular wings with 
round and sharp edges, it was found!’ that the cross 
flow theory would satisfactorily predict the lift and 
pitching moment on these wings for any side-edge 
shape only if the leading edge were round. For the 
case of a sharp leading edge, it was found that the ex 
perimental lift and moment on the wing was substan 
tially greater than expected, and the indications were 
that additional separation must be taking place. This 
conclusion was borne out by pressure distributions which 
showed that on a sharp leading edge, the flow sepa 
rated at small angles of attack, but subsequently re 
attached to the upper surface at some distance down- 
stream. This point of reattachment depended on angle 
of attack; as the angle of attack increased, the point 
of reattachment moved further aft on the surface. 
On low aspect ratio wings it was found that so long as 
the flow reattached to the airfoil, the lift continued to 
increase nonlinearly, and when this reattachment no 
longer took place, the airfoil appeared to be stalled. 

In comparing the experimental pressure distributions 
over the sharp leading edge with the theoretical pres 
sure distributions of linear theory, it was found that the 
overall effects of this leading-edge separation was to 
eliminate the peak negative pressures near the edge and 
to result in high negative pressures over a greater por 
tion of the wing. As a consequence, the total lift on 
the wing was then increased due to the broadening of 
the negative-pressure region, and the center of the lift 
was moved aft on the wing. Since the cross-flow theory 
refers to linear theory for its basic term, it is desirable 
to visualize these changes due to leading-edge sepa 
ration as additions to the potential-theory pressure dis 
tributions. Consequently, for the present purposes, 
this additional lift is treated as though it stems from 
an additional pressure distribution which can only add 
to and not subtract from the potential-flow pressure 
distribution. Using this concept, the center of the vis 
cous lift appears to act well back on the wing and far 
distant from where the actual physical phenomenon 
takes place. 

In correlating the effects of sharp leading edges with 
the cross-flow theory, it was found empirically in refer 
ence 10 that good lift correlation could be obtained 
if an additional a* term were incorporated in the ex- 
pression for lift coefficient, and good moment correla- 
tion would be obtained if this additional lift were treated 
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as though it acted a distance 6c(dCy dC_)o aft from the 
leading edge. The expressions for the lift and pitching 


moment become 


C; = (dC; da)yp a+ Cp,a* + la? 5 
Ca = (dCy da)va — (& c)Cp,a? + 
O(dCy dC; la 6 


where Cp, is the cross-flow drag coefficient of the sid 
edges. It can be seen that this additional lift term js 
in keeping with the cross-flow concept by noting that 
the cross-flow drag coefficient for an infinite plate js 
2, or the cross-flow drag coefficient per edge should by 
Ls 

In examining sharp-edge delta wings and a certain 
class of sharp-edged tapered wings, it was found that 
the cross-flow theory based on side-edge effects (Cp 
2) agreed well with experiment, and that there was ap 
parently no effect of the leading edge—that is, all edges 
acted as side edges. In considering all tapered wings 
falling between a rectangular wing and a delta wing 
it might be anticipated that there was some leading 
edge sweep angle or taper ratio at which the leading 
edge ceased to incur a separate and distinct separation, 
and began to act as a side edge. An attempt at de 
termining this critical plan form was undertaken by 
correlating the available experimental data with the 
cross-flow theory. Jt was found that there was littl 
experimental data available for the sweep-angle range 
of interest, and consequently a series of tapered wings 
for investigating this phenomenon was incorporated 
in this experimental program. ‘These wings were de 
rived from a basic .R = 1.5 sharp-edged rectangular 
wing by adding delta-shaped noses to the basic plan 
form. Four such wings were tested with leading-edge 
sweep angles of 0°, 10°, 20°, and 30°. Fig. 17 illus 
trates these wings and the pertinent physical character 
istics of each are listed in Table 1. It will be noted 
that all four wings had the same bevel angle normal to 
the leading edge. 

The force measurements obtained from the wind 
tunnel tests were, at first glance, completely contrary 
to the previous results, in that the rectangular wing 
results correlated best with Cp » (no leading-edge 
separation) while the 30° swept wing correlated best 
with Cp, 3 (full leading-edge separation However, 
the pressure data for all wings indicated that leading 
edge separation was present, but did not begin until 
the wing was pitched to an angle of attack of 3° to 6 
It appeared that this might be the source of the discrep 
ancy between theory and experiment since the theor) 
assumes separation to begin immediately at a > U 
Since the bevel angles on these wings were sizable 
half angle of 7.2° normal to the edge) it can be see! 
that leading-edge separation could not begin until the 
local flow at the leading edge was inclined at an angle 
greater than the bevel angle. It should be noted that 
the side-edge separation can begin immediately 10! 
a > 0° since the bevel is normal to stream direction 
In order to take into account the effects of this delayed 
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lge separation, the cross-flow theory was for 


leading-e¢ 


mulated as 


0 dCy dC. a 


where ag is the geometric angle of attack at which 
iding-edge separation first occurs, Cp is the cross 


flow drag coefficient of the side edges, and (2% /c) is the 


wing centroid of area. In order to estimate the sepa 
ration angle of attack, it was assumed that separation 
ccurred when the local angle of attack exceeded the 
bevel angle. Taking the induced angle of attack to be 


rR dc 


then 


da)ya, the separation angle of attack is 


ac TT (] rR (dC, da)yas i 


1 


where ¢ is the half angle of the leading edge measured 
parallel to the air stream. 

The above forms of the cross-flow theory have been 
used in correlating the lift and moment data obtained 
for the four wings, and these correlations are shown in 
Figs. 1S, 19, 20, and 21. Examination of these figures 
shows that good correlation is obtained for all four 
wings and it appears that the leading edge continues to 
ict as such for all leading-edge sweep angles up to 
30° (A 0.7 It should be noted, however, that the 
above equations, when applied to tapered wings, 
assume the side-edge effect to exist over the entire 
wing and the superposition of an additional leading- 
edge effect. If there is some critical sweep angle 
at which there is no effect of a leading edge, it 
would be expected that the experimental data would 
rather suddenly cease to correlate with the above ex- 
pressions and would then correlate with the original 
cross-flow formulations. The fact that this sudden 
lack of correlation was not encountered with these 
four wings could show that if a critical sweep angle or 
taper ratio does exist, it does not exist in the range 
covered by these wings. On the other hand, it is per- 
haps more logical to expect that such a critical param- 
eter does not exist and that there is effectively a gradual 
degeneration of the leading edge into a side edge as the 
sweep angle is increased. At the limiting taper ratio 
of zero, this would then require that the nonlinear 
viscous effects on delta wings be a function of sweep 
ingle or aspect ratio. The data previously discussed 
lor delta wings (see Figs. 3 and 8) might be interpreted 
as showing some slight aspect ratio effect, but at best, 
it can only be said that such an effect is indeed small. 

[tis apparent that the data from these four wings are 
not sufficient for assessing the validity of these two con- 
cepts because of the small variations in viscous effects 
with plan form, and because of the limited range of 
plan forms considered. It would appear that the only 
possible way to determine this question would be to 
test a complete series of tapered wings ranging from a 


rectangular wing to a delta wing. If such a series of 


LOW ASPECT RATIO WINGS 929 





harp-ed 


wings were to be tested, it would be important to en 
sure that the leading edges were similar in order to 
minimize such effects as encountered here—i.e., de 
layed leading-edge separation. 

The lift due to leading-edge separation has often been 
considered only for the case of sharp leading edges. 
Since an aerodynamically round and sharp edge differ 
only in that the separation occurs at higher angles of 
attack and that the effects of the separation are of 
different magnitudes, it would be logical to expect that 
a wing with a rounded leading edge should exhibit char- 
acteristics similar to those of a sharp leading-edge 
wing. Again using the notion of an empirical Cp, for 
the leading edge for thin rounded leading edges, the 
additional viscous term in the lift and moment should 
be proportional to (1/2)a*, while that for a thick 
rounded leading edge should be proportional to (1/4) a’. 
If such an effect for rounded leading edges does exist, 
it might be expected that this effect would not be in 
evidence until some angle of attack corresponding to a 
separation angle of attack for the leading edge. For a 
sufficiently thick leading edge on a relatively high as 
pect ratio wing, this effect might never be reflected in 
the forces due to the counter balancing effect of perma- 
nent-flow detachment on the tip portions of the wing. 
However, the effect might be expected on very low 
aspect ratio wings which appear to suffer least from 
permanent-flow detachment at small angles of attack. 

An effect having all the characteristics of this lead- 
ing-edge phenomena occurring on wings with rounded 
leading edges has been noted in another phase of this 
experimental program, as well as in other data.*! 
However, both sets of data are for tapered wings rather 
than rectangular wings, thereby casting some doubt on 
the possibility of attributing the additional nonlineari- 
ties to leading-edge separation. The effect has also been 
noticed in Scholz’s* data for rectangular wings, but 
since these wings had retreating side edges, this cannot 
be accepted as proof that the effect is entirely due to 
leading-edge separation. It would seem that in order 
to establish the existence of this phenomenon, it will 
be necessary to conduct a systematic investigation on a 
number of rectangular wings with various leading- 


edge shapes. 
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INFLUENCE OF EDGE SHAPE 
Retreating Side-Edge Effects 

For the purposes of this investigation, a wing is said 
o have retreating side edges if there exists a chordwise 
station on the wing where the span is greater than the 
trailing-edge span. It can be seen that according to 
this definition many of the wings used today have 
retreating side edges. 

Since these retreating side edges form a portion of 
the trailing edge, it must be decided whether or not the 
Kutta condition is satisfied along these edges before 
theoretical treatment is possible. Jones*’ treated the 
case of a wing with sharp retreating side edges, and 
in applying the Kutta condition along these edges, found 
the lift to be a function only of the maximum span and 
the lift on those portions aft of the maximum span to 
be zero. More exact theories show that if the Kutta 
condition is satisfied along the retreating side edges, 
the net lift on those same portions will be small and 
may be either positive or negative. 

For rounded retreating side edges, it was concluded 
in reference 10 that the Kutta condition does not appear 
to be satisfied, and that Munk’s* slender-body theory 
might lead to a better approximation for the aerody- 
namic characteristics. On this basis it was concluded 
that the wing lift would be a function only of trailing- 
edge span, and that the net lift on the area outboard 
of the trailing-edge span would be zero. For an esti- 
mation of the pitching moment contribution, slender- 
body theory was used to calculate the centers of pres- 
sure of the wing with tips, and of the wing defined by 
the trailing-edge span. This calculation showed that 
the effect of the retreating side edges was to cause a 
shift in center of pressure forward from that given by 
linear theory by an amount 

x i £ b 
ch | b.. 


dx — l (Q) 


where } is the local wing span and #,... is the trailing- 
edge span. This center-of-pressure location correlated 
with the available experimental data, but since these 
lata were for a particular type of edge shape, it was 
decided to test a series of wings with extreme retreating 
side edges. 

The retreating side-edge wings used in this program 
are illustrated in Fig. 22. Two basic wings of identical 
plan form were tested in this program—one with a thin 
round leading edge and sharp side edges, and one with 
thick round leading and side edges (see Table 1). In 
order to assess the importance of round and sharp re- 
treating side edges, provision was made for the insertion 
of 2-in. wide plates along the retreating edge on the 
round-edge wing. This wing was then tested twice; 
once with the rounded retreating edges and once with 
the sharp plates installed. 

The results of these tests when correlated with theory 
were at first somewhat conflicting in that the data from 
the thin sharp-edge wing agreed well with theory, but 
the data from the thick round-edge model, with and 
without the plates, showed relatively poor agreement. 
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This was surprising in that it was anticipated that the 
only difference in the data from the sharp-edge model 
and the round-edge model with plates would be in the 
value of Cp, in the cross-flow term. In addition it was 
found that the tip correction apparently overestimated 
the tip effects on pitching moment. These results 


lead to a re-examination of the retreating side-edge 
theory as well as an examination into the validity of 
the experimental data. 

In examining the theory it was noted that the tip 
correction always overestimated the destabilizing tip 
effects when correlated with the existing experimental 
data. This tip correction was based on slender-body 
theory which for a rectangular wing predicts all posi 
tive lift to act at the leading edge. Realizing that in 
the actual case the center of pressure is well aft of the 
leading edge, the tip correction was altered to take this 
into account by assuming the negative lift predicted 
by slender-body theory to remain unchanged in magni 
tude and position, and the positive lift to be of the same 


magnitude but acting at the point predicted by lifting 












Retreating side-edge wings 


Fic. 22. 
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surface theory. The lift on the wing was then as orig- 


inally proposed in reference 10 


Cr = (S./S) (dC. /da).a + Cp,a? (10) 


but the pitching moment was altered to become 
Cu = (x/c) (S/S) (dC_/da).a — (#/c)Cp.a? (11 


where (dC, da), is the lift curve slope and SS, is the area 
determined by the trailing-edge span, 7 c is the centroid 
of the entire wing area S, and x c is given by 


b ( b, y dCy | ™ ( b ) 
= — + dx — | (12 
Ps Bes dC, c { De. 


where },, is the maximum wing span and (dCy dC; 
is the linear center of pressure location predicted for the 
wing with sharp edges. 

After examining the experimental data, it was con 
cluded that a sizable portion of the conflicting results 
could be attributed to flow detachment from the blunt 
trailing edge of the thickest wing. It will be noted 
that the trailing-edge angle on the round side-edge 
while the same angle for all 
other wings in this program was only 14 It is then 
logical to expect that the flow would detach from the 
thick wing much sooner than from the other wings. 
The effect of this flow detachment is particularly notice- 


wing was more than 28 


able in this case since for the wing with rounded re 
treating side edges there is negative lift over the trail 
ing-edge area while for the same wing with the sharp 
retreating side-edge plates, there is positive lift over 
the same regions. In order to estimate the effects of 
this flow detachment, it was assumed that 7 S of the 
beveled trailing edge was under flow detachment, re 
sulting in the complete loss in lift on that area. For 
the rounded retreating side edges, this loss in lift was 
calculated in accordance with slender-body theory by 
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using a new trailing-edge span corresponding to tha; 
where detachment was assumed to occur. Such a ca] 
culation for the same wing with sharp plates attache 
is possible but very laborious since only the more exac 
theories predict the positive lift shown by the pressuyr, 
distributions. Since the flow detachment assumptions 
did not warrant such calculations, the loss in_positiy 
lift was estimated by finding the experimental linea; 
slope and comparing this with the theoretical lifting 
surface result. In both cases, the change in pitching 
moment was estimated by applying the lift loss at th 
centroid of the assumed flow-detachment area. 

The correlation between the experimental data and 
the theory with flow detachment correction applied 
are shown in Figs. 23 and 24 for the round side-edg 
wing with and without the sharp plates. Fig. 25 is th 
correlation of the sharp-edge wing data with the cross 
flow theory, and the correlation is seen to be good y 
to high angles of attack. Consequently, it can be con 
cluded that the cross-flow theory holds for the shar 
edge wing and that apparently the Kutta condition is 
satisfied on the sharp retreating side edges. 

The correlations of Figs. 23 and 24 are also seen to 
be good, but owing to the correction for flow detach 
ment, it cannot be said that these data establish the 
exactness of the rounded retreating side-edge correction 
The fact that one round-edge wing configuration (with 
plates attached) differs from the sharp-edge wing only 
along a portion of the side edge which cannot alter the 
linear characteristics certainly lends support to the pos 
sibility of flow detachment. Even if this possibility is ac 
cepted, the calculations to account for this detachment 
can only be considered as crude approximations. _ If one 
wishes to discount the possibility of flow detachment 
occurring on the thicker wing, an indication of the 
effect of rounded retreating side edges and of the worth 
of a tip correction can be realized by comparing Figs 
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INFLUENCE OF EDGE SHAPE 


93 and 25. This comparison shows that by rounding 
the edges, a 30 per cent loss in linear lift was incurred, 
and the wing became, at best, neutrally stable. 

In order to examine the validity of the new tip cor 
rection more accurately, the linear center-of-pressure 
lata and the linear lift curve data of reference 4+ have 
been compared with predictions based on the new tip 
correction and the original lift correction in Fig. 26 
It can be seen that the new correction yields a good 
estimate of the edge effects, especially for aspect ratios 
ibove 1. The lack of more exact correlation at ),. « 

2 can be traced, in part, to the difficulty in meas 
uring the linear moment slope at very small angles 
of attack. It should be noted that in spite of the inac 
curacies of the tip correction, it does serve to predict 
the neutral stability near aspect ratios of 1, while 


if 


linear theory predicts acceptable stability. 


CONCLUDING REMARKS 


It has become apparent that when dealing with the 
erodynamic behavior of low aspect ratio wings, espe 
cially at high angles of attack, the shape of the wing 
edges is a significant parameter. As has been shown, 
the viscous contribution to the lift and pitching mo 
ment may be of a magnitude roughly equivalent to the 
potential contribution or there may be little or almost 
no viscous contribution, depending upon the edge 
shape. 

Further research providing additional data on edge 
influences would be desirable and is to be hoped for. 
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General Air-Frame Dynamics of a 


Guided Missile 


BARNEY L. BYRUM* anp E. RUSSELL GRADY? 


Chance Vought Aircraft, Ine. 


ABSTRACT a 


A set of equations which represents the general air-frame 


dynamics of an intermediate range guided missile is derived and w* 
presented in a form suitable for the application of differential A} 
analyzer and digital computer techniques. The equations pro MV 
vide a somewhat realistic description of the motion of aerodyna - 
mically controlled winged vehicles and are applicable to terminal k 
accuracy studies and dynamic stability studies where higher K 
order effects are of interest. h 
The motion is referred by means of a set of modified Euler . 
angles to a frame of reference fixed on the earth. The modified H 
Euler angles are generalized coordinates that orient the coordi 
nate system whose motion is described. It is shown that equa- Ay 
tions describing the motion of body axes are, in general, more Op! 
useful in the three-dimensional case than those describing the N 
motion of the relative wind system 
It is assumed that the missile is a rigid body of constant mass, V 
that the earth has zero curvature and is not rotating, and that oF » 
the gravity vector is constant in magnitude and direction. If W or V 
the decrease in mass due to fuel consumption is not negligible Y 
during the time of interest or if jet vanes are used in combination 
with the aerodynamic control surfaces, the equations are amen 
able to the necessary modifications 
8 
SYMBOLS 
X,:(X1, Xo, Xs) = principal axes of inertia or p-system, 
orthogonal (sense indicated in Fig. 1) g 
X,:(X1, N, X3’) = skew axes or s-system (sense indicated 
in Fig. 2) ro 
X,:(X1', Xo’, Xs’) = inertial axes or i-system, orthogonal 
(sense indicated in Figs. 1 and 2) b 
Agila, Ye) = relative wind axes or r-system, orthog- S 
onal (sense indicated in Fig. 2) 
Ae; pe) = body axes or b-system, orthogonal (sense 
. indicated in Fig. 1) 
uy, = unit vector along the k-axis q = (1/2)pz,? 
€, = unit vector along the sth skew axis © 
7 = kinetic energy m 
QO, = generalized forces 
o, 0, ¥ = modified Euler angles' or generalized a 
coordinates (q,) A 
’, 0, Vv = time rate of change of the generalized 
coordinates 
o = angle between «and Us (Ug =ahia=u L = (—u,)L 
6 = angle between m and XyXy plane 
(uo = €)5 € = Uy 
y = angle between uy’ and the projection of D ati mgs \D 
u, on the X,-X.-plane (uy = ¢); 
ED — 
. , Y = u,} 
w = angular velocity of the p-system 
w,:(p, q, 1) = components of @ in the p-system 
P.O, & = time rates of change of the w, ; 
Received April 14, 1954. TP, = uz'1 
* Senior Systems Engineer, Aerophysics. Now, Nuclear 
Design Specialist, The Glenn L. Martin Company 2 7 
+t Senior Systems Engineer, Guidance. mg = uz’ mg 


angular velocity of the r-system relatiy¢ 
to the p-system 

angular velocity of the r-system 

the matrix A 


applied moment 
linear momentum 
time rate of change of linear momentum 


angular momentum 


time rate of change of angular momen 
tum 

direction cosine between the axes &, | 

time rate of change of a,, 

intersection of the X2X;3-plane and the 
X,'X2'-plane 

velocity of the missile 

velocity of the wind 


velocity of the missile relative to the air 


angle of attack,? measured from Y’ to 


the normal projection of V, in the 
X,X;-plane 


angle of sideslip,?, measured from the 
normal projection of V, in the X,X;- 


plane to V, 

angle between Y,; and X’ (lies in Y,X 
plane ) 

perpendicular distance from _ missile 
c.g. to thrust axis 

total wing span 

total wing area, except in case of a hinge 
moment derivative when it is that 
portion of the hinged surface lying aft 
of the hinge line 

dynamic pressure 

mean geometric chord 

average mass of the missile, including 
equipment and fuel, during time of 
interest 

aerodynamic reaction, which resolves 
into the conventional aerodynamic 


forces and moments 

lift foree, acts through the c.g. perpen 
dicular to Vy, and parallel to the X,X 
plane 

drag force, acts through the c.g. parallel 
to V and oppositely directed 

sideforce, acts through the c.g. and forms 
a right-handed orthogonal system 
with L and D 


engine thrust, directed along the thrust 
axis that is parallel to but not coinci- 


) 


dent with Y’, Fig. 2 


force due to gravity 
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"gst = total rolling moment about YY’, a body 

axis passing through the c.g. parallel 

to the thrust axis, Fig. 1 

C48! = total pitching moment (excluding that 
due to thrust) about }” (parallel to 
X22), a body axis passing through the 
c.g. perpendicular to the ,Y;-plane, 
Fig. 1 

St = total yawing moment about Z’, a body 
axis that forms a right-handed ortho 
gonal system with XY’, Y’, Fig. 1 

= density of air at flight altitude 
I = time rate of change of moment of inertia 
(1 

J, = generalized momentums 

{4,3 = components in the /-system of the 

time rate of change of } Pe 


(1.0) INTRODUCTION 


A SET OF EQUATIONS which represents the general 
air-frame dynamics of an intermediate range 
guided missile is derived and tabulated herein. The 
equations provide a somewhat realistic description of 
air-frame motion and are applicable to terminal accu- 
racy studies and dynamic stability studies where higher 
order effects are of interest. 

In order to guide a long-range, nonhoming type 
missile to a given point in space, its motion must be 
referred to the primary inertial system—i.e., the mean 
position of the fixed stars--or to a set of axes moving at 
a constant velocity with respect to the primary system. 
The motion of intermediate range missiles may be re- 
ferred to a frame of reference fixed on the earth, how- 
ever, as if the earth were a secondary inertial system. 
Although this is not a logically correct procedure, it can 
be shown that the detriment to terminal accuracy due 
to such an assumption is small (relative to that caused 
by equipment limitations and incomplete theory of 
certain physical phenomena). 

The equation of motion may be referred to inertial 
axes through any set of body axes, Fig. 1, or through 
the relative wind axes, Fig. 2. Just which coordinate 
systems are most suitable for investigating the general 





«ml 





Fic. 1. Illustration of coordinate systems (missile in climb 
attitude). 


motion of a guided missile is still an open question. 
Principal axes of inertia yield the simplest rotational 
equations, while relative wind axes yield the simplest 
translational equations. In the two-dimensional case, 
the present state of theoretical aerodynamics favors 
the relative wind axes, since the aerodynamic relation- 
ships are simpler in this system than in any system of 
body axes. Unfortunately, all moments of inertia and 
products of inertia become functions of the angle of 
attack (a) and the angle of side slip (8) when the motion 
of the relative wind system is described. If a and 8 
are less than approximately 4°, the corrections are so 
small that it is sufficiently accurate to regard the rela- 
tive wind axes as identical with the body axes. Since 
the corrections are not generally negligible during a 
terminal maneuver, however, equations describing the 
motion of body axes are usually more desirable in the 
three-dimensional case. 

The mechanics of a rigid missile in three-dimensional 
space leads to six special, second-order, nonlinear 
differential equations: (a) a set of three equations 
representing the rotational motion of an axis system 
with origin at the center of gravity of the missile and 
(b) a set of three equations representing the transla- 
tional motion of the origin of this axis system. The 
derivation contained in this paper is based on La 


grange’s equations of motion 
(d dt) (OT Og;) — (OT Og } = Q 


The salient assumptions in the derivation are: 

(a) The missile is a rigid body of constant mass. 

(b) The earth has zero curvature and is not rotating. 

(c) The gravity vector is constant in magnitude and 
direction. 

In order to note any mathematical simplifications 
and to facilitate procurement of data on the aerody 
namic relationships, in general, two sets of equations are 
derived: (a) a set of equations describing the motion 
of the principal axes of inertia (and associated missile 
mass distribution) in the secondary inertial system 
and (b) a set of equations describing the motion of the 


relative wind axes (and associated missile mass dis 








Fic. 2 Illustration of coordinate systems 








536 JOURNAL OF THE 


tribution) in the secondary inertial system. The 
motion of each coordinate system is referred by means 
of a set of modified Euler angles to a frame of reference 
fixed on the earth. This fixed frame of reference will 
hereinafter be called the inertial or 7-system. The 
modified Euler angles (¢, 6, W), Fig. 


orient the axes 


2, are generalized 


coordinates that whose motion is 


described in the inertial system. For infinitesimal 

rotations, these generalized angles become the con 

ventional roll, pitch, and yaw angles. 
Prior to launching, the X;,, Xx, 3 


established in the missile by means of gyros 


axes, Fig. 1, are 
and are 
assumed to remain parallel to the corresponding axes 
of the inertial system throughout the time of interest.” 
The error-sensing device most often used in present-day 
aircraft having autopilot control— free or displacement 
gyros—generally gives pickoff indications at arbitrary 
plane attitudes for small deviations from these atti- 
tudes, which are functions of the direction cosines 
relating the body and inertial systems.” 4 

The ever increasing successful application of differ- 
ential analyzers and digital computers as aids in the 
solution of equations describing the behavior of com- 
plicated physical systems,’ the development of three- 
dimensional fluid flow theory, the evolution of the 
inertial guidance system, and the potential capacity 
of the guided missile both as a lethal weapon and a 
scientific tool justify the work presented in this paper. 


(2.0) DYNAMICS OF THE p-SYSTEM 


Equations that describe the motion of the principal 
axis system are derived and presented in a form suit- 
able for the application of simulation and digital com- 
puter facilities. 


(2.1) DERIVATION BASED ON LAGRANGE’S EQUATIONS 


With reference to Table 2 


Zo b £9 
Myf = Apsy (Wsf (1) 
I a t 1 t ) 
psf = Apss 1 ps (=) 


Hencet 
p = © — Wsin 0 
q = Ocos¢@+ W sin ¢ cos 0 
r = —Osin ¢ + V cos ¢ cos 6 
©’ = p+ (qsingd +rcos ¢) tané 
0 = qcos¢?—rsing 
Vv = (qsin @ + r cos ¢) sec 6 

In general, the rotational kinetic energy of a rigid 


body distributed in a given coordinate system may be 


expressed as 


* Under a given set of conditions, analysis will indicate what 
refinements are warranted 
t The singular point (@ = 72/2) does not present any serious 


restrictions in practical applications. (See reference 5.) 
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T = (1/2)2,m,v 


where the origin of the coordinate system coincides wit} 


the center of mass of the body,f and m kth iner 


mental mass, w ngular velocity of the system, an 
r, = vector from origin of the system to ™m 

When the rotational kinetic energy is that of the 
system, 


TX, ut 


and it follows from Eq. (4) that 


With reference to Eqs. (3) and (5), it may be shown that 
the equations describing the rotational motion of th 


p-system are 


0, IT; if) 
Qs = Il. cos @ — IIs sin @ 7 

QO, = —I, sin 6 + Ha. sin @ cos 6 + II; cos ¢ cos 6 
\ 


where 


I, = Pi, + qr(ls - J, 

I], = Ql. + pr, — 71 

IIs; = RI3 + pq — I; 
The generalized forces, Section (5.0), are defined as 
follows: 


Hence 


‘OP = NM, | 

Qe = Mz cos @ — A sin } (10 
Q, = —M, sin 6 + (1\/2 sin @¢ + Ms cos ¢) cos a) 
where 


M, = (C, cos o + C, sin a)gSb 
M. = Cage — Tyr 





M; = (C, cos ¢ — C, sin o)qSb 


Eqs. (6), (7), and (8) show clearly that Lagrange's 
equations essentially express the fact that the applied 
torque equals the time rate of change of angular mo- 
mentum. 

The equations describing the translational motion of 


the p-system are given by 
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t Since a homogeneous gravitational field is assumed in t 
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present derivation, the center of mass coincides with the center 


of gravity 
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TABLE | 
Tabulation of Certain Relationships 


. 2 ; es ae ee 
at Lyaudn = Lava I = Ww 
i ’ l,m,n = 1, 2,3 (cyclically) 
lA~Amszn 
w = (> W sin 0) u, + (8 cos @ + V sin @ cos A) 4+ 
(W cos 8 cos @ 0 sin @)u 
'" - 2 — 
4a _ aw 1 , 
p p l,m,n = X, Y, Z (cyclically) 
lAmsAzAn 
Q= ts Qoyty) (da/dt) + (dp/dt)uz 
“Ss = » * 
Law 7 —_ has Wy 
: ‘ Lmna2wX,Y.Z (cyclically) 
~~ HN =~ Hn 
oa >) ; 
w - hae it 
r 
." es % 
La = 54 a = 0 


t The relationships (a), (b), (ec), and (d) follow from conditions 
f orthogonality 


or, as an alternative, 


izr-l \ at a. \ tos (7 
A iy iiry i Ary 1Qir§ YQirg X 
\ tint es i 
Qing (Rey = Airy VP (12) 
where 
‘ hae ' evige ‘ 
» Lirg Fe = )Qpis ) Ft (13) 


gives the relationship between the forces in the two 


systems 


2.2) PROCEDURE FOR SOLVING THE EQUATIONS 


The equations describing the rotational motion may 
be written in a form more suitable for solution by dif- 
ferential analyzers and digital computers. Since the 
generalized momenta 


J, = OT 0g, = A, + J(Q, + OT Oq,)dt (14) 
Eqs. (6), (7), and (S) may be written 
Jg=Agt+ SQ + 1/2 (Is — I3) (W? sin 26 cos? 6 — 
0? sin 2¢ + 20W cos 6 cos 2¢)] dt (15) 


Jy = Ag + S1Q¢ — 1,0W cos 6 + (1 2)W? sin 26 X 
(I; — I, sin? @ — I; cos? g@) — (1/2)X 
(1. — 1;)0 W sin 2¢ sin 0] dt (16) 


J, = Ay + JQ, dt (17) 
where the A, are constants of integration, and 


J, = 11( — W sin 6 
Jy = (1, cos? @ + I, sin? 6/0 + 
(1°2) (ls — Is) sin 2¢ cos 0 


Jy I, sin? @ + J. sin @ cos” 6 + I; cos* @ cos? 6) VY 
1,® sin @ + (1/2) (Ils — I3)O sin 2¢ cos 6 


The generalized forces are given by Eq. (9). 

Assuming I, = I, (which is a quite reasonable as- 
sumption in many cases), Eqs. (15), (16), and (17 
simplify to a system that can be handled more effi- 
ciently : 


37 
l ; 
@ = I (A, +hwWsin 6+ fQ, dt) (18) 
l 
dO dt = (1/Iz) [Qe + (OT 086) ] 
(1/Iz) [Og — ¥ cos @ (J, + IW sin #)] (19 
— — [Ay + sin 0A, + SO, dt) + SQ, dt] 
I» cos: @ ‘ 
(20 


Beyond this point one may proceed in a variety of 
ways. The procedure suggested is somewhat as fol 
lows. Eqs. (15), (16), and (17) or (18), (19), and (20) 
may be solved simultaneously with Eqs. (25), (26), 
and (27) for a, 8, VV," and the instantaneous missile 
attitude (¢, 6, YW) throughout the time of interest. The 


relationship V, = V’, + IW will then yield the V,” or 
l’,’, and application of (d) or (f) in Table 2 will give the 
lV,’ from which the trajectory of the missile may be 


computed. 


(2.3) COMMENTS ON WIND 


The effects of a changing wind are considered, and a 
means of solving for a, 8, and I’, is indicated. 

In the presence of an arbitrary wind, the relative wind 
vector (]’,) is defined to coincide at every instant with 
the x-axis of the r-system. There is a time delay, of 
course, between the application of a wind on the body 
system and the acquisition of this wind by that system. 

With reference to Table 2, 


Ve} = faut’ | V,' (21) 
tA} = haus’ (Voit + tairf’ (Ao 
or 
dit’ }Aytt = (1/m) } Ff = {A} - 
lait’ }V,'t (22) 
Since 
ay = et Wy 
and 
{Ag} = foul’ {W9 + fan}! {Aoi} 


Eq. (22) may be written 


(1/m) \F,t = $A, + fant’ jAut — 


ot (23) 


Eq. (23) gives the components of the acceleration of the 
missile along the axes of the r-system when the missile 


is in the process of acquiring the wind. 


(2.3.1) Resolution of a, 8, and V, 


Eq. (23) may be utilized, in conjunction with the 
equations describing rotational motion, to determine 
the angles of attack and sideslip, their respective rates, 
and the air speed for utilization in the aerodynamic 
forces and moments. 

With reference to Table 1, Eq. (23) may be written 
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1 (| \4 | a da/dt = (1/V,") [a\(a + o) + As + (pB — g)V? H 
«Fy? = <4 0) + fait’ [Aw + {T} < 0 (24) (1/m)F,] 97 
mI | o lo 
3 where 
where x 2 wo 1] 
41; Le Upiftu ie 
| 0 —w,* w,* | 
Ti = aie 0 —w,* (3.0) DYNAMICS OF THE 7r-SYSTEM 
ey , * E , : 7 ‘ 
~ me 0 ) Equations that describe the motion of the relatiy, I 


F, = mgaz, + T. cos a cos 8 — D wind system are derived, and a number of simplifying 
assumptions that may be valid in a given case are me; 





F, = mgay, — T, sin 8 cosa + Y a 
: ae es tioned. These equations are generally more difficult 
F, = mga3, + T. sina — L . ; : = 
to manage than those resulting when the motion of g The 
body axis system is considered. The variable moments be in 
* ada - - ‘ 
Wy = Ary dt a ph Wy py of inertia, characteristic of the relative wind system ind I 
are discussed in Section (4.0). gener 
dp 
w,* = - » WpLpz ‘ ) 
dt p (3.1) DERIVATION BASED ON LAGRANGE’S Eguations 
, tae ; : : Letting the generalized angles (¢, 6, w) orient th 
Linearizing in and neglecting products of a and 8, Eq. 2 4 . Bi Phe 
4) vield r-system in the inertial system, the equations describing usuall 
(24) yields : ‘ : . i Boo 
i the motion of the r-system are easily formulated. f ine 
At = A, + AB — Axa + c) — (1/m)F, (25) Since the generalized angles now orient the r-system, sideral 
it may be shown that the equations that describe the than I 
i 2 F a _ yz = rz i ‘ i ‘ “ ' ® : 
dp/dt = (1/V,*) [AiB — Az — pV (a + o) —rV, + rotational motion are of the same form as Eqs. (6), (7 sidesli 
(1/m)F,| (26) and (8), where the //, are written as follows: ey ae 
he ass 
a onsid 
ABLE 2 
Coordinate Transformations apa 
: —— : : OE : : = : ; ; proxin 
Suppose that the missile has assumed an arbitrary orientation in the inertial system. The direction cosines relating the coordinat i 
systems employed in this paper may then be derived from unit sphere transformations or methods of projective geometry. The basi such ¢ 
transformation matrices are: nd pi 
(a) {Xp} = {ape} {Xe} (for the components of a contravariant vector sufficie 
where $1den 
j ] 0 sin 6 he 
1dps) = 0 cos @ sin @ cos 6 the 
lo sin o cos 6 cos @ aa F 
(b) [Not ™ UGeel (Xe) for the components of a covariant vector; Section (5.0) lations 
where ems 1S 
l 0 0 
Dye} sin @ tan @ cos @ sin @ sec @ 1() 
cos @ tan @ sin @ cos @sec @ 
(c) iXp} fdpr} {X Let 
with re 
where 
; ments 
cos (a@ + a) cos ps -sin 8 cos(a@ +a) sin(a@ + a) ) i 
lay} = sin 8 cos 8 0 lesired 
sin (a + ¢) cos p sin(a + ¢) sin8 cos(a +a)! uanti! 
(d) {X>} = fap} {X3} associa 
where lunctio 
; ; in dy 
j cos #@ cos ¥ sin y cos 6 —sin @ _ : 
ldpi| = (sin @sin@cos yy — sinycos@ sin@singsiny + cosg@ecos yy sin @ cos #d leserib 
sin 8 cos @ cos ¥ + sin ¢ sin ¥ sin 6 sin Y cos @ — sin @cos yY cos #écos @ | tionshi 
(e {Xp} = {aps} {Xo} } taneou: 
where with re 
j cosa 0 sing rhe 
lan} = ) 0 l 0 pass 
—sing () cose i 
B), fi 
(f) {Xs} = {air} {X-} nition: 
where 


{Girt = fapi}’ {apr 
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GENERAL AIR-FRAME 
H. « (Plz — Qlzy — Ri.2) + 
pr OF sy — 12 a) RHP, = Bie FH te) + 
r(plry, — iy T fie. (2S) 
H PI,, + Ql,, — RI,z) + 
(=f 25: Gly — Flys) + 
b1,. + gly: — rl + r(pl.,. — gi, — ri, (29) 
H PI - OI,. + RI.) + 
= i gV,- + rT + p(—pl, + gi, — 71 + 
g(— pl. + dy + riz: (50 


The p, g, 7, P, Q, Rin Eqs. (28), (29), and (30) are to 
as applying to the 7-system; the I's 


be interpreted 
} and the 


ud I's are to be determined from Eq. (37 


generalized forces are given by 


The maximum values of the products of inertia are 
usually only a small percentage of those of the moments 
of inertia. The maximum value of £ is usually con- 
siderably less than that of a. Since 1; is usually greater 
than I, it is noted that for extremely small angles of 
sideslip, 1, is greater than 1), 1, is approximately equal 
to In, and I, is less than I;.. As in the body axis case, 
the assumption that I, = I. 1s often valid and results in 
onsiderable simplification of the equations. Further 
simplification results when a, 6 remain less than ap 
proximately 4° throughout the time of interest. In 
such cases, the variations in the moments of inertia 
ud products of inertia are so small that it is usually 
sufficiently accurate to regard the relative wind axes 
isidentical with the body axes. 

The equations describing the translational motion 
(11), where A,, k,, 

and F,. The re 


lationship between the applied forces in the two sys 


of the y-system are given by Ea 
ind F, are to be replaced by A,, k 


» 


tems is given by Ea. (13). 


1.0) TRANSFORMATION OF MOMENTS OF INERTIA 


Let (/) denote an orthogonal system that is rotating 
with respect to a rigid body. Suppose that the mo 
ments of inertia and products of inertia of the body are 
desired in the rotating system. In such a case these 
quantities become variables, since mass increments 
issociated with ‘points’ in the moving system are 
lunctions of time. This is the situation encountered 
in dynamic stability studies where the equations 
describe the motion of the relative wind system. Rela- 
tionships are required which are functions of the instan 
taneous direction cosines orienting the rotating system 
with respect to an axis system fixed in the body. 

The moment of inertia of the body about any axis 

passing through the origin of an orthogonal system 
B), fixed relative to the body, follows from the defi- 
nition: 


bo 


1, = > m(r.xu, (* 
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where 
“uy = Yapitty unit vector along | 
B 
r. = NX, uy (vector from origin of the B-system 
B 


to m,.) 
X zp: (&, 9, £) (axes of the B-system 


It may be shown that 


I; = > linda dpe 
iB 

where 

1,B=8é,7y, 3 

, jl;A=B - 
e = i) 
(-1;A #B 

and I4zn = Iga moments of inertia when A B and 


4 me 


products of inertia when A # B. Thus, in general, 
the moment of inertia of a rigid body about an arbi 
trary axis is a quadratic function that contains both 
the moments of inertia and the products of inertia as 
expressed in a coordinate system wherein the body 1s 
fixed. 

Expressions for the moments of inertia about the 
axes of the system are found by repeated applica 
tion of Eq. (33), but the determination of the products 
of inertia in this system requires an additional equa 
tion. Let A, uw, v denote the axes of (/) and let I,, de 
note the products of inertia in the system; then by 


definition 


In terms of the B-system, 


& * > aBX p! 
B 


T; S a, BX p’ 


ra 


B 
and Eq. (54) becomes 
|. = oF m,.0.<4X 4° A,BXp' A,B =é,9,¢ 
1B 
= = >. lp pp; — I.,(d:.4,; + ¢.,.G T 
B 
Tar(Qplee + Bre Apr) H ler (Cee Ape TH Age Ae, 309 


If the axes of the B-system are principal axes of 
inertia, the products of inertia vanish in the body 


system, and Eqs. (33) and (35) simplify to 


where 
é,7T =A, HW, 1 
e j= lje=r 


30 


By repeated application of Eq. (56), the following 
relationships result for the moments of inertia and 


products of inertia in the 7-system: 
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l,-, = 11ai;° + Todo," + 1303,7 
Ly = Tidiy? + Teae,? + 9s,’ 
(.. = liai.? + Iede,? + 03," 
= —(1,Q1;01, + Iod2,d2, + 1303,03,) 
— (1,a),41, + Ied2,02, + 131332) 
= —(1a).04, + Tede.do, + 1303-03, 


(5.0) TRANSFORMATION OF MOMENTS 


The transformation matrix }6,,{ is derived, and it is 
shown that the components of the applied moment 
transform in this manner. The conventional tensor 
notation is employed. 


(5.1) DERIVATION OF | Onsf 


Let Y denote an arbitrary vector, then 


Xa X% = X,¢ a, B, y = p-system (38) 
and 
. =e = 2 1, j, R = s-system (339) 
where 
Xo ete (contravariant transformation) (40) 
From Eqs. (38) and (39), 
X“(€,°€") Y,(€° -€ 
X%’ =XG” i,7 = «¢«-¢ = fl;a = 7 
lO;aF# ¥ 
x? = X¥G” 
or 
a at # 
and 


Similarly, 


and 
X; = GiyX! 
Now, substituting in Eq. (40), 


1aB ve ’ , 
G™X loos Aaja" X j 


or, in matrix notation, 


(GY I Xet = (daily (GY [XG (41) 
and, since, 
[G™i—! = {G")’ 
SX pf = pai (G7 5 (G4 Xi (42) 
where }a,;| = }a,,{, |G! is a unity matrix, and |G} 


is to be determined. 


Since 
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{ | Q —sin 4) 
[Gi} = 0 0 
Pipi 6 O | ) 
A ' 
G” = : G cos- @ 
G 
and thus 
\ sec” 6 Q sin @ sec | 
‘G4 = 0 0) 
bcin 6sec?@ O sec” @ | 
Hence 
Dass = Vat (GY G4} = Ldpsf 13 
and therefore 
IXot = pout Xe (covariant transformation 


(5.2) FORMULATION OF THE GENERALIZED FORCES 


The generalized forces 1.e., the O,--are defined as 


the components of the applied moment in the s-system 
The manner in which the components of the applied 
moment transform is desired. 


By definition 
AM rm I] 


With ry and F defined in the p-system, .\/ has a covariant 


description, 


M = M,& = M,eé’ = M’ «, M,=M 
AM,.( eX e;) : 
= M*é,. (t,. 4, k; evclic 45 
Vv IG; 
With the following substitutions, 
VG cos 8, } €: 1 ae Ce 


and collecting coefficients of the e,, .1/ is written 


M = Ae + e(sin @ tan 0 1, + cos @M2 + 
sin @ sec 0.1/;) + (cos ¢ tan 0.1/, 
cos @ sec 6.1/3) = Me, (46 


sin @.\/. + 


which yields the transformation 


where } ./7,{ = | QO}. 
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Propagation of Stall in a Compressor 


Blade Row 


FRANK E. MARBLE# 
California Institute of Technology 


ABSTRACT 


par 


observations on compressors, in 


Recent experimental 
ticular those of Rannie and Iura, have clarified some features of 
the phenomenon of stall propagation. Using these observations 
isa guide, the process of stall in an airfoil cascade has been char- 
acterized by a static pressure loss across the cascade which in 
creases discontinuously at the stall angle, the turning angle being 
from this simple 


flected in only a minor Deductions 


model yield the essential features of stall propagation such as 


way 


dependence of the extent of stalled region upon operating condi 
tions, the pressure loss associated with stall, and the angular 
velocity of stall propagation. Using two-dimensional approxi 
mation for a stationary or rotating blade row, free from interfer 
ence of adjacent blade rows, extent of the stalled region, the total 
pressure loss and stall propagation speed are discussed in detail 
for a general cascade characteristic. Employing these results, 
the effect of stall propagation upon the performance of a single 
ixial compressor is illustrated and the mechanism of enter- 
The essential 


stage 
ing the regime of stall propagation is discussed 


points of the results seem to agree with experimental evidence 


INTRODUCTION 


[' THE RATE AT WHICH FLUID passes through a com- 
pressor blade row is decreased while the rotative 
speed remains constant, the angles of attack of the com- 
pressor blades eventually exceed that for which un- 
separated flow is possible. Individual blades then 
behave in a manner somewhat resembling stalled iso- 


The 


nature of compressor operation in the stall regime is of 


lated airfoils and the compressor is said to stall. 


considerable importance not only because of the attend- 
ing loss in efficiency and pressure rise, but also due to 
the fact that the onset of stall is accompanied by 
pulsation and oscillating blade loads which may lead 
to structural failure of the The complexity 
and general irregularity of the flow during stalled oper- 
ation have, in the past, discouraged a detailed study. 
have 


blades. 


investigations 
demonstrated the 
ippears that the blades of a given blade row do not 
stall simultaneously, but initially only a few blades 
stall. Furthermore, it is observed that the stalled re- 
g10n ol 
same group of blades but that the region moves about 
the blade row at an appreciable fraction of the rotor 


Recently several experimental 


some regularities in process. It 


a blade row does not always consist of the 


The name stall propagation or rotating stall 
In addition to the rea- 


velocit y. 
Originates from this behavior. 
sons mentioned above, stall propagation is of impor- 
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tance through its possible connection with compressor 
surge and because the appearance of the more or less 
regular phenomenon may indicate that the compresser 
is operating near its surge line. 

During the past few years experimental observations 
at several research institutions and manufacturing 
establishments have indicated stall propagation in a 
variety of compressors. Unfortunately, most of this 
information exists only in private communications. 
It is probable, however, that the regularity was first 
recognized by Emmons together with his co-workers at 
Harvard University, and some of their findings have 
recently been reported in reference 1. Of the avail- 
able experimental results concerning stall propagation, 
the most detailed and complete are those of Rannie 
and Iura,* performed at the California Institute of 
Technology on a large three-stage axial compressor. 
Knowledge of these experiments has been the principal 
guide in the present formulation of the problem. 

The results of Rannie and Iura indicate that two 
types of stall propagation occur which are, at least 
superficially, different; in one the stall covers the entire 
blade length, in the other only a portion of the blade 
length. As_ the throttled, a _ single 
stalled area is observed which covers only the root or 
With further throttling 


ce AN presst rious 


tip regions of a few blades. 
two, three, or more similar regions appear, uniformly 
spaced about the periphery. The propagation speed re- 
mains remarkably constant 
Finally, with further decrease in flow, the pattern 
a single region extending over the blade 


as throttling proceeds. 


changes to 
length; this change is accompanied by a small decrease 
in propagating speed. This type of pattern persists 
until the flow is nearly shut off, the peripheral extent 
of the stalled region increasing as the flow is reduced. 
Some information on experiences with other machines 
has been reported by Huppert and Benser.’* 

The basic features of the stall propagation phe- 
nomenon seem to be quite clear. When the blade row 
operates at nearly its stall angle, any small disturb- 
ance will cause one or more of the blades to stall. The 
stall reduces the blade force and decreases the local 
pressure rise across the blade row so that the flow will 
be retarded near the stalled region and diverted to 
either side. Depending upon the geometry of the 
blade row, the flow angles upstream of the blade row 
will be increased on one side and decreased on the other. 
Where the flow angles increase, neighboring blades 
stall and the stalled region propagates in that direc- 


tion. Where flow angles are decreased, blades already 
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stalled may be relieved so that the stalled area retains 
a definite length. The rate of propagation depends 
upon some characteristic time associated with the stall- 
ing process. This characteristic time is related to the 
process of boundary-layer separation and establish- 
“stalled” flow pattern from an unstalled 
and 


ment of a 
flow pattern. 
Iura clarify this point. 
moval of every second blade from one or more of the 


further results of Rannie 


They observed that the re- 


Some 


blade rows had only a slight effect upon the speed of 
stall propagation. Consequently the influence of time 
required for boundary-layer separation must be rather 
small, since the propagating speed would be inversely 
proportional to the blade number if this were the govern 
ing phenomenon. Therefore the speed of propagation 
must be associated with the time required for the exist- 
ing pressure difference across the blade row to trans- 
form the flow pattern. This time is not strongly de- 
pendent upon the blade geometry or the blade gap. 
On this basis an elementary model for the stalling proc 
ess was constructed. Then it was not a difficult matter 
to deduce the resulting details of stall propagation. 

It is a pleasure for the author to express his deep 
gratitude to Professor H. S. Tsien who supplied the im- 
petus to start the investigation and the encouragement 
to finish it. The initial formulation of the problem 
and the requirements of the model for the stall process 
were suggested by Dr. Tsien. Professor W. D. Rannie 
likewise contributed greatly to the present picture of 
stall propagation through illuminating discussion of 
his experimental work. It was Dr. Rannie who first 
pointed out that considerable clarity could be gained 
by departing from the conventional formulation ac- 
cording to airfoil theory of cascades. Professor W. R. 
Sears, Cornell University, has been most helpful; first 
in communicating the results of his analysis prior to 
publication, and later through stimulating discussion. 


THE PuHysIcAL MECHANISM OF STALL PROPAGATION 


An airfoil exhibits two modes of operation: the 
normal mode, many features of which are well described 
by the theory of ideal fluids, and the stalled mode which 
may be described roughly by the ideal flow about a 
boundary consisting of a portion of the airfoil and a 
constant pressure stream line that bounds the stagnant 
wake. The pressure associated with the stream line 
is, in this case, equal to the undisturbed static pressure 
of the free stream. When airfoils are arranged in a 
cascade or lattice, as they are in a compressor blade 
row, they possess similar normal and stalled modes of 
flow. 
duce the flow velocity from upstream to downstream so 
that for the normal mode of operation a static pressure 
Thus when one air- 


In a compressor cascade, the blades act to re 


rise occurs across the blade row. 
foil of a cascade is stalled for some reason, the situation 
is somewhat different from that of an isolated airfoil 
because now the associated with the free 
stream line is that far downstream of the cascade which 
Since the 


pressure 


is considerably in excess of that upstream. 
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high pressure now penetrates to a point near the lead. 
ing edge of the stalled airfoil, the pressure in the entir; 
vicinity is increased. As a consequence the velocity 
approaching the channel formed by the stalled surfac 
and the next blade is reduced, as shown in Fig. 1, and 
the quantity of fluid passing through the channel js 
proportionately decreased. 

Now if, during normal operation of the entire cascade. 
one of the airfoils stalls, the flow pattern corresponding 
to stall develops through intermediate stages indicated 
in Fig. 2. As the separation point moves forward 
from the trailing edge, the wake begins to develop 
and spreads downstream with a velocity less than the 
local free-stream velocity. The growth of the stall 
region gradually reduces the flow through the passage 
by diverting the fluid ahead of the cascade. The 
time elapsed during transition from unstalled to stalled 
velocity fields consists of that needed to separate the 
boundary layer and that required for the existing pres- 
sure difference across the cascade to deform the flow 
pattern. The latter is a simple inertial effect, while 
the boundary-layer separation is more complicated; 
actually the processes are not strictly independent. 

The mechanism of stall propagation may be seen by 
observing, from Fig. 2, that, during successive stages 
in development of stall for one airfoil, the angle of 


When 


this airfoil reaches its critical angle of attack, it too 


attack on the airfoil just above it is increased. 
stalls and the process continues as shown in Fig. 3. 
Just how soon the adjacent airfoil begins to stall de- 
pends upon the normal angle of attack and the cascade 
geometry as well as upon the rate with which the 
field is modified by the airfoil 
Furthermore, the angle of attack is 
below the stalled one 


neighboring flow 
originally stalled. 
decreased upon the airfoils 
Therefore as the stall region grows upward, airfoils far 
below will eventually become unstalled and the process 
will settle down to a steady propagating stalled region 
of constant length. The actual length depends, in 
particular, upon the proximity of the normal angle of 
attack to the stall angle. 

The time required for the stall of one blade to deform 
the flow sufficiently to stall the next airfoil is a most 
significant item, and the factors that affect it may be 
isolated by an approximate analysis of dynamic equilib 
rium of the flow while the airfoil is stalling. Suppose 
that during normal operation a static pressure rise 4) 
takes place across the blade row. If the pressure rise 
is small, it is related to the change Aw in velocity from 
inlet to outlet by Ap = Now when the 


blade stalls, the fluid passes through the stalled blade 
There- 


— pwdw. 


channel with negligible rise in static pressure. 
fore, assuming the static pressure field downstream of 
the blade row to be only slightly affected by the stall, 
a pressure difference Ap is available for retarding the 
from far upstream to flow 


fluid which approaches 
Since the blade 


through the stalled blade channel. 
chord is neglected, the only characteristic length im the 
problem is the gap g (or length of stalled region), thus 
the pressure gradient that retards the flow is Ap/Ké; 
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where x is a constant of proportionality. At the 
initial instant of stall, where the approach velocity w 
does not vary with distance upstream of the airfoil, 
the rate of change of w is given by Ow/Of = —(1/p) 
(Ap/xg). An estimate of the time 7 required for the 
stalled pattern to develop may be obtained by assum- 
ing the transformation to continue at this rate. Now 
since there is no pressure rise through the blade row in 
stalled conditions, the flow approaching the blade 
row must undergo a change in velocity almost equal 
to Aw in order that it may discharge to the pressure 
downstream of the blade row. Therefore, approxi- 
mating Ow/Ot by Aw/r, it follows that Aw/r = 
—(1/p)(Ap/xg) so that employing the relations be- 
tween Ap and Aw, 


WT - = x 


If the rotative speed of the compressor is QR and m is 


the number of blades, then 


wr/g = [w (QR)] [(2QRmr)/(mg) | 


where w) (QR) is usually about unity and for the present 
consideration is constant. But mg = 27R and the 
angular velocity of stall rotation is w = 2m /mr, since 
mr is the time required to stall the total number of 
blades. Thus approximately wr/g = Q/w, the ratio 
of blade speed to that of stall propagation. Then 
since wr/g = « the ratio of stall propagation and rotor 
speeds is constant, w 2 = I «, and in particular is 
independent of the blade gap as indicated in experi- 
ment. 

It is usual that several blades are involved in ob- 
served stall propagation so that the blade spacing is 
not the significant dimension but rather the length of 
the stalled region. Consequently it is appropriate 
to simplify the model to an actuating line across which 
the velocities and pressure may change discontinuously. 
It is convenient to think of this simplification carried 
out by allowing the chord and gap to vanish simul- 
taneously while their ratio remains fixed. Then it 
remains only to assign realistic conditions on the 
change of flow properties across the line. 

To follow the convention found convenient in com- 
pressor practice, let the local inlet flow angle 6, be the 
independent variable, Fig. 4, in describing the char- 
acteristics of the cascade. The performance of the 
cascade is completely defined by prescribing (1) the 
local turning angle or the outlet flow angle and (2) 
the local static pressure rise across the cascade. The 
discharge angle @) is a well-determined function of the 
inlet angle, and, in the absence of losses, the static 
pressure rise is then fixed also. For considering the 
flow with respect to the blade row in question and 
accounting for the fact that the velocity Ul’ normal 
to the cascade is equal ahead and behind the cascade, 
the Bernoulli relation gives simply Ap = (1/2)pl” 
(tan? 8, — tan*® 8). By virtue of the dependence of 
82 upon @,, the pressure rise across the cascade is deter- 
mined. Now when a relation between Ap and 8, 
different from this is assigned, it can be possible only 
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because the total pressure changes. Consequently 
choosing Ap(§;) is equivalent to assigning a value of 
the total pressure loss at each value of the inlet angle 

During stall propagation the entire blade row js 
operating sufficiently close to the angle of stall that 
the cascade characteristic need be approximated near 
this value; let 8,* be the inlet flow angle for which 
the cascade stalls. Locally the turning angle is 
linear function of the inlet angle as shown in Fig. 5, and 
there is experimental evidence to show that it is not 
profoundly influenced by the stall. Then it may be 
characterized by a discharge angle 8,* and a slope of 
magnitude a just prior to stall. The pressure rise 
may likewise be approximated by a linear function of 
the inlet angle below the stall angle, but it changes dis- 
continuously at 8,*. According to the discussion of 
the stalling mechanism, the static pressure rise across 
the cascade may be approximated reasonably well by 
assuming it to vanish above the stall angle, to have a 
value Ap* at 6,*, and to have a slope 6 just prior to 
stall. This approximation to the pressure rise char- 
acteristic is shown in Fig. 6. 

Such characterization of the cascade performance 
and stalling process differs from that employed by 
Emmons or by Sears in either of his analyses.*° The 
differences are not of great importance so far as the 
unstalled cascade characteristic is concerned, for any 
pair of performance parameters which depend upon 
the upstream flow angle is adequate, although not 
equally convenient. For example, the performance 
of a cascade may be expressed by giving the lift co 
efficient C, and the drag coefficient Cy, as functions 
of the angle of attack. Properly these coefficients are 
based upon the geometric mean of the upstream and 
downstream velocity vectors and the angle of attack 
measured with respect to its direction 6,,. Then if 
F, and F, are the blade force components in the x and 
y directions, respectively, referring to Fig. 7, the lift 


and drag coefficients are just 


C, = (—F, sin B,, + F, cos Bn) 
[(c 2)pl*(1 + tan’ 8, 


Cp = (F, cos B,, + F, sin Bp) 
[(c/2)pU7(1 + tan? Bn)] 


and the angle of attack is simply 
1 
Bn — ay = tan [(tan 8B, + tan ~)/2] 


where a is an arbitrary angle of zero lift. An elemen- 
tary momentum calculation shows that F, = pl’g(! 
tan 8, — U tan 8), while the axial force results from 
pressure difference and is F, = —Ap g, where 4 
is the static pressure rise across the blade row. Now 
if C, and Cp are specified as functions of the angle of 
attack, then clearly the expressions for F, and F, 
given above determine 8. and Ap as functions of 8). 
Therefore this viewpoint is equivalent to that used 
in the present analysis. Using the blade channel 
approach, the outlet angle 6. and the total pressure 
loss Ap, are specified as functions of the inlet angle 41. 
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This is closely related to the specification employed in 


the present analysis, for clearly 
Ap = (1/2)pl/*(tan? 8, — tan? B,) Ap, 


and hence Ap(,) and 6.‘3,) may be obtained directly. 

It is interesting to compare the mechanisms of stall 
prescribed by Emmons! and that suggested in the 
present analysis. Emmons utilizes the fact that stall 
restricts the flow through a stalled blade channel and 
recognizes that the development of steady stalled flow 
from the initial separated flow requires a time that is 
the governing factor in the rate of stall propagation. 
However this time is not calculated from the details 
of the nonsteady motion, and consequently no stall 
propagation speed may be determined. Emmons and 
his co-workers seemingly overlook the fact that the 
characteristic time is proportional to the peripheral 
extent of the stalled region and therefore draw the 
erroneous conclusion that stall propagation speed is 
proportional to the wave length. It must be recog- 
nized, however, that Emmons’ mechanism contains 
elements that, if pursued to their logical conclusion, 
would produce propagation speeds essentially in agree- 
ment with those of the present analysis. The fact that 
Emmons et al. do not account for the influence of flow 
downstream of the blade row is a significant but not a 
decisive shortcoming. The ‘amplitude’ of stall and 
hence its peripheral extent and associated pressure 
loss are essentially nonlinear effects that, because his 
mechanism is not employed to full advantage, Emmons 
does not determine. 

In his initial treatment, Sears! employs the approach 
of nonstationary airfoil theory, assuming an irregu- 
larity in lift near the stall angle but neglecting losses 
associated with stall. He assumes, moreover, that 
the actual stall lags a certain arbitrary phase behind 
the occurrence of critical inlet angle, this lag pre 
sumably being associated with the time required for 
boundary-layer separation, satisfying the Kutta con 
dition or some other real fluid phenomena. It would 
seem that the most serious defect of this mechanism is 
its failure to account for the sharp increase in drag 
associated with stall. At the suggestion of Rannie, 
Sears improved his cascade characteristic in reference 
5 to account for drag or pressure loss at the stall. He 
employed the conventional “‘channel theory’ with the 
reasonable restriction that variation of discharge angle 
with inlet flow angle is neglected; the concept of stall 
phase lag, however, was retained. Since Sears solves 
the nonstationary flow problem for both cascade pre 
scriptions, the velocity of stall propagation is deter 
mined except for influence of the arbitrary phase lag. 
Under special circumstances described later, the stall 
propagation speed computed under this channel theory 
agrees with the result developed herein. In both 
cases, however, Sears linearizes the cascade perform- 
ance in the neighborhood of stall and consequently 
is unable to determine values for the extent of stall 
region or the pressure loss associated with stall. 
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GENERAL THEORY OF STALL PROPAGATION IN SINGLE 
BLADE Rows 


Consider a cascade formed by an annular clement 
of radius R from a rotating blade row. During normal] 
operation the flow relative to the cascade is represented 
(Fig. +) by aconstant component UL’ of axial velocity and 
constant flow angles 8; upstream of the cascade and 8, 
downstream. Corresponding uniform values of static 
pressure Pp, and p, exist ahead of and behind the 
blade row, and, if 8; is close to the stall angle 3;*, the 


pressures satisfy the relation 


po — py = Ap = Ap* + b(B; — Bi*) (1/2)pU’ 
During stall propagation deviations from the uniform 
state occur which depend upon both location and time 
The problem may be reduced to one of steady-stat. 
flow, however, by referring to a coordinate system x, J 
moving with the velocity wR of the stalled region. Let 
V be the undisturbed tangential velocity component 
in this system; the quantities u, v are the velocity 
perturbations and p the pressure perturbation asso 
ciated with stall propagation. A subscript | will de 
note conditions ahead of the blade row and a subscript 
2, those behind. 

The correspondence between inlet and outlet flow 
angles necessitates that perturbations A@; and Aj, of 
the inlet and outlet flow angles satisfy the relation 


AB> = aAp, | 


so long as the inlet angle is near the stalling angle. In 
prescribing the static pressure jump across the cascade, 
it is necessary to distinguish between stalled and un- 
stalled regions. Denoting the compressor circumfer- 
ence by 27R on the y axis, the’stall will be assumed to 
cover a fraction a of the circumference extending over 
the interval —ar < y R < az. This pattern is re- 
peated with a period 27 so that, as indicated in Fig. , 
the actuating line is stalled in the regions (27 — a)r < 
y/R < (2n + a@)aw and unstalled elsewhere. There- 
fore the perturbation pressures /,(0, y) and p,(0, ¥ 
immediately ahead of and behind the cascade, satis!) 


the condition 


p20, vy) — pi(O, v 


(1/2) pl? 
Ap* y 
~ oe b(B; = Bi*)s (2n —a)r< < 
(1 2) pl . R : 
(2n  @ 4 
bAB; : elsewhere 


Due to the fact that the undisturbed flow fields ahead 
of and behind the cascade are initially uniform, the 


pressure perturbations satisfy the Laplace equation 


(0°p,/Ox*) + (0°p,/Ov") = 
(O7p2/Ox?) + (O°? p2/Ov" 0) 


Furthermore, the pressure perturbation due to the 
stall must vanish far upstream and assume, at most, 


a constant value far downstream. Then for a pre 
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SINGLE scribed state of steady operation the cascade perform- function but that it is indeed the harmonic conjugate of 
ance given by Eqs. (1) and (2) determines completely P and, further, that this useful relationship holds so 
the pressure field for the corresponding state of stall long as the vorticity in the field is small. 

lement propagation. Now it is required further that the inlet Far upstream of the cascade both pressure and angle 

normal flow angle exceed 8,* in the stalled regions and be less perturbations vanish and hence Eq. (5) may be inte 

=~ than 8;* in the unstalled regions. These conditions grated directly to 

7 mav be met only when the stall propagation speed 1 . 

-— Pr wR and the length 2a7R of the stalled region are prop- 2[(m/U) + (/U) tan i] = —P, 

Static erly chosen. Thus if there is a condition of stall Noting that 

7 the propagation which corresponds to a given set of inlet - " - 

> the conditions, the peripheral extent of the stall and the Opi/On = —(1/2)pU%(06,/00 
rate of stall propagation are determined. Eq. (6) also integrates in a similar manner to 

pl” The flow angle deviations that are involved in the oie } i 

~ boundary values of the problem must be found from 2[(%/U) — (m/U) tan 3,] = 6 

form the pressure field. When the undisturbed flow angle is The local flow velocities u; and 7; are then simply 

ume # = tan! (I’/U), then the disturbed flow angle 3 + Pa a ; . 

a Ad is related to the axial and tangential velocity per- RE al} = i 

a, J turbations as 2(1 + tan? 3;) (v,/U) = —P, tan’, + 6 (12) 
Let 

onent tan (# + Ad) = (V + 2)/(U + wu) In order to reduce the problem to one of steady mo- 

locity and since « and v are small tion, the reference frame has been taken fixed in the 

asso rotating stall pattern, but when dealing with inlet 

Il de- (1 + tan? &)Ad’ = (v/u) — (u/U) tan 8 (4 and outlet angles of the blades it is necessary to con- 

script These velocity disturbances are related to the pressure sider angles that pal wae with — ~~ Oe 

blades themselves. Since the variation in flow angle 


distribution through the equations of motion which, 
flow for the present purpose, are most conveniently written 
MP2 Of along and normal to the streamlines of the undisturbed (1 + tan? B)AB = (v/U (u/U) tan B 


with respect to the blade is just 


| flow. Let the coordinate along the streamline be o 


. ae . it is convenient to write in general 
and that normal to the streamline, v, as indicated in : : 








(] 
Fig. 9. The equation of equilibrium along the stream- (1 + tan? B)AB = 
In line is (1 + tan? d)Ad — (wR/U) (u/U) (12 
cade, U(d/d | »tan @) = —20/d 5 : f 
1 un- oi wey ee) = oer (9) where account has been taken of the fact that 
nfer- and has an obvious integral which is the familiar Ven 8 - Utena é «= off 
ed to linearized Bernoulli relation. The equilibrium equa- 
over tion normal to the streamline is where w is the angular velocity of stall propagation. 
S re- , , In particular, the variation of inlet angle to the blade 
| pU(0/0c) (v — u tan 3d) = —Op/or (6 pe ; " 
io. 8 row is, using Eq. (10), 
r< and, by reference to Eq. (4), is clearly related to the ' P 
a ee ‘ . . ae 2 sec 8; sec J = )) s — dv t 
1ere- angle deviation Ad. In fact, Eq. (6) is just SOe fe SHE PLA Ae) a — %1) 
) 6,(0) cos (8; — 3%) (13) 
, 9 ea ‘ 
ea pU*(1 + tan? 3) (0/d0a) (Ad) = —Op/dv 
ust) Now the condition for the pressure jump given by 
Furthermore, the continuity relation for this coordinate Eq. (2) may be written explicitly as 
system says simply that 
P(0) — P,(0) = 
(0/Oc) (u + vtan 3) = —(0/dr) (v — uw tan #) [  Ap* 7 
. = + b(8; — 6,* 
and consequently Eq. (5) reads (1/2) pU? 4 
2 pl*(1 + tan? #) (0/drv) (Ad) = Op /da (Qn — a)e < y/R< (Qun+ a)r 
o : (14) 
hen two dependent variables may be defined 
. b r wR 
ead P=hp [(1/2) pU?] (4 , cos® B, |} (1 + LD sin 3; cos 3))6,(0) + 
he ee - L 
the 6 = 2?(1 + tan’ 3)Ad \R ] 
w ’ « 
; ; cos? 3 P,(0) |; elsewhere 
which reduce the equations of motion to J 
3 OP Oc = 06 /Ov (S) Downstream of the cascade the flow may be treated 
. inasimilar manner. In addition to satisfying the above 
the OP /Ov = —086/00 (9 ae ‘ ‘ , , 
condition on pressure jump, the discharge angle vari- 
st, Chese are the familiar Cauchy-Riemann relations which ation A@. must bear the relation to the inlet angle 


re- State not only that @ (and hence Av’) is an harmonic prescribed by the cascade characteristic, Eq. (1 
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Then utilizing Eq. (12) and the fact that the axial ve 
locity is continuous across the cascade, Eq. (1) may be 


written 


cos? 82) 00) + (wR UL) cos? 3;[P1(0) + 6,(0) tan 3) ] 
= ( COSs- By a0 + wR cos a, [P, 0) + 
6,(0) tan & ]{ 15 


Summarizing the problem, the pressure and angle 
perturbations vanish far upstream of the blade row 
while the pressure perturbation p.( and the angle 
perturbation Ad.(©) may at most have constant 
values far downstream. Taken together with Eqs. 
(14) and (15), which respectively express the physical 
conditions that the pressure jump and the cascade 
discharge angle are those prescribed by the cascade per 


formance, the problem is completely specified. 


STALL PROPAGATION WITH SIMPLE CASCADE 
CHARACTERISTICS 


The details of stall propagation are particularly 
simple when the cascade discharge angle is unaffected 
by the inlet angle and when the pressure rise across 
the cascade is zero or Ap* accordingly as the region 
in question is stalled or unstalled. In high solidity 
cascades this condition may be nearly realized; in the 
present formulation it corresponds to the case when 
a = b6 = 0. Then the pressure jump condition is 
simply 
P.(0) — P,(0) = 

Ap* ae 
| = eis Ce < | 1 
(1/2) pU? R (16) 
| 0 : elsewhere 


Since P; and Ps» are both harmonic functions, it 1s 
convenient to commence construction of the solution 
by obtaining the harmonic function that assumes 
the values P.(0 P,(O) along the y-axis, is regular 
elsewhere, and approaches at most finite constant values 
far ahead and far behind the blade row. This is 
easily done by recalling the ideal flow field associated 
with an infinite row of sources or sinks. The com 
plex potential of a unit source at the origin is just 
[1 (27)] log s where s = x + (ty is an arbitrary point 
in the plane. The imaginary part of this general po 
tential, the stream function, represents the flux of 
fluid across a curve connecting the point s with any 
fixed point of the plane. If s moves upward along the 
y-axis at small negative values of x, the flux increases 
by half the source strength as 2 passes the origin, the 
change being discontinuous for sufficiently small 
values of x; and thus /m(1,/7) log z possesses a unit 
jump at the origin. Then clearly the desired har 
monic function may be constructed by a superposition 
of such jumps; more precisely it may be constructed 
by superposing a negative unit jump at each point 
y = 1(2n — a)r and a positive unit jump at each point 
y = 1(2n + a)r. The positive jumps, considered sepa- 
rately, correspond to an infinite vertical row of sources 
with spacing 27R, and the resulting complex potential 
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is known (see, for example, Lamb*) to be 
1 mw) log ;sinh (1 2 sR) + tar}! 


The potential representing the negative unit jumps is 
entirely similar so that the desired harmonic function js 


jsinh (1/2) [(z R lat | | 
F(x, vy) =" Im log 2 
; T (sinh (12) [(s/R) + iar]! 


) 


( ytan(1 2 v R) — ar]| 
tan : 

1 ' tanh (1/2) (x/R / 

jytan (1 2) [(vy R) + aw)) 


tan 
| tanh (1/2) (x/R { / 


Phe function F(x, vy) is odd in x; in fact 


: jl) 
F(O-. y at 
while 
: pom 
*(() 
, : 1 of 


and similarly 


Now the perturbation pressures P; and P» can differ 
from a multiple of F only by an harmonic function 
that is regular except at certain points of the cascade 
Its form is determined by the remaining condition t 
be imposed, Eq. (15), on the cascade discharge angle 
From the fact that P and @ are harmonic conjugates, it 
is clear that the supplementary harmonic function be 
ing sought can differ from the conjugate of F by, at 
most, a constant. Denote this function G(x, y) such 
that 


OG Ov. = OF Ov, OG Ov = = OF On 


Then the pressure field may be written 


P, = A[F(x, y) — a] + BG(x, v | 
[ Ap* 
P, 4| [Fix, v) + a] - 
1/2 pl J . 1s 
ap* | 
BG(x, y a : 
. 1 2)pl 


and the field of flow angle perturbations as 


4, = —AG(x, ¥) TF B|F s, 3h = a | | 
Ap* 
o = — Y— —~AB 4 9 
(1/2) pU? : | 


B[F(x, vy) + a] + aC 


where A, B, and C are constants to be determined. 
Now the angular velocity of stall propagation w and 
the fraction a of the compressor circumference which 
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‘< stalled for a given inlet angle are determined by the 


erturbation of inlet angle which is induced by the stall 
+celf. Referring to Eq. (13) and substituting the 
values for P; and 6;, the inlet angle perturbation may be 


writtel 
sec sec AG, A cos (8 v - 
B sin (3 v,)1G(0O, v) +4 
{ sI 3 7} + B cos (8 - wo F () _ a) 
20 
t s1mce 
Ee 4sin (1/2) [(2/R lat | | 
G(x, 4 R/ log ( 
r (sin (1 2 ¢/R) + tar 


t is clear that G(O, vy) and hence, in general, 49,(0, y 
have logarithmic singularities at the edges of the stalled 
region, v R =aTr Therefore, the only possibility 
f satisfving the condition that 8; = 3,* at these points 
is that the coefficient of G(O, v) vanish in the expression 


for the inlet angle perturbation. Therefore 


A cos (3 v = Bsin (8; — Vv) 


A(1 + tan 8, tan 3,) = B(tan B, — tan 3) (21 


uid employing this to simplify Eq. (20), the expression 
for the inlet angle perturbation becomes simply 


A|F(0 a = a] 


2(tan 8, — tan vd, 

A a 
F(O \ - i 22 

2 wR l 
The determination of constants A and B requires 
imposing the condition on cascade discharge angle 
given by Eq. (15). In the present example where 
0, Ag 0 and consequently 

(UY) + wR U’) cos 3, [P, 0 + 6 Q) tan 3; | = 0 (23 


Substitution of the proposed solutions [Eqs. (18) and 
19)] gives, due to the linear independence of F(O™, v 
nd G(0, v 


wR l A : 2 — tan- vy —_ 
1 + tan 8, tan d,)|JB = 0 (24 


2(1 + tan? d, (1 + tan 8, tan 3;)]A + 
wR/U)B| = | Ap*/[(1/2)pU?]} 
(1 + tan? d,; 25 
+ tan? a))C = [(wR/U)A] — 
2(1 + tan? 3,) — (1 4+ tan 8, tan 3;)] (26 


It appears immediately from Eq. (26) that the down- 
| ; | 

Stream flow angle Avd.(« is unperturbed by the 

stall. Comparison of Eqs. (21) and (25) gives directly 
that 

é ‘ oun A 9° 

A = (1/2) {Ap*/[(1/2)pl?]} 27 

There remain two homogeneous relations, Eqs. (21 

and (24), to be satisfied and this may be accomplished 


only if their coefficients are related—that is, if there 
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exists a relation between the inlet angle and the speed 
of stall propagation. Dividing these two relations gives 


simply 

wR U = 1+ tan-d 
or in terms of the inlet angle relative to the blades, 
wR/U = (1/2 1 + tan? 8) tan 8,] ese 23, 28 


As indicated while discussing the physical mechanism 
of stall, the propagation speed is independent of the 
peripheral extent of the stalled region in conformity with 
the experiments. Since the critical item in determining 
the propagation speed is the inlet angle to the cascade, it 
seems reasonable that the propagation speed is inde 
pendent of the stream deflection in passing through the 
blade row. It is not surprising that this is identical 
with the result obtained by Sears,’ when the phase 
lag associated with his assumed pressure loss vanishes 

The extent of the stalled region must be determined 
through consideration of the actual distribution of inlet 


flow angle which is now known explicitly to be 


AB, = |1/[4(wR/U)]} [F(O-, vy) — @] 
fap*/[(1 2)pl?]} (29 


The total angle of attack 8, + Ag; must equal 3,* at 
the edges of the stalled region, and, since F(O~, y) is 
discontinuous at these points, the proper limiting value 
is its mean, which is equal to 1 2. Therefore the con 
dition A8,(0, + axrR B,* — @, determines the frac 


tion of stalled blade region to be 


50 


So far as the dependence of a upon 0; is concerned, it is 
predominantly the factor 3,* — 8, which affects the 
result, the variation of csc 23; being small under usual 
circumstances. It is clear that the stall will first appear 
for angles below that which would cause the cascade to 
stall uniformly. On the other hand, the inlet angle 8 
must exceed 3,* by a certain amount before the entiré 
blade row is stalled. It must be kept in mind, hew 
ever, that the solution becomes inaccurate long before 
this condition is reached since the linearization rests 
on the assumption that a < 1. 

The propagation speed wR UL is shown in Fig. 10 as a 
function of the inlet angle. It 1s of interest to note that 
the propagation speed is particularly insensitive to 
variations of inlet flow angle in the neighborhood of 
8, = 2m 4, a value that is usual for normal axial com- 
pressors. This is certainly related to the observed fact 
that the propagation speed does not vary markedly as 
the compressor is throttled. The manner in which the 
stalled fraction a increases with inlet angle is shown in 
Fig. 11 for the particular values of 3,* = 7/4 and 
Ap* /[(1/2)pU?] = 0.5. It appears that a rather 
narrow range of inlet angle is involved in passing from 
the first appearance of stall propagation to a com 
pletely stalled blade row. While this is probably real- 
istic near the onset of stall, it is certainly in error as 
the stalled region grows. The pressure rise Ap* has 
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a strong influence upon the angle at which stall first 
appears as well as the range of angles over which stall 
propagation is possible. The angle 6, for incipient stall 
propagation is decreased for a heavily loaded blade 
row. Physically this says simply that the ‘‘relief’’ on 
other blades due to a region of stalled blades is greater 
when the stall involves a larger loss of pressure. To 
illustrate this influence, the variation of a@ with £; is 
shown in Fig. 12 for values of Ap*/[(1/2)pU?] = 0.5, 
1.0, and 2.0 with the inlet angle 8,* at which stall occurs 
equal to 45°. 


STALL PROPAGATION WITH GENERAL CASCADE 
CHARACTERISTICS 


The results of the previous section are somewhat 
modified when the pressure rise across the cascade and 
the discharge angle from the cascade depend upon the 
inlet angle near the stall. The pressure field ahead of 
the blade row is of the same form as before 


P,(x, y) = A[F(x, y) — a] + BG(x, y) (31) 


but the downstream pressure field must be modified in 
order to satisfy the general condition on pressure jump 
given by Eq. (2). In the stalled region the pressure 
drop is now | Ap*/[(1/2)pU?]} + 6(8, — B:*) so that 
the pressure field downstream of the cascade could be 


chosen as 


P.(x, Vy) = (| Ap*, [1 2) pU?]! + b(p, = 6,*) —_ A) x 
[F(x, vy) + a] + BG(x, y) — 
a(} Ap* [(1 2)pU?}} + b(B, = 8,*)) 


were it not for the fact that the pressure rise across the 
cascade is proportional to the inlet angle perturbation 
in the installed region. Now since the formal expres- 
sion for the upstream pressure field is identical with 
that employed in the solution for simple cascade char- 
acteristics, the expression for Af; is general and hence 


P,(0, vy) — P,(0, vy) = (6/2) [A/(wR/U)] (1 — @) 


in the unstalled region. The appropriate expression 


for P.(x, y) is therefore 
Ap* 


P(x, vy) = = 
1 (L/2)e0* 


+ b(B, = B1*) —_ A| x 
A 


) 
—- (l—a) xX 
2 (wR/U) . 


[F(x, y) + a] + BG(x, y) 
Ap* 


b(B — B,*) (32) 
(1/2)pu2* ” 


[F(x, vy) + 1] —a | 
inasmuch as F(0*, y) + 1 vanishes in the stalled region 
and is equal to unity elsewhere. The corresponding 
flow angle variations are then 


6, = —AG(x, y) + B[F(x, y) — a] (3: 


j | = s 
o= ’ 
: (1/2) pU? 


b A 
( 
2 (wR/U) 


b(B; — B,*) — A| G(x, vy) + 
B[F(x,y) + a] + 


| i a)G(x, y) +aC (34 


where, as before, C is a constant to be determined. 
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The condition on the outlet angle, given by Eq. (15 
now gives, upon substitution of the explicit forms fo; 
the pressure and angle perturbation fields, three rel, 
tions to determine the constants A, B, and C. 


(tan 8; — tan #,) [1 — a(cos? B;/cos” B.'|A 
,2(1 + tan? 3,) — [1 — a(cos? B, cos? g 
(1 + tan J; tan B,){B = 0 (33 


Ne, b (1 — a) , cos? 2 
a4 ( + tan‘ d,) — — ¢ - 
\ * 2 (wR U) _ aa ] ‘. 


; { 
(1 + tan ?; tan Bi) A + (tan 6; — tan d,) X 


cos” B Ap” 
(: = ) B = (1 + tan? d)) r 
cos~ Be (1/2 of 
b(B; — a) | Sin) 


C = (tan 2, t. 


] 

|. 
i. 2 
»2(1 — tan? 3) — [1 — a(cos? B,/cos? Bo 


(1 + tan #, tan 6; iB 0 (37 


tan #,) [1 — a(cos? B,/cos? po 


where as before the constant C = 0 so that the flow 
angle far downstream is unperturbed by the stall 
Furthermore, it follows by comparison of Eq. (21) with 
Eq. (36) that the constant A is just 


| Ap* 
‘ P - + b(B, — B,*) 
2 (L(1/2)pU? 
Eaten 
(3S 
1 (wR/U 
which reduces to the previous result [Eq. (27)] when 


b= @ 
done for the simple cascade characteristic, the corre- 


Now by dividing Eqs. (35) and (21), as was 
sponding result is obtained for the speed of stall propa- 
gation 
(wR/U)? = 14+ tan? dj + 
2a cos” B,/(cos? By — cos? B;)] X 
(1 + tan #, tan J 

Upon elimination of #, the relative inflow angle, the 
propagation speed becomes simply 

wR/U = [1 + a(cos? B,/cos? B2)] ese 8; ay 
Since a is, in general, a positive quantity, it appears 
that the influence of discharge angle deviation is to im 
crease the propagation speed over that which would be 
obtained with a fixed outlet angle, a = 0. 


This value of the constant A, which has been ob 
tained in Eq. (38), allows computation of the inlet 


angle perturbation according to Eq. (22 Thus 
Ap* 
~ + b(6:——:*) : : 
AB; = (1/2) pU? [F(O-, vy) — a} (40 


(wR/U)+(b/4) (1—a@) 


and using the same reasoning as that employed in the 
problem of the simple characteristic, the fraction ©! 
the blade row periphery which is stalled is 
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l m (1/2\pU? 
a= $1(3,* — By) : 
Z Ap* 
| , Cc 1 2B, b } 
= @ = ese 28, + (41) 
L cos” Bo Ss 
where terms of higher order in 6;* — 8, have been de- 


leted. Eqs. (39) and (41) give stall propagation speed 
wR/U and the peripheral extent of stall a as functions 
of 8; alone for fixed values of a, b, and Ap*/[(1/2) pl]. 

Since the value of a is usually small, the variation 
of outlet angle is not a strong influence on stall propa- 
gation speed. Since cos* 6,/cos? 8 is usually less 
than unity, it is not probable that the propagation 
speed will be augmented by more than 25 per cent in 
this manner. The corresponding small influence upon 
the peripheral extent of the stall region is to decrease 
the range of inlet angles for which stall propagation 
may be expected. Similarly the inlet angle where stall 
propagation will first occur is increased by an increase 
in the cascade discharge deviation, a. The most in- 
fluential factor in determining the extent of stall regions 
for a given inlet angle is the cascade pressure rise Ap*, 
as was found to be true with the simple cascade char- 
acteristic. The influence of the slope of pressure 
characteristic at the stall is illustrated in Fig. 13 where 
the variation of a with @; is shown for the same values 
of Ap*(1/2)pU" used in Fig. 12, with the exception 
that now the slope is taken equal to 6.0. A large slope 
of the pressure characteristic, such as that shown, may 
exert considerable influence in reducing the range of in- 
let angles for which stall propagation is observed. 


PERFORMANCE OF A STALLED COMPRESSOR STAGE 


The results of the analysis of stalling performance of 
a cascade may be employed to estimate the performance 
of a simple compressor stage, Fig. 14, which consists 
of a rotor following an inlet guide vane, which will be 
considered infinitely far upstream. Then referring to 
Fig. 15, the absolute outlet angle from the guide vane is 
y. It is convenient to discuss the compressor per- 
formance in terms of the two conventional parameters, 
¢, the flow coefficient, and y, the pressure coefficient. 
These are defined as 


@ = U/(QR) (49 
Y = Ap,/[(1/2) p(QR)?] 1! 


where @ is the angular velocity of the rotor and A 
is the mean increase in total pressure across the stag 


Then the compressor characteristic is given by plotting 


¥ versus ¢ as indicated in Fig. 16 where the perform 


ance is shown for normal operation up to the stall point 


that would be reached if all blades stalled simultane 
ously. 

Now according to the velocity diagram the inlet angk 
3, to the rotor is 


tan 6, = (1/¢) — tan y (44 


so that the value ¢* of the flow coefficient for which the 
normal stall occurs is @¢* = (tan 8, + tan y)~. Fyr. 
thermore, the total pressure rise Ap, across the rotor is, 


for no losses 1n the blade row, 


Ap, = pQR(U tan 8B, U tan Bs) 


or in terms of the pressure and flow coefficients 


y = 2¢ (tan 8, — tan B 


which holds up to the stall angle. Thus specification 
of the flow coefficient and pressure coefficient is equiva- 
lent to specifying the inlet and discharge angles from 
the rotor. 

Referring to the compressor performance of Fig. 16, 
it is clear that, as the flow coefficient @ is reduced, a 
point is reached where stall propagation is first pos- 
sible. This is just the condition that a 0 and, 
assuming for the moment the simplest cascade char- 
acteristic where a = 6 = 0, it follows from Eq. (30) that 


Ap*/[(1/2)pU?] = 8(6,;* — B) ese 28, 


Now according to Eq. (44) 6:* — 8; (cos* B,/¢° 
(¢ — ¢$*), and, consequently, the value ¢; of the flow 


coefficient for which stall propagation first appears 1s 


Ap* 


(46 
| (1/2) pl 


%: = ¢*|1+ - (1 — ¢* tan y) 


Actually, of course, the stall can not start with a re 
gion of arbitrarily small peripheral extent but must 


involve at least one blade channel. Under usual cir- 


cumstances this represents a negligible influence on 


The total pressure loss involved due to the propa- 


gating stalled region is simply equal to the static pres 
sure perturbation far downstream of the blade row, the 
change in kinetic energy resulting from distortion of the 
flow field being of second order. Thus, according to 
Eq. (18), the pressure field perturbation far downstream 
is P3(o) = —a} Ap*/[(1/2)pU?]} and hence the changé 
Ay in pressure coefficient due to stall is just 

Ay ~ —a¢*| Ap/[(1/2)pU?]} i 
The pressure loss is thus simply proportional to the 
fraction a of the circumference which is stalled, so that 
utilizing Eq. (30) of the previous calculation 
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Ay = [—2¢*/(1 — $* tan y)] (¢ — @) (48 


Therefore, as the compressor stage is throttled, stall 
propagation first becomes possible at a flow coefficient 
¢,, as the flow coefficient is reduced farther, the pressure 
coefficient falls below that which would occur in the 
ibsence of stall by an amount proportional to the re- 
duction of the flow coefficient below ¢;. To the pres- 
ent approximation this reduction continues until the 
1, at which 


periphery 1s entirely stalled—that is, a = 


point the loss in pressure coefficient 1s simply 
r2}t 
pl*); 


\ctually the small perturbation analysis employed is 


¢?{ Ap*/[(1/2) 


invalid when a becomes of appreciable size, and losses 
computed under these circumstances are at best of the 
proper magnitude. 

For the general cascade characteristic the flow co 


efficient for incipient propagating stall is 


. jf 1 ‘ Ap* 
¢,=o* {1+ (1 @”* tan y) = a 
(4 (1/2) pU? | 
cos” B,* 9 sie 
l+a rp b sin? 6,* (49) 
cos? Bo* 32 f 


while the change in pressure coefficient Ay due to stall 


becomes, neglecting terms of small order, 


2o* cos” B,* 

Ay * l+a a > _ 
1 — ¢” tan 7 cos* Bo 
b 


(50) 
2(1 + tan? B,*) 
the principle effect of variations in pressure rise and 
cascade discharge angle due to changing inlet angle is, 
under normal circumstances, to increase the pressure 
loss somewhat. 

To illustrate the influence of propagating stall upon 
the compressor performance, consider the arbitrary 
compressor characteristic given in Fig. 16. Under con- 
ditions of uniform flow the stage stalls at a rotor inlet 
angle of 8; = 8,* = 51° 30’ and at a corresponding flow 
coefficient ¢ = ¢* = 0.56. Assume further that the 


guide vane discharge angle is y = 28°, while a = 0.25 
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Fic. 18. Rate of stall propagation for compressor stage 


mately those values for the guide vane and first rotor 
near the root of the compressor employed in the ex 
periments of Rannie and Tura using vortex blading. As 
the compressor is throttled, the first possibility of stall 
propagation occurs, according to Eq. (49), at @ 0.6 
Then as the flow is further reduced, the pressure co 
efficient drops below the normal characteristic by an 
The 


resulting performance with propagating stall is shown 


amount that is easily calculated from Eq. (50). 


as a broken line in Fig. 17. 

Now suppose, as is usually the case, the flow through 
the compressor is reduced by adjusting a throttle down 
stream of the compressor while the compressor speed 
is held constant. Then the pressure rise Ap, imposed 
upon the compressor is proportional to U*, or in terms 
where the 
There 


of the compressor parameters y Ko’, 
coefficient x depends upon the throttle setting. 
fore the states of normal operation and propagating 
stall which correspond to the same throttle setting lie 
“throttling” that 

It is probable, there being some 


along a_ parabolic curve such as 
indicated in Fig. 17. 
indication from the 
that the stall propagation does not appear as soon as 


the flow coefficient ¢; is encountered, but, rather, the 


experiments of Rannie and Iura, 


transition from normal operation to stall propagation 
takes place as the result of instability of the former at a 
flow coefficient somewhat less than ¢;.__ In this case the 
transition must be to that state of stall propagation 
connected by the appropriate throttling curve. 

The values of stall propagation speed for the com 
pressor stage follow directly from Eq. (39) using the 
relation between flow coefficient and rotor inlet angle 
given by Eq. (44). Since the speeds of propagation 
are usually observed relative to the fixed blade rows, 
the significant quantity is (2 — w)R, and the ratio of 
absolute angular velocity of stall propagation to the 


rotor angular velocity is just 
1 — @(wR/U 


This speed of propagation for the compressor stage 1s 
shown in Fig. 18 at the performance conditions of Fig. 
17. The stall propagation speed observed by Rannie 
and Tura is about 0.3, and it is of interest that the value 


of propagation speed computed near the onset of stall 
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that is, at ¢ = 0.6 —is in remarkable agreement with 
their value. In view of the facts that a two-dimen- 
sional theory has been employed and that the inter- 
ference of neither the guide vanes nor the stalled blade 
rows downstream of the first rotor is accounted for, 
the agreement must be considered, in a certain measure, 
fortuitous. Although the magnitude is _ essentially 
correct, the only discrepancy seems to be the varia- 
tion of calculated propagation speed with flow co- 
efficient, while the experimental value appears con- 
stant over a wide range. This feature, however, 
is simply a result of the linearization. For whereas in 
the linearized analysis the flow approaches the stalling 
blades at nearly the undisturbed upstream flow angle, 
the real flow approaching the blades is strongly dis- 
turbed by stall and seems to adjust itself to almost a 
constant flow angle, regardless of how severely the com- 
pressor is throttled. Consequently, as was clear at the 
outset, accuracy of the calculation can be expected 
only near the condition of incipient stall. 


CONCLUDING REMARKS 


From the experimental data that are available at the 
present time, it appears that the model upon which 
the foregoing analysis is based is a realistic one and 
contains the essence of the phenomena that have been 
observed. Based upon this model, the features of 
stall propagation that is, the extent of stall region 
the pressure loss associated with stalled operation to- 
gether with its attending effect upon performance, as 


well as the speed of stall propagation arise naturally 
and simply. Discrepancies in comparison with ex 
perimental results on propagation speed cannot b; 
interpreted as inadequacy of the model until either g 
more detailed experimental investigation of the mecha 
nism is undertaken or the results are extended t 
situations more nearly approximating those on which 
observations have been made. Extensions appear 
logical in two directions: the effect of mutual in 
terference between adjacent blade rows in relatiy; 
motion and stall propagation in a three-dimensional 


actuator desk. 
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An Application of Variation Methods to the 
Computation of the Deformation of a 
Cantilever Plate Subjected to Nonuniform 
and Aeroelastically Induced Loads 


EUGENE E. COVERT* 
Massachusetts Institute of Technology 


SUMMARY 


The equations that govern the shape of the deflection surface 
f.acantilever plate subjected to an arbitrary transverse load are 


solved in terms of an influence surface which is approximated by 


the Rayleigh-Ritz procedure. If the loading is also dependent 


ipon the deflection, this procedure results in an integral equation 
The deformation surfaces of rectangular plates with one edge 
lamped and the remaining edges free, for both uniform thick 


ness distribution and double wedge thickness distribution have 


been calculated and are compared with measurements. The 


increment in the slope of the lift coefficient caused by the deforma 
tion was calculated for a wing of aspect ratio unity flying at a 


Mach Number of 7 2 


SYMBOLS 


= coefficients in series expansion 
= semi-span 


= semi-chord 


D = plate modulus = &f*/[12(1 vy?) 

E = Young’s modulus 

h = semi-thickness distribution 

C,, Cc. = Busemann’s pressure-distribution coefficients 

Cee = slope of the lift curve 

,£n = mode function 
K = kernel function = p(x, yo) (Ow, /Oy 
k = constant of proportionality resulting from solution 
of Eq. (18) 

V = Mach Number 

n = direction normal to plate edge 

5 = direction along plate edge 
= plate thickness at maximum thickness point 
= dynamic pressure = (1/2) pV? 

At = pressure difference across surface 

P, p = loading distributions 

( = constant part of loading 


= deformation dependent part of loading 
V = flight velocity 
= plate deformation due to a unit point load 
= plate deformation due to loading p(x, y) 
plate deformation in > direction 
approximate solution 
\ = x coordinate of unit load 
y coordinate of unit load 
a aerodynamic angle of attack 
plate energy 
6 variation operator 
v = (0°/Ox?) + (0?/dOy?), Laplacian operator 
ir density 


v Poisson's ratio 


Received August 2, 1954 


* Aeronautical Engineer, Naval Supersonic Laboratory 


INTRODUCTION 


iy THE DESIGN of high-speed aircraft and missiles, the 
wing and tail surfaces are frequently of a high 
solidity ratio and can be analyzed in terms of thin plate 
theory. This analysis is used by the designer to com 
pute the stresses in the surfaces and by the aerody 
namacist to compute the increment in the aerodynamic 
characteristics due to the elastic deformation. In this 
application, the plate is subjected to mixed boundary 
conditions (clamped at one edge and free at the re 
maining edges) which are such that exact solutions of 
the governing equations are rarely found. The most 
common engineering approach is the application of 
the minimum energy theorems and the calculus of vari 
ations, and here several alternate methods of getting the 
approximate solution are practical. Fung'? and 
Luce* * have used the Ritz procedure to find approxi 
mate solutions for cases with uniform lateral loading, 
pure torsion, and uniform shear located along the edge 
opposite the clamped edge. Reissner and Stein® have 
applied a semidirect method that consists of substi 
tuting a set of approximating functions which are known 
in one direction and undetermined in the other direc 
tion into the energy integral. The first variation of the 
energy integral leads to a family of linear differential 
equations whose solution defines the latter set of ap 
proximating functions. This approach has been used 
successfully for torsional and uniform lateral load 
distributions.* ° 

The solutions to the plate bending that have been 
cited are all for fairly simple loading systems. While 
the methods are applicable to problems where the load 
ing is nonuniform, the procedure is cumbersome, and 
due to its nature, the solution cannot be expressed 
algebraically fora completely arbitrary load distribution 

This difficulty is eliminated herein by computing the 
plate deformation due to a unit load using the same 
variational methods that the previous investigations 
employed. Once this solution is known, the problem of 
computing the deformation for an arbitrary nonuni 
form loading reduces to one of integration of the product 
of the load distribution and the unit load solution over 
the plate. This approach is then generalized to include 
the case where the load is not only nonuniform, but may 
also depend linearly on the deflection or its derivatives. 
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This latter problem arises in estimating the aeroelastic 
influences upon aerodynamic characteristics of flexible 


structures. 
ANALYSIS 


Assumptions 


The analysis that follows contains within it the 


usual assumptions of thin plate theory namely 

(1) the thickness of the plate 1s small compared to the 
other dimensions of the plate, and 

(2) the deformation of the plate is small compared 
to the thickness of the plate. 

The first assumption implies that both the stress in 
the Z direction (see Fig. | for the coordinate systein 
and the shear stresses about the X and Y axes are so 
much smaller than both the stresses in the X and Y 
direction and the shear about the Z axis, that the former 
stresses may be neglected. The second assumption 
implies that the only stresses are those due to the bend 
ing. Since the deflection surface is anticlastic, it 
would be expected to set up tension in the midplane 
and thereby increase its rigidity, but the small defor- 
mation assumption limits the deflections so that these 
tensile forces are negligible; therefore it also follows 
that the midplane is unstretched. The assumptions 
allow for the effects of two-dimensional stress—i.e., 
the deflections are reduced by the lateral extension. 
The latter statement is, in some cases, equivalent 
to replacing the Young's modulus (/£) by //(1 — v?). 


Formulation of the Problem 

Since the governing differential equation is linear, a 
solution may be constructed by superposition of suit 
able primary solutions.’ In other words, if a known 
solution for a unit point loading at the point x», vo is 
W(X, VY, Xo, Yo), then the solution for an arbitrary load- 


ing [P(x, y)]is 
a hed 
“u(x, y) = } 704(X, V, Xo, Yo)P(xo, Vo)dxodyo (1) 
If the loading is dependent upon the deformation in 
some manner, which may be written 
P = u(x, y) + v(x, y)f(w) 


where f(z) is either w or any one of its derivatives, then 


the solution 1s 
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._,; 
w(x, y) = } W(X, V, Xo, Vo)U(Xo, Vo)dxod yo 
Jd S 
. . 
W(X, V, Xo, Vo) V(X, Vo)f(w)dxydy 
Jd S ; 


Further, if the first integral on the right-hand sid 


is bounded, and if w(x, ¥, Xo Yo), V(X0, Vo) and f(w) ar 


bounded over the region of integration, it can be shown 
that the integral equation may be solved by repeated 


substicution and the series so obtained is convergent 
and converges to the solution.’ It will be noted that 
if f(w) depends upon w and its derivatives, Eq. (2) be 
comes a fairly complicated integro-differential equation 

The problem now becomes that of defining the fun 
tion w(x, y, Xo, V This will be done by application oj 


) 


the minimum energy principle! * 7 


a eee 
67 6 | | h*(y) ) (Vex “+ 2(1 v) X 


9) 


' 
H a 


O'w O°w || P 
- —— | dedy — | wP(x, y 
exe" Gy" a ; PS : 


, Ow g I n 
| M, dS — | Ww (0. _ mi ‘ as| 0 (3 
. On Os | 


dxdy + 


S 


where /(y) is the normalized semithickness distribution 
[h(y) = h(—y)], W,, Q,, and J/,, are respectively the 
bending moment, the shear, and the torsional moment 
at the edges, and the last two integrals on the right 
hand side of the equation are line integrals about the 
edge of the plate.’ If the first variation vanishes and 
if the following inequalities are satisfied, 


0<r<07 


and if the deflection is such that each spanwise line of 
the deflection surface intersects the X-Y plane only at 
the root, then the necessary and sufficient conditions 
for a “weak minimum’” are satisfied. A weak mini 
mum is adequate, for the deformation and its slope are 
small, and the loadings are such that the first derivative 


of the deformation surface is continuous. 


Approximate Method of Solution 


There are at least four methods of approximately 
satisfying the minimum energy requirements |Eq 
(3) ]; however, the Ritz method will be used herein and 
the sequence of approximating functions (see Appendix 
will be chosen in the manner prescribed by Fung’ 
i.e., let the solution be of the form 
wx, vy) = wol(x, vy) + DD AmnWm 


haw hae 


where 

Winn = fm(x)gn(y) 
and f,,(x) and g,(y) are the complete set of mode shapes 
of a vibrating beam having the proper thickness dis- 
tribution and the proper boundary conditions. Eq. 
(6) is substituted into Eq. (3) and the resulting equa- 
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AN APPLICATION OF 


tious for the A,,, are found [using éw = (Ow OA») X 
5A4,,,], such that the integral, Eq. (5 


)) is chosen to be the beam bending 


is zero. The 
trial solution (% 
ind torsion solution for a load distribution typical of 
the class of loadings being studied with / replaced by 
E/(1 y This will reduce the correction terms as 
much as possible, so that the number of terms to be 
carried is small. This procedure leads to a set of mz & n 
linear algebraic equations in the mm X a» unknowns 
A which generally has the additional property (i 
the trial solution is not too poor an approximation 
that the terms in the main diagonal (4,,,,) are pre 
dominant. Thus the procedure given by Timoshenko’ 
can be used to compute the 4, This procedure will 


be illustrated below with m l,z~=3. Let the system 


pe 


ayAy T (pA T 3A 1, = bi | 
(oA 1) T (9A 12 T (193A 1, = bo 


anAi T ApA tT a Ady »,\ 


where the a,; are known constants found by the integra 
tion. Now under the assumption of the importance 


of the main diagonal terms let 


A," = b) ay (S) 
and 
Aj." => (bd, (1,A 3)") Qo» (9) 
A,3° = (bs (131A y)' 132A 32") /A33 (10) 
and by resubstitution into Eq. (7) a second approxima 
tion 1S 
A)! = (d, QypA 1" 143A 3) ay, | i 
A})! = (be ~ AA}! - dle A 3° (do § 


and so forth. The rate of convergence is dependent 
upon the relative unimportance of the coupling terms 
and in the cases tried, the third iteration was adequate. 


EXAMPLES 


Calculation for Particular Static Loadings 


It is desirable that the method outlined above be 
checked by experimental data, and such data had been 
obtained by the M.I.T. Aeroelastic and Structures 
Laboratory.” Fig. 2 shows the experimental and cal 
culated influence surfaces for two load positions on a 
flat plate, and Fig. 3 shows the same results for a double 
wedge cross section. The approximate solution shows 
the correct curvature variation, but overestimates the 
centerline deformation in the flat plate case and under 
estimates the centerline deformation of double wedge 
by about 8 per cent. It will be noted in Fig. 3 that 
the agreement is poor at the edges. Some qualitative 
estimates indicate that the convergence is slower at 
the edges, due to the taper. It is expected that in- 
creasing the number of terms would improve the agree- 
ment in the calculation of the influence surfaces. Fig. 
4 shows a comparison of the deflection surface of a 
uniformly loaded double wedge and as is expected, the 
integration improves the overall accuracy. 
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On the basis of these results, it appears that the 


method is reasonable. 


Calculation of Lift Due to Elastic Deformation* 


Now the method will be applied to the problem of 
computing the aerodynamic induction effects upon the 
loading, deformation, and the aerodynamic coefficients 
themselves. 

Symbolically, the process may be written as follows 
Let the local angle of attack distribution be a,, where 


Ow Ov 13 


and assume that the load at each point is proportional 


to the local incidence 
Ap q pix, v) [a + (Ow Oy) ] l4 


It will be observed that the derivative of w appears in 
Eq. (!4), rather than w. Since the basic solution, Eq. 
(2), can be written in terms of (Ow) (Oy), as well as 
w, and since it represents a simplification to write the 


solution in this way, then w(x, v, vo, Vo) will be replaced 


by (Ow) /(Oy) (x, V, Xo, Yo), So that Eq. (2) becomes 
Ou ” OW} 
a P(xo, ¥ dxod yo + 
Oy JJS Oy 
gs OW; Ow zs 
pi Xo, Vi dx. dy » 
JJS OV OVo 


or, letting 


and 


K(x, ¥, Xo, ¥ P(Xo, Vo) (Ow, Ov 


Ow O%> wi Ow . 
a + I K(x, x0, ¥ dxody 16) 
Oy Ov JJs "OM : 


and upon applying the resubstitution technique? 


* After the work described below was completed, in English 
Electric Company, Naval Projects Division Report L.A.t.038, 
entitled ‘‘The Theory of the Aeroelastic Load Distribution an 
Deformation of a Lifting Wing at Supersonic Speeds’ of Oct. 6, 
1952, by B. W. Bolton-Shaw came to the attention of the author 


Dr. Bolton-Shaw’s method of approach is essentially the same as 


that given below; however his solution for the particular problen 
of a supersonic wing is finally written in the form of a difference 
equation involving an experimentally determined influence sur 


face 


tT If the cross sections come to a sharp edge and the loading Is 
finite at the edge, the slope of the deflection surface is infinite and 
the resubstitution will not converge. In fact such a curling has 
been observed at the NSL during wind-tunnel tests on an adia 
batic flat plate with a carefully honed leading edge. Generally, 
this effect will not be encountered, and there are several obvious 
methods of avoiding this singularity—i.e., compute the kernel 
for a finite thickness at the edges 
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AN APPLICATION OF 


, Ow» al ; Ow» ae 
Ort = a a } K dx dy i p B I x 
OV ¥ Jd S OV t=2e S 


K \, vy, Xo, Vo } | K ADs ay » X2, Ve Pees x 

-_ Ow . , = 

BASe-y, Vi Ke. 9 dx dy; ...dxody 14 
. ov 


so that after substituting in Eq. (14), the increment 
in the pressure coefficient due to the flexibility is 


A Ap Ow OW» 
P\%X, ¥ ap(x, Vv) T 
qd i Oy a Ov 
"Cc .. Cs 
I K dxodyo + .. (1S 
Ja s OVo 


[he change in the aerodynamic characteristics can 
easily be found by integrating the local loading incre 
ment. Inspection of Eq. (1S) indicates that the effects 
of the static elastic deformation on the aerodynamic 
characteristics is linear in the angle of attack and dy- 
namic pressure and involves a function of 1 D that 
here appears as a power series in (1 PD)", since A in 
volves 1D. 

In linearized supersonic flow, p(x, y) may be ap 


proximated by 4/V J/* — 1 in the two-dimensional 
region, and 4(6 — x) (c — y) in the tip region. The 
ratio Cr, C, was computed for an aspect ratio 


a “rigid 


unity flat wing at J = y/2 to the order of (1 D)* to be 


C k 
e 14 q 


: where k = 
ULaris (t/c)?10" 


0.319 for steel wings 19 
1.055 for aluminum wings 

rhe largest source of error in the above calculation of 
C;, is due to the fact that the linearized pressure dis 
tribution does not place the center of pressure at the 
proper point. The accuracy of the calculation could 
undoubtedly have been improved by using Busemann’s 
second order term 

Ap . = C\ a+ (Ow Oy r C2 a+ (Ow- Ov 

ind retaining only terms of order (Ow) (Oy), so that 


Ap q Cra + Cra? + (C, + 2Cra) (Ow Ov 


If this were done then the structure of Eq. (18) would 


be 
Ap Ap ‘ OWls 
= a(C; + * a T 
J flex Qrigid Oy 
. j , ; Ps SE 
C, + 2C2a) (CG, + Cra K dS + 
S Ov 


which is nonlinear in a. 


CONCLUDING REMARKS 


rhe examples have all dealt with rectangular plates 
that are either flat or have a double wedge cross section. 
The generalization to arbitrary plan forms can be 
effected, usually, by the introduction of a coordinate 
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system which is geometrically similar to the wing plan 
form, and by writing the energy equation in this co- 
ordinate system. It is not difficult to make the variable 
change ; however, the form of the equation is fairly com- 
plicated in the general case. The details for specific 
plan forms are given by Fung? for a sweptback wing 
and Stein ef a/.® for delta wings. 


APPENDIX 


The complete set of functions that were used to make 
the calculations shown in the figures are listed below 
for the applicable cases.!. It will be noted that the 
semispan of the surfaces is 6 and the semichord is unity. 
The terms w,,, in Eq. (6) are written w,,,(%, Vv) = 


fin(X)2,(¥)! ? as shown below 


Case 1. Uniform Thickness 


] > me 


rm I 
(x) = — cosh —— cosh 
} 
. , Pm . Pm* 
o» \ sini sin 
h f) 


D ) ) 
where 
cosh P, tT COS Pp, 


sin ~p,, + sinh p, 
and the p,, are the roots of 
cosh p cos p + 1 0 


fm(x) represents the clamped-free vibrating-beam mode 


shape. 
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cosh gy, COS gny + COs g, cosh g,y 


ny) 
V cosh? g,, 


+ cos” Qn 


if m is even, and the g, are the roots of 


tan g + tan hg 
qo 

When x is odd 
sinh gy, sin gyy + 


£rAAY) = 
V sinh? q,, 


and gq, is defined as before. Note 
L(y) = 1's 


) 


gly) = (3/~ 


0 
= 0 


sin g, sinh g,Vv 
— sin’ g, 


that 


~ 


2)y 


Case 2. Rectangular Double Wedge Surface 


fm(x) are the same as in Case 
results! 
G,Ji[2¢n \ ( l 
V(1l- 


£n(¥ 


1. Using Nicholson's 


if 7 is even and G, is a normalizing factor where 2g, are 


the roots of /2(2q,) = Oand 
V 2-29, l 
G,, = q ° 
V (2qn)? — 1 Si(24n) 
For the odd 
— SIiP2anrV 1 -— |i|)] 
Gn\V) = G, 
Vi(1l- v1 
J\(2qn) Mi[2¢,V (1 — y)] 
11(2qn) Vv (1 Vi 
and 
2n(—y) = —Bnly) 
The g, are the roots of J2(2q,)J2’(2q,) — Je'(2qn)12(2qn) 
= ()and 
‘ IoleJ;? ae JoJol;? z JI = J ToS I 
G, = + 


9 


I, 64 2¢n 


and the argument of the Bessel fun 


rf 


) 


(4 1)/2 


ction is 2g,, J, is the 
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Bessel function of the first kind, mth order with rea] 
argument, /, is the Bessel function of the first king 


nth order with imaginary argument. 


Note that here g(y) = 1 and gi(y) = 3y. The 


fm(x) and g,(y) have been made normal and orth 
onal—.e., 


/ 
| Fim(x)fn(x)dx = 6,” 
0 


od | 
~-. SL mm n 
R(v) omn(Wgn(widy = 4, 
- : a (O mtn 


where &(y) is the appropriate weighing function. For 
Case 1, k(y) = 1. For Case 2, k(y) = (1 y 
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Analysis of Peak-Holding Optimalizing 
Control 
H. S. TSIEN* ano S. SERDENGECTI?* 
California Institute of Technology 


SUMMARY 


[he peak-holding optimalizing control is analyzed under the 
yssumption of first-order input linear group and output linear 
group. Design charts are constructed for determining the re 
quired input drive speed and the consequent hunting loss with 
specified time constants of the input and output linear groups, 
the hunting period, and the critical indicated difference for input 


drive reversal 
INTRODUCTION 
INVENTED BY C. S. 
Their 
In almost 


nae CONTROL WAS 
Draper, Y. T. Li, and H. Laning, Jr.! 
basic idea can be summarized as follows: 
all engineering systems, within the restrictions of oper 
ation, there is an optimum state of the system for per 
formance. For instance, in an internal combustion 
engine, within the restriction of producing the load 
torque at the specified speed, there are optimum set 
tings for the manifold pressure and the ignition timing 
for minimum fuel consumption. Another example is 
an airplane under cruising condition; then under the 
restriction of engine cruising r.p.m. and assigned alti- 
tude, there is an optimum combination of trim setting 
and engine throttle for maximum fuel economy or maxi- 
mum miles per gallon of fuel. But more important 
than the existence of an optimum operating state is the 
fact that 
exactly predicted in advance because of the natural 


the optimum operating state cannot be 


changes in the environment of the engineering system: 
In the case of the internal combustion engine, it is the 
changes in the temperature and the humidity of the 
air; in the case of the airplane, it is unavoidable changes 
in the aerodynamic properties of the airplane and the 
engine performance with age. Therefore if the purpose 
is to operate always near the optimum state in spite 
of the “drift’’ of the system, then the control device for 
the engineering system must be so designed as to search 
out automatically the optimum state of operation and 
to confine the operation close to this state. This is the 
basic idea of optimalizing control. 

The application of Draper's optimalizing control 
to the general cruise control of airplanes was discussed 
by Shull. Shull emphasized the possible elimination 
of extensive flight testing of new airplanes for perform- 
ance determination, because the optimalizing control 
will automatically measure the performance whenever 
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the airplane is flown. This in itself would constitute 


a great saving. But moreover, in critical circumstances 
such as flight through icing atmosphere, the ability of 
the optimalizing control to extract the best perform 
ance of a radically changed system (through ice deposi 
tion on the airplane) could be of utmost importance. 
There are two fundamental problems in the theory of 
optimalizing control. the the 
dynamic effects of the controlled system on the per 
The other problem is the 
rhe two prob 


One of problems is 


formance of the control. 
elimination of the noise interference. 
lems are somewhat interrelated, because if large devi 
ations from the optimum state or the optimum oper 
ating point and hence large loss can be tolerated, then 
the noise interference will not be critical. The basic 
design aim of optimalizing control is to have the smallest 
loss or to operate as close to the optimum state as 
possible without the danger of having the control 
misled by the noise interference. Both of these prob 
lems were considered by the original inventors of opti 
malizing control. The noise problem is essentially the 
problem of detection of a sinusoidal variation under 
heavy random interference, a subject of much current 
research. The purpose of the present paper is to solve 
completely the first problem of dynamic effects under 
the assumption that the dynamic properties of the con 
trolled system can be approximated by a first-order 
linear system. We shall begin with the brief review 
of the operating principles of an optimalizing control 
of the peak-holding type a type least affected by the 


noise interference." ° 


PRINCIPLE OF OPERATION 


The heart of an optimalizing control system is the 
nonlinear component that characterizes the optimum 
operating condition of the controlled system. For sim- 
plicity of discussicn, it is assumed that this basic com 
ponent has a single input and a single output. For 
the time being the dvnamic effects will be neglected 
and the output is assumed to be determined by the 
instantaneous value of the input. Since there is an 
optimum point, output as a function of input has a 
maximum at the output \y at the input x), as shown in 
Fig. |. It is convenient to refer the output and the 
input to the optimum point and put the physical input 
as x + x and the physical output as y* + yw. The 
optimum point is then the point x = y* = 0. The 
purpose of an optimalizing control is then to search out 
this optimum point and to keep the system in the 1m 


mediate neighborhood of this point. In this neighbor 
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hood, the relation between x and y* can be represented 


as 
y* = —kx? (1) 


where & is a characteristic constant of the controlled 
system. 

The operation of a peak-holding optimalizing control, 
neglecting the dynamic effects, then would be as 
follows: Say the input x is below the optimum value 
and is thus negative. The input drive is then set to 
increase the input at a constant rate. At the time 
instant | (Fig. 2) the input changes from negative to 
positive and passes through the optimum point. The 
output y* is thus maximum at the time instant | and is 
decreasing after the instant |. Now if an output sens- 
ing instrument is so designed as to follow the output 
exactly when the output is increasing, but hold to the 
maximum value after the maximum is passed and the 
output starts to decrease; then there will be a dif- 
ference between the reading of this output sensing 
instrument and the output itself after the time instant 
1. This difference is shown in the lower graph of Fig. 
2. When this difference is built up to a critical value 
c at the time instant 2, the input drive is tripped and the 
direction of the input drive is reversed, but still at the 
same constant rate as before. After the instant 2 
then, the input decreases and the output increases till 
a maximum in output is again reached at the time 
instant 3. At time instant 3, the input, of course, 
again passes from positive to negative, and the indi- 
cated difference between the output sensing instrument 
and the output itself again builds up. At the time 
instant 4, the difference reaches the critical value c 
again, and the input drive direction is again reversed. 
At the time instant 5, the input + becomes zero again 
and another maximum of the output is reached. The 
period of input variation is thus the time interval from 
the instant | to the instant 5, and the input, when 
plotted as a function of time, consists of a series of 
straight line segments forming a saw-tooth variation. 
The period of output variation is the time interval from 
the instant | to the instant 3, and the output, when 
plotted as a function of time, consists of a series of 
parabolic ares. The periodic variations of input and 
output are called the hunting of the system, and the 
period of output variation is called the hunting period 
7. The period of input variation is thus 27°. 

The extreme variation of output A (Fig. 2) is called 
the hunting zone. If a is the amplitude of the saw- 
tooth variation of the input (Fig. 2), then due to Eq. 
CT), 


& = ka (2) 


The difference between the maximum output and the 
average output of the hunting system is called the 
hunting loss D (Fig. 2). Because of the fact that the 
output is a series of parabolic arcs, 


DD = (1/S)4.= (1L/3)ka* (3) 


For this idealized case, the critical indicated difference 


c between the output sensing instrument and the out 
put itself is equal to A, the hunting zone. It is they 
clear from this discussion that in order to reduce the 
hunting loss for better efficiency of the system, ong 
must try to reduce the hunting zone or the amplituc 
of input variation., Unfortunately the critical ind; 
cated difference is also reduced by such modification 
and a limit is set by the noise interference on the proper 
tripping operation of the input drive. 

The dynamic effects are so far neglected. But in am 
physical system, this is not possible because of the ever 
present inertial and damping forces. The output y* 
given by Ey. (1) has to be considered then as the fic 
titious “potential output’ but not the actual output 
y measured by the output indicating and_ sensing 
instrument. y* is equal to y only when the period 7 
of hunting becomes extremely long. The relation 
between y* and y is determined by the dynamical 
effects. For the conventional engineering systems 
these dynamical effects are determined by a linear 
relation. For instance, in the case of an internal 
combustion engine, the potential output is essentially 
the corrected effective pressure generated in the engine 
cylinders, while the actual output is the brake mean 
effective pressure of the engine. The dynamical effects 
are here mainly due to the inertia of the piston, th 
crankshaft, and other moving parts of the engine. 
For small changes in the operating conditions of the 
engine, such dynamical effects can be represented as a 
linear differential equation with constant coefficients 

Since the reference level of input and output is taken 
to be the optimum input ., and the optimum output 
yw, the physical potential output is y* + yo and the 
physical actual output is y + y. Thus the relation 
between the physical potential output and the physical 
actual output can be written as an operator equation 


y+ vw = Fi(d/dt) (y* + yo) } 


where /) is generally the quotient of two polynomials 
in the time differential operator d/dt. In the language 
of the Laplace transform then F,(s) is the transfer 
function. Let the linear system which transforms the 
potential output to actual output be called the output 
linear group. Then F,(s) is, specifically, the transfer 
function of the output linear group. Ry implication 
however, when the dynamical effects are negligible or 
when s = 0, the potential output is equal to the actual 
output. Therefore 


F,(0) = 1 5) 


Since the optimum output y) only varies extremely 

slowly by the drift of the controlled system, during a 

time interval of many hunting periods yp) can be taken 

asaconstant. Then the condition of Eq. (5) simplifies 
Eq. (4) to 

y = Fi(d/dt)y* (6 

In a similar manner, let «* be the ‘‘potential input 


that is actually the forcing function generated by the 
optimalizing control system but not the actual mput 
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x. Itisx* that has the saw-tooth form shown in Fig. 2, 
but not x. The relation between x* and x is deter- 
mined by the inertial and dynamical effects of the in- 
put drive system. This input drive system can be 
called the “input linear group” of the optimalizing 
control. The operator equation between the potential 
input «* and the actual input x is 


x = F(d/dt)x* (7) 


F(s) is thus the transfer function of the input linear 
group. Similar to Eq. (5), the meaning of potential and 
actual inputs implies 


FO) = 1 (8) 


Thus a simple representative block diagram of the 
complete optimalizing control system can be drawn as 
shown in Fig. 3. The nonlinear components of the 
system are thus the optimalizing input drive and the 


controlled system itself. 


FORMULATION OF THE MATHEMATICAL PROBLEM 


The general relation between the input x and the 
output y is determined by the system of Eqs. (1), (6), 
and (7), with the potential input «* specified as a saw- 
tooth curve with period 27 and amplitude a. Let 
w, be the hunting frequency defined by 


w = 2r/7 (9) 


then «* can be expanded into a Fourier series, 


a Sa ox (—1)” : wot 
= = Sm (2H = 1) 
7” n=0 (2n si is 2 
. n 
Sa = Dea! [(2 + 1)/2 |iwot 
= (e = 


T° n=0 (2n + Dis 21 


{[(2n+1)/2 |iwot 


) (10) 


Therefore by using Eq. (7), the actual input x is given 


by 


Sa (—1)’ 


By using Eqs. (11) and (16), the actual output y js 

given by 

LGa*k ~~ - (—1)° 
a _X 


rT! 1-0 mao (2a + 1)? (2m + 1 
a 4 
_ .  fzararil. 
) Pol(n + m + 1)tax|F lwo} X 


Fy{(n — m)ia | X 


Bi - 


By comparing Eqs. (11) and (12), it is seen that the in 
put has half the frequency of the output. This is, of 
course, to be expected from the basic parabolic relation 
of input and output as specified by Eq. (1). 

The average of the actual output y with respect to 
time ¢, being here referred to the optimum output 5, 
gives directly the hunting loss D. Equation (12) shows 
that this average value is the sum of terms with ” = m 
from the second and the third terms of that equation 


Therefore, using Eq. (5), 


32a7k | 2n + I 
D= 2. F tun ) X 
ri =) (2n + 1) ( 2 ) 


This equation can be easily checked by observing that 
when the dynamic effects are absent, F; = 1, then the 
series can be easily summed and D = (1 3)a°*k as te 
quired by Eq. (3). Equation (13) also shows that the 
average output and hence the hunting loss are inde 
pendent of the output linear group. This agrees with 
the one's physical understanding: Only detailed time 
variation of the output is modified by the dynamics of 
the output linear group. In the case of an internal 
combustion engine, the average output specifies the 
power of the engine. The dynamics of the output 
linear group is determined by the inertia of the moving 
parts. The power of the engine is certainly independent 
of the inertia of the moving parts. 

Equations (11) to (13) fully determine the perform- 
ance of the optimalizing control system once the 
values of a, k, and w) are specified and the transfer 
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functions / and F,(s) of the input linear group and 
the output linear group are given. The following sec- 
tions give the detailed calculations and results for the 


case of first-order input and output groups 


FIRST-ORDER INPUT AND OUTPUT GROUPS 


[he frequency w) of the optimalizing control is 
usually low, and the important dynamic effects come 
from the inertia in the input and the output linear 
sroups. Then these linear groups can be closely ap 
proximated by first-order systems. In other words, 


their transfer functions are 


Fi(ww) = 1 (1 + twr (14 
F lw l \ l i lw T { ] 5 } 
where 7; and 7 are the characteristic time constants 


of the input linear group and the output linear group, 
respectively. It1is evident that these transfer functions 
satisfy the conditions of Eqs. (5) and (8). 

By substituting Eq. (14) into Eq. (11), the actual out 
put 1 is given by 
Sa 1)? f eli2nt 2 Jia 
wr? 9 2t(2n + 1 + (2n + 1)t(ay7;/2 


(16 


When the summation is carried out, Eq. (16) yields the 


following equations for the input xv: 


r t/T)/(7i/7 
= [ T Ti € 
7 / ZT cosh (7/27; 


for — : (17a 
») 


To calculate the actual output y, both Eqs. (14) and (15) have to be substituted into Eq. (12) Le. 
ET IAR. eas. ems - cs 
S S x 
us i mga (2n + 1)? (2m + 1 
t ? l)w 
j e 
{1 + (a + m + 1)tworo] [1 + (2 + 1)t(wor;, 2)] [1 ” l)t(ayr; 2 
piinan 
1+ (n m)tooto) [1 + (2m + 1)i(wor; 2)] [1 2m l)t(a@y7; 2 
- " 
(n m)twyto| [1 (2n + 1)t(wor;, 2)] 1 2n l)a(ayr; 2 
, , 1) wol 
é { 2] 
ll — (mn + m + 1)taor (2n + 1)4(wor; 2)] [1 ” l)i(wpr, 2)]9 
By changing the summation indices, Eq. (21) can also be written as 
17°2.\ 2p aan = eo 
a ‘ oe 
7: gs (1 + 1Sw T\ 
s l 
5 \ 
nog (2n + 1)2[(2m + 1) — 2s]? [1 + (2n + 1)t(wor; 2)] [1 — [(2n + 1 2s ]i(wyr; 2); 
1 i 
> | é l ) 
(1 isapto) 29 (2m + 1)2[(2m + 1) — 2s]? [1 — (20 + 1)t(wor;/2)] }1 + [(22 + 1) — 2s]t(wors/2)} 
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and 


i ® ( 4 T e! ' ] 
x= —-NT|—-{[{14 + 
1 1 T cosh (7° 27;) J 


for l7b 
») /7 ») 


where .V is the constant input drive speed —1.e., 
A 2a 1 1S 


By using these equations, the variation of actual input 
x with respect to time can be calculated for any speci 
fied data. Examples of such calculations are shown 
in Figs. 4 and 5 for 7; 7 0.1 and 7; 7 0.4, respec 
tively. Both show the expected effect of rounding 
off of the sharp corners of the saw-tooth curve and a 
time delay. It is of interest to note that while the 
delay is almost equal to 7, itself for small 7, 7, the de 
lay is less than 7; for larger 7; 7 
With the first-order transfer function of Eq. (14 


the hunting loss given by Eq. (13) becomes 


29725 
D = 32a*k a | 


hast << p ») ool 
rT? 4=09 (2" + ) + 2n + 2 |7w*T it 
(19 


By carrying out the summation, Eq. (19) gives the 


hunting loss as 


D = (N°T°R 12) [1 27,/7 
24(7,; 7)? tanh (7° 2r 20 


Fig. 6 shows a dimensionless plot of this equation. 
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SNe) OS. (—1)*~ "eo" 
74 n=o (2% + 1)4{1 + (20 + 1)? (wor: /2)?] s=1 (1 + tsar) 
o [(22 + 1)? + 457] [(1 + tenzis) + (wor; /2)? (20 + 1)2] + 8(wo7i/2)2s2(2n + 1 
=o (2a + 1)*[(20 + 1)? — 487]? [1 4+ (wor; 2)? (2n + 12) [CL 4+ tworis)? + (wor;/2)? (20 + 1)?] 
l tswol 
= (=o ; = 
> zi, x 
=< (1 — 1t5w79) mpers 
((27 + 1)? + 4s?] [C1 —- tanris) + (wor, / 2)? (2n + 1)?] + S(wo7; 2)2s2(22 + 1 | 
ae 
(2n + 1)2[(2n + 1)? — 452]? [1 + (cpr, /2)2 (2m + 1)2] [(1 — tworis)? + (wor; 2)? (Qn + 1)2]f O°? 
The last two summations in Eq. (22) are complex conjugate of each other, thus 
ST?N°k us aay ee 
y= e. % + 2R1 > > & 
Z { ) a st ‘ { 9 ° )\2 - Ea <2 \ fw 
T nao (2m + 1)4 [1 + (2n + 1)? (wor al s=1 (1 + t5a070) n=0 
[((2m + 1)? + 45?] [C1 + tworis) + (ao7;/2)? (22 + 1)?] 4 S(wor;/2)2s2(2n + 1 } 
99 
» » ae it > > . , 9 } -0 
(2n + 1)?[(2m + 1)? — 4s?]?[1 + (aor; 2)? (2n + 1)2] [C1 4+ taoris)? + (wor; 2)? (2n + 1)2)8 


where RI means the real part of the expression following it. In order to carry out the summation with respect to 


the index 7, Eq. (23) is resolved into the following partial fraction form: 


ST?N7k l .~ (-—1) aaa 
y= : — : > + 2Ri > x 
i 0 (2n + 1) [1 + (27 + 1)° (Wy T 2)? } 1 | 1 + tsaT 
l > l . (wy7;/ 2 \- | 
(452(1 + iswor “29 (2a + 1)? 2(1 + tsuo7;)? [1 + ts(a7;/2)] »xo [1 + (wor; 2)7(2n 4 
(QT 2 —~ 
S 
2(1 + tsey7;) [1 + ts(a@o7; 2)] ,=0 [A + tworis + (wr; 2 n+ | 
+ t8a97; + 4(wor;/2)?2s?7] A 2 "a ‘ 
[1 t tSw97; + 4(wo7;/2)"s s l , 2 - | 94 
bs?(1 + tsay7;)? yao (22 + 1)? + (22s)? 1 + tseor;) soo [((22 + 1)? + (42s)2]79 


By using the summation formulas given in the Appendix, the sums with respect to 7 can be evaluated and the result 


is, noting that tan rs = 0 for integer values of s, 


ST-°N°k 7 T- fouti\* wT foot: \° T 
y= — — ( + tanh + 
° 14 96 S 2 | 2 WoT 


(—1)*~ "eS (x? l T (wor ;/2)* tanh [7 (wo7;) | | . 
ORI ; . " +h Bh : 95 
1 (1 + aisworo) 14 (482) (1 + tswor $ (1 + iswori)? [1 + ts (wor) /2)]9, 


25) is again resolved into partial fractions in order to carry out the summation with respect to s, viz., 


ST?N°k [ = i (2) T (* ) r | 
y = , = + tanh + 
= 1 OG \ 2 | 2 WT | 


| (wo7o/2)? j 2(wo7;/2)*? tanh [2 / (wo7,) | | b> —] 5 ia 
: Te IX : T 
2 \[(wo7e 2) — (wor /2)] U[ (wore 2) — (wor; 2) [wore — (wot; /2) ] 1 (1 + tSw7 
2 lise 
(a7 ;/2)* 2(wo7;/ 2)? tanh [7 /(wor7;) ] ) (—l)" ‘e T [wor WT 
a RI >> ; ; x 
[(anro 2 — (ay7;/2) | [ (wor 2) — (WoT 2) } J 1 | 1 + 1say,7 2 z = 
lL» iseoot lise 1 iswot 
(—1) te 7 — (=<) € (we7;/2)* tanh [2 /(wo7;) | a -] e 
RI >» + — RID : - MZ 
# | , § [2(wot) 2) — (wor; /2)] 1 (1 + ts(wo7i/2) | 


I 


(wore 2) — (an7;/2)] 1 (1 + tswor 


2(wor; 2)1 tanh [7 /(wo7;) | (—] e | ; 
RI >> | 265 
The result of carrying out the summations in Eq. (26) and simplifving the expressions is, 
T2N2 [ jl (=) (; 1 >) ri | l (5) | TiTo ay] 
, = 20 ae a i pees ne nay ie aan oe a 5 
: | \2\7 7 T/\7 2\7 T° € ah 
UT (rT 


l (7o/T)? j 2(7,/T)* tanh (7°/27;) ) Le 
(To 


2 T— 7:17) U[(r0/T) — (7:/T)] [2(00/T) — (7:/T)] sinh (7° 27 
jt ‘ (7,/T)? \ w/e (7,/T)' e7 (2 T)Mai/1 )] has 
T)] cosh (7 /27;) [2(r70/T) — (7;/T)] cosh*(T/27 i 


(7 ((mo/T) — (4:/T))9 (re/T) — (4; 
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ile U 
or ve | 3G) - (: 4. ta (5) +15 (5) + iT + (%) + To | Te a “ 
y L tar rr’ Ff T 2 r) T? T r °F aF’ 
| ( (7o/T)? j 2(7,/T)* tanh (7°/27;) a Pr tiabledieche 
2 [((m/T) — (7:/T)] UlCre/T) — (7i/T)] [2(90/T) — (7:/T)] 1 sinh (7/270 
lf (7,/T)? Si (r,/T)%e8 VT em (7,/T)* , ' )| (7m 
7 trol) — (7:/T J ((7o/T) — (7:/T)] cosh (T/2r (2(70/T) — (7;/T)] cosh? (T/27;)/ J ™ 
for 1/2 S$ t/T & 3/2. 
In Eqs. (27a) and (27b), there are apparent singularities whenever 7 7) = 7, 7 and 27 / 7, 71; that is, the 


value of output y seemingly cannot be determined for these values of time constants. However this is deceptive 


By using a simple limit procedure or by direct evaluation of Eq. (25) for these two cases, it can be shown that this 


isnot the case. For example, for 7 l= %/1T 
rewre | — 3 (5) +2(3)G) - 33) +1L- G) - 3) 
2T?N7k < — 2 ee 2 Et eae =) + 
\ ee 7 7 24h2 7 Se 
T t | ; i 3 i l ‘— . T 
( \( ) a |e UDMA TY sech + (“*)e UDMTT) osch + (=) ” headed ” sech? (28 
T/\1 8 2r, 2\7 2r, 2\T/ 27,4 
for —(1/3) & 8/T S 1/2, and, for 2%,/7 = -,/f, 
rvs (0) —2()() 205) = LGN0) #209) #1) oa? 
y= —2Z7°NR - + : i 2 ‘tnh +- 
lo\r 2\7)\r] 7 4\7r T/\T r nw 
rs [(2t/T)/(ri/T) ] (=) ( t =e ((¢/7 talialded 
ctnh ——- i — 2s = £— = (29 
Ir cosh? 7/27; 7 7 1 cosh 7/27, | 


for —(1/2 S £/7 3 1/2. 


An analysis for the continuity of Eqs. (27a) and (27b) at t/ 7’ = 1/2 shows that the values of y and its derivative 
with respect to ¢ are the same at ¢/ 7 = 1/2 whether they are computed from Eq. (27a) or from Eq. (27b). 
Now the computation of the potential output y* can be accomplished by letting 7/7 = 0 in Eqs. (27a) and (27b). 
Thus, 
Dene Let T t l/r t r oan Fe b falta” rows 
y* = 21°N*k | — i eet FF Be 4 Bae tees |g = _ & a (30a) 
2a 1 1 i 7 7 T/ cosh(7'/27;) ZX cosh?(7°/27;) 
for — 2) S$ #/T S 1/2 and 


voarwed (+ (NC) 205) -G)LE- 


; _ 
(" ie ri/1 (2) t/1 wey) 
ee I i _ ee — (30b) 
7 cosh (7/27;) 2\T/ cosh*(T/2r,) § 
for 1 2 t/z = 3/3. 


These expressions check with the result of direct calculation of y* by Eqs. (1) and (17). 

Figures 7 and 8 show the dimensionless plots of actual output y and potential output y* for the particular values of 
™ Tand7, 7. In these figures it is clearly seen that the dynamic effects not only decrease the output of the system 
but also introduce a time lag and lower the maximum output of the system. Figure S with 7, 7 = 0.4, 7% 7 = 0.6, 
has the maximum value of y almost at the very instants of input drive reversal points, ¢ 7 n+ (1/2 This is 


indeed an extreme case 


DESIGN CHARTS 


From the principle of operation of the peak-holding optimalizing control, it is seen that the most important quantity 
to be specified for its design is the critical indicated difference c between the reading of the special output sensing in- 
strument and the output itself. By definition, c is the difference of the maximum of the actual output y and the value 
ot y at the tripping instant of the input drive. The instant of reversing the input drive is typified by ¢ 7 = 1 2. 
If the corresponding instant of maximum y is /*, then the critical indicated difference c is calculated as 

c = y(t*/T) — y(1/2 31 
by using any one of Eqs. (27), (28), or (29). Since the instant of input drive reversal must come after the instant of 


maximum output, f* 7 < 1/2. 
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To determine /*, one may use the condition of zero slope—1.e., dy dt = 0. Then Eq. (27a) gives 


f T ry] to/1 j 2(7:/T )* tanh (7/27 at 
& / Til Ms./F - (7,/T)] Ul[(ro/T — (7;/T)) [2(r0/T (7, 1 A sinh (7/27 
" = (CC \ r,/1 i r:/1)? ¢ 
{ eo 4 to/1 r, T)\9 T/1 r,/I°)| cosh (7'/27 
r;/1 é ai 
0 (32 
2(70/T r,/I')| cosh*(T'/2r 
[his transcendental equation for ¢* 7 may be solved by iteration. For instance, for small 7 7 and 7, 7, only 
terms within the first brackets are of importance, then ¢*/7°~ (7) + 7;) 7. This is already recognized by Draper 
ind co-workers. The complete results of calculation are shown in Fig. 9, which shows that ¢* 7 is almost only a 


function of (7) + 7;) 7 with minor modifications from the parameter 7) 7;, the ratio of characteristic times of the 
output linear group and the input linear group. Values of t*, 7 beyond | 2 are not shown, as clearly then the maxima 
{the output will occur after the corresponding input drive reversal points and proper operation of the control will be 
lifficult if not imps ssible. 

With /* 7 determined, Eq. (31) gives c by substituting Eq. (27a). However the specified quantities of an opti 
malizing control are k, the characteristics of the controlled system, and 7;, 7, the characteristics of the linear group 
From considerations on the noise interference, the designer can make an appropriate choice of the period 7’ and the 
critical indicated difference c for input drive reversal. Therefore the quantities that the designer wishes to know, 
iter he has the values of k, 7;, 7», 7, and c, are .V, the input drive speed, and D, the hunting loss. Thus the result 


of calculation with Eq. (31) should be written as follows 


TN ({ ("| (“)(" VG ) (ro/7;)? (7;3/T 
ea ens eee a 2 be + | li “ x 
/c/R | 4+ \7 | Tj] \7r (4 5 [(7o/7;) — 1] sinh (7;/70)(7 / 27 


2(7;/T) tanh (7'/2r | s/T)/lae/es) (xs/7 . 
= 7 c (e e 4 
Ul (rm 7) — 1) [2(70/ 7) -— 1] 
2(7,/T j| ¢* (e,/7 IT) a /T l (7,/1 ] T/22,) 
aie ) cae a e ape’ T € ( 
[(7» 7i) — 1] cosh (7'/27,;) (7 (to/7i) — 1] 2 (t%/T | 
, Qe/T)/(ri/T Pi 1/2 
(7:/T)*(e é 
, 33 
[2(70/7i1) — 1] cosh? (7/27 
When .V is determined, Eq. (20) then gives the hunting loss D. 

Figures 10 and 11 are the design charts for peak- signer should improve his input drive system so as to 
holding optimalizing control computed from the equa- reduce the constant 7; 7) is, however, a constant 
tions of the preceding analysis. Figure 10 gives of the intrinsic characteristic of the controlled system, 
TN Vc k as a function of 79/7; with (rm + 7,)/T as due to, say, the inertia of the moving parts of the sys 
parameter. Figure 11 gives relative hunting loss D « tem. 7») is thus not at the disposal of the designer of 
again as a function of 7) 7; with (7) + 7,) 7 as param- the control system. However, suppose there is a com 
eter. The peaks of curves near 79/7; = 1 indicate a pensating circuit between the output y and optimalizing 
sort of resonant effect between the input linear group input drive unit (Fig. 3), such that the effects of the out 
and output linear group. The hunting loss for fixed put linear group is completely compensated. Then 
tT; + 7) Zand c is smaller for 7) 7; away from unity. the effective signal for input drive reversal is not the 
For fixed 7;, 7», and c, clearly the way to reduce the actual output y, but the potential output y*. In other 
hunting loss is by increasing the period 7. words, the value of 7) is made to be effectively zero. 


Even if complete compensation is not achieved, the 

CONCLUDING REMARKS effective value of 7) can still be greatly reduced. For 

1 ie ; difficult cases then, such a compensating unit should 
rhe present analysis gives the necessary input drive : ass 
; ; certainly be added to reduce the hunting loss. This 


speed .V and the hunting loss D for any specified hunt a fi 
: will be just a minor complication when compared with 


ing period 7, time constants 7; and 7 for the input oe ; 2 a 
the additional equipment required for satisfactory 


linear group and the output linear group, and the ey aie 
| ae rs pee = noise filtering. 
chosen critical indicated difference c. 7 and ¢ are 

hxed by considerations on the noise interference. The . 
na ' PG Appendix 
analysis shows that whenever the hunting period is rela 
tively short with respect to the time constants 7;, 7», és ; : 
} Bee : : ['YPICAL SUMMATION FORMULAS 

or whenever (7; + 7 )/7 is relatively large, the hunting 

loss will be large, especially when 7; and 7» are nearly Rl and Im mean, respectively, the “real part of” 


equal. To avoid such unfavorable condition, the de- and the “‘imaginary part of”’ the expression following it. 
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SYMBOLS 


enclosed area, (in.* 

modulus of elasticity, (Ib./in 

modulus of rigidity, (Ib. ‘in 

applied load, (Ib 

coordinate directions 

tangential direction 

normal direction 

displacements in a z directions respectively 
normal and tangential stresses (Ib. /in 
Poisson's ratio 

A ingle of rotation per unit length (rad 

. stress function 

My twisting couple (in.Ibs 


moments and products of inertia 


ars tere = SS ¥ 4A, tive = Sfp 9 44 tay = Spx a 


I PREVIOUS PAPERS it has been shown that a modification 
0 


f the conducting sheet analogy proposed by 


ind Soroka’ and the Gilberts‘ affords a rapid and con 


electrical 
Waner 
venient method for the solution of problems related to plane 
tate of stress and torsion. The advantage of the modified pro 
edure lies in the ease with which electrical potential gradients 

be imposed along a boundary of any complex shape by 
means of a resistance paint as compared to the procedures re 
quiring the establishment of a prescribed current density dis 
In view of this advantage, other 
e<., 
read 


tribution throughout the field 
problems associated with previously proposed analogies 
the Griffith and Taylor analogy®’—can now be solved mort 
This is particularly true for the combined bending and 
beams 


torsion of 


} 


holes) for 


having a multiply-connected cross section 


which the location of the shear center by means of 
the shear flow approximation is no longer valid 
To demonstrate this point, consider a cantilever beam of arbi 


Cartesian coordinate 


trary shape. The origin of the x, y, 3, 

system is taken at the centroid of the fixed end with x and y 
lirected along the principal axes of the cross section. The free 
end, at , is subjected to a twisting couple J; and a load P 
cting at the point x = a, y = 0 in the direction of the y axis.” 
With the svmbols defined above, it is shown in reference 67 that 


differential equations of equilibrium of forces are satisfied by 


eglecting body forces and putting 
0 0 = 7 =z () 1 
o —p _ x /I 2 
A restricted class of problems can be treated for other load distributions 
1 as those disc ed in reference ¢ 
+7 : 
'o put the results appearing here into the notation used in reference 
ve introduce the following transformations; V/ ?, b& G E/201 
a = 6,Wr/I, P/Ir2 


571 


The Editorial Committee does not hold itself 
and 
Od P L 
T = t = 2 — Gy, 
Oy zz \1 +2 ; 
Op 
%y = — t+ GOx 3 
Ox 


where (x, y) is a harmonic function 


The boundary condition for ¢(x, y 


do P ( ss v 
= _- y- 
ds 21 22 i+ ) 


on any bounding curve I, is 


, dy . dx 
t Géy t+ GOx 
d .) ds 


=I 
or alternatively in the Dirichlet form 
. 
P od v r l 
o = a* dy = - y> + - GAx? + y*) + 
lca d Fi O(1 + ») I.“ 2 
1=1,2 , N 5 


where C; is a constant of integration and the curves I, are 
described as positive if the enclosed area is to the left of an ob 
server moving in the direction of description. If the region of 
integration is simply connected then we may let the constant ( 
be zero. If the solid is multiply-connected we can still set C 
equal to zero on the external boundary I’, but we cannot assume 
that C 


for these curves are determined from the condition { 


vanishes on internal boundaries Ir, The constants C, 


do 
ads =U 1= 2, 3, e N (6 
r, adn 


Finally the applied moment, 1/;, can be calculated from the 


uM, = lf (xTy2 — Yrzrz)dA + Pa 7 
JSR : 


which by virtue of expressions 


relation 


3) can be put into the form 


N 
Pe. 2 , ] 


WV, = 2 if Q dA OX ay =, 7 dx T 
JiJR i=1 ry 
—~ ( —— J } 
— LIuus + a «6((8 
ar, (4 » fe ~ Lm) 


The value of AV, is readily computed once the distribution of @ 


GA Izz + Ivy 


has been determined 
In the analog problem, potential gradients required by condi 
a facsimile 


tion (4) can be easily imposed on the boundary of 


conducting sheet model. The function ¢(x, y) is then obtained 


by a simple mapping of the field in the interior of the model 


For multiply-connected regions the constants C; appearing in 


expression (5) can be evaluated as follows. Substitute the series 
representation 
OX y S } Q 
‘ Ln 
\ 
with S ms, = 1 into the requirement (6 his leads to the 
La 1 
matrix set 
~ The method for demonstrating expression (t outlined in the appendix 
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_ | io ("] 
ly, a doy Ms = )0 


(10 
ay; ay ayy he ¥ 
Ww here the coefticients a are expressed as 
de i223, ,N 
a ( ds 1 
I dn s = ], 2, , Nv 


Phe harmonic functions ¢, are the solutions of a series of boundary 
value problems in which C; is zero and the coefficients C; = 6 


(7 = 2,3 N 


are arbitrarily assigned. Thus for each b,, there 


corresponds an a,, given by expression (11). The weighting 
terms m, can be found from Eq. (10) by a simple inversion pro 
vided the square matrix a is nonsingular. Using the values so 
obtained, it can be shown that the coefficients C; are determined 


from the relation 


\ 
Cc, = ¥> mb (Gj =2 3, \ 12 
Las 
As an example, we write for V = 2 
Wy +t Ax Ml = () 
my, +m. = 1 
which vields 
a de 
a = Mle 

ad = ¢ ad» = @ 

and 
a = 0 Co = myba + mobo» 


The solution of the problem of (a) shear center, (b) torsion, 
and (c) bending can be obtained from Eq. (7) by putting @ and 
Vi, P, and AJ; equal to zero, respectively From case to case 
only a slight modification of the boundary conditions is involved 


APPENDIX 


To derive relation (6) in the main body of the note we con 


sider the expressions for displacement as follows: 


P ~ 
= () ~ ) > + po] — cs) (x? — y? — Oy (I) 
2EI 3 ' ? 


w = O(x, y) - 


where ¥(x, y) is a harmonic function Introducing the conju 


gate harmonic function ¢(x, y) as 


GA dy /Orv) = Of/Oy GW Oy/dy) = —(0¢/dx) (IV 
and making use of Hooke’s laws, it can be shown that these ex 
pressions are compatible with expressions (1)-(5) of the main 


text 
By virtue of Eq. (IV) we can write 
dop/dn = —GO(dy/ds) (V) 


Solving for dy /ds from Eq. (III) and using Green’s theorem for 


the plane, we are led to the result (6) since w is single valued 
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Turbulent Boundary Layers on Nozzle Liners 
R.E. Meyer 


Department of Aeronautics, University of Sydney 
February 28, 1955 


a OF THE GROWTH of displacement thickness ay 
usually based on two groups of simplifying assumptio 


Firstly, a Prandtl Number of unity and absence of heat transfer 





at the wall are assumed. Secondly, the velocity profile is 
proximated by a definite power law and the skin friction by som 
empirical relation Ruptash,! for instance, suggests four differ 


ent relations for the skin friction, two of which are 
Cr = 0.0262 | wi /( pimyx 


and 


ér = 0.045 Me ( py lt) 6" . 9 


where suffixes 1 and w indicate values taken in the mainstream 
and at the wall, respectively, and 6’ denotes the total boundary 
layer thickness. Both Eq. (1 
calculate by his method* the growth of displacement thickness or 
rhe thicknesses 
Moreover 


1 used to 


and Eq. (2) have beet 


a certain nozzle liner for Mach Number 2.41 
found at the end of the liner differ by 55 per cent 
Eqs. (1) and (2) lead to different laws of variation of displa 
ment thickness with stagnation pressure! 

The purpose of this note is to draw attention to a simpl 
method of boundary-layer estimation that avoids the assumy 
tions of the second group. If the approximation p/7 const 
is added to the first group of assumptions, the [Illingworth 
Stewartson transformation*®: * would apply if the boundary lay 
was laminar. It is plausible, and experimental evidence tend 
to confirm,’ that this transformation is also relevant to th 
turbulent boundary layer. Assume, then, that it applies 

The stagnation conditions of the correlated incompressibl 
layer may be chosen equal to those of the nozzle and the two 


main streams are then related by? 


and 


ds/ds; = (Ty)/T))81-0/G7-2 = 
1 + (1/2) (y — 142] 37-! 2) (4 


~ 7 


where Ul’ denotes the incompressible mainstream velocity and 
s the length measured along the wall; suffixes 7 and ‘) refer to th 
incompressible layer and stagnation, respectively Che growth 
of the incompressible layer so defined may be estimated with 


confidence. Following Tetervin,®> its momentum thickness 


is given by 


(6,04 +2]9+) = [6;,,U," +2] 9+) + 


PSs 
(n + 1)kv Z UU ds 0 


where a suffix v denotes values at any one fixed station and k and 


n are constants chosen such that 2k( l/@;/v))~” closely approx! 


mates the incompressible skin friction coefficient for the range ol 
, . a° ly 
(U@;/vo) occurring. Since nozzle boundary layers art usually 
thin at the throat, the values of k and n are usually close to ! 


The shape-factor, /7, lies between 1.4 an A 


g 


and 1/4, respectively 
1.6 and is assumed constant. 
The Illingworth-Stewartson transformation also implies the 
relation 
pu dy = pouidy 


and with 6’ 0 at the throat 


* With velocity profile : 
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ess 
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cnesses 
reover, 


splac 


simple 
sump 
const 


worth 


layer 


*SSIDIt 
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READE 


between the distances, y, from the wall, and the relation 

u/u, = u;/l (7 
at corresponding stations. For 
heat 


the velocity profiles, 


Prandt! Number of unity and no transfer, moreover, 


T/T; = 1 + (1/2) (vy — 1)M?7[1 — (u/u 


and hence, the compressible displacement thickness is 
: 7 u 
= )ay 
, T\ l 
R 
1/2) (vy — 1) 47 ' a | 
}5; + (1/2) (vy — 1)M7(0, + 6;){ 
6), Since // 6;/0;, and by Eqs. (4 5 S 
E,6 + (n + I)R( w/e { Bd 
with 
1 + (1/2) (4 1)M2]~ (y+0/@y—2) 9 
}H + (1/2) (vy — 1)M*%(1 + A) 
y [1 + (1/2) (y — 1)APIA—-8V/E7r—2 


7/5, H = 3/2, and n 1/4 


Fig. | shows E and B for y 
Note that 6 is proportional to the [—/(m + 1)]th power of the 
stagnation pressure; when m = 1/4, this agrees with the predic 


tion obtained from Eq. (2) by Ruptash’s method.! 
For comparison, Eq. (9) was also evaluated for the liner men 


tioned above, again with 6 Q at the throat. The results’ 
obtained from Eqs. (1) and (9) are shown in Fig. 2. That ob 
tained from Eq. (2) agreed closely with that from Eq. (9), 


which may be regarded as further evidence in favor of agreement 
tween the predictions of the Illingworth-Stewartson trans 


formation and experimental skin-friction distributions, within 


the present limits of measurement. For fuller comparison 
Fig. 2 shows also the results obtained from Eq. (9) for 10% = 5.2 


nd 3.2 in. at the throat It would appear from these curves 


that the correct value of 6 at the throat cannot differ greatly from 
5.10°%in. in this case 


rhe nearly linear variation of 6 with x, from near the throat 


to the test-rhombus, is due to the variation of B and E with / 


see Fig. 1); for any conventional nozzle they balance so that 

is not far from constant, at any rate up to J/ 2.5, if x is 
easured from the throat This constant, 6(/)//, can be esti 
ited quickly from Eq. (9), with 6, = 0, s» = 0, and s = /, the 


ength of liner from the throat The local rate of growth, dé/dx, 
varies more markedly but it can be shown by the method of refer 
ence 6 that the difference dé/dx — 6(/)// on the supersonic part 


f the liner corresponds to Mach Number variations of only about 
and ambient stag 


All are based on a stagnation pressure of 110 p.s.i-g 


tation temperature 
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0.003 in the test-rhombus, most of which are due to the values 
of this difference not far from the throat 
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Ws: 


governed by the transonic equation 


CONSIDER a compressible two-dimensional inviscid flow 


— — ] a | r¥? T Puy 0 l 


where ¢ is the perturbation velocity potential, y the ratio of the 
specific heats of the gas, a* the critical speed of sound, the sub 
script notation for partial differentiation is used, and the velocity 
directions are given by 


components (U,V) in the positive (x,y 


U=at*-—¢ V=-»9, 2 


Eq. (1) holds under the assumptions 
a ¥y ol ol 
(1 <1 II <1 III < 
a* a* oy Ox 


The transonic Eq. (1) can be put into the normalized form 


2 + & = () 5) 
upon application of the transformations 
x = K-~(2/9E yay nN "Hy + I } 


The title flows arise from an investigation of the possible in 
Eq. (3). To 


consider the 


variant (similarity) solutionst of determine the 


functional form of such solutions, we invariance of 
Eq. (3) under the following continuous one-parameter group of 


transformations: 


T Those solutions of partial differential equations whose functional form 
is such that the number of independent variables in these equations is re 


duced by one. See reference 2 for precise definitions 
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where a (#0) is the numerical parameter and m, n, p (#0) are, 
as yet, undetermined numerical constants. A complete set of 
absolute invariants of Eqs. (5) is given by 


e pim 


ee grey é ¢ (6) 


Now, the normalized transonic Eq. (3) will be constant con 
formally invariantt under the group of transformations (5) if 
m, n, and p satisfy the relation 


n/m = (3/2) — (1/2) (p/m) (7) 
Thus, with g = p/m, Eq. (3) will, by Theorem 3 of reference 2, 


be reduced to an ordinary differential equation by the substitu- 
tions 


eg = Ef(n) n = EA/2D)@—dy (8) 
The resulting ordinary differential equation is 
Igf + (1/2) (q — 3)n(df/dn)] [g(¢ — 1)f + (5/4) (¢ — 1) X 
(¢ — 3)n(df/dn) + (1/4) (q — 3)?n?(d?f/dn?)] + 
(d*f dn?) = () (Y) 
All of the invariant (similarity) solutions of the transonic equa- 
tion under the group (5) are then characterized by Eq. (9) 


T See footnote on page 573 
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In order to determine g, we must impose additional (side 
conditions on the flow. For instance, if we require that th 
velocity component U be unchanged (absolutely invariant 
If, on the other hand, it is required 


] 


under Eqs. (5), then g = 1. 
that ¢ be absolutely invariant under Eqs. (5), theng = 0. Addi 
tionally, the side condition that the transonic equation itself b¢ 
absolutely invariant under the group of transformations (5) im 
poses the restriction g = 3/2. 

Our present investigation, however, will be restricted to that 


class of flows, characterized by Eq. (9), for which g = 0—ie 
¢ = ¢. Under this restriction Eq. (9) becomes 
27? faf an + 45n?f,? — 8fon = 0 (10 


where the 7 subscripts on the f’s denote ordinary differentiation 
with respect to 7. On letting ¢ = f,, and interchanging the 
roles of the independent and dependent variables, we obtain 
from Eq. (10) the following ordinary differential equation for 
nin terms of €: 

45n?¢?(dn/de) + 27n°¢ = 8 (11 
whose solution is found to be given by 

7 = [Cig 9/5 +4. (2 3)¢ 113 (12 
where C; is a constant of integration. 


Further considerations will now be restricted to those flows for 


which C; = 0 in Eq. (12). Since ¢ = f, we then find that 
l : : tie # ; 
eo =f(x) = —-97°+Q=- + (13 
3 Byt+1y? 


where C; is a constant of integration and the expressions for the 
velocity components U’ and V are 


: ] = : > 2 s&s 
Ui = @* 1+ : V=- - (14 
y Fb - syTt1y$ 
In a polar coordinate system, with @ measured with respect to the 
y-axis, such that 
x =rsin@ y = rcosé (15) 


the expressions (14) for the velocity components Ul’ and V be 


come 
2 a 
U = e@*{ 1 + tanto) ; = == tan? @ (16) 
7 Sy +1 
We immediately note that for this flow the sonic line (W = a*) is 


straight and located on the y-axis. Also, the velocities U and V 
are singular on the x-axis. Furthermore, the velocities U and 
V are constant along any one ray issuing from the origin 

The stream lines for the flow defined by Eqs. (16) are deter- 
mined by the differential equation 


dy/dx = V/U i7 


Upon substituting Eqs. (16) in Eq. (17), this differential equa 
tion may be integrated by standard methods to give the following 
expression, in the polar coordinates (15), for the stream lines: 
3 4tan?@+ 3 

exp* tan™! > 

 LaV6ly + (5/8) 2V/6ly + (5/8)] § " 
— cos 6[2 tan? @ + 3 tan? @ + 3(y7 + 1)]"/4 
where C is a constant of integration. Note that r/C is a dimen 
sionless quantity. When @ = 90°, Eq. (18) becomes r 1.180C. 
The stream lines given by Eq. (18), which are symmetrical about 
the line @ = O°, are shown, for C = 0.5, 0.7, and 1.0, in Fig. 1 

The velocities U, V, and W = [U2 + V? and the angle 

a = tan~! V/U for the class of flows given by Eq. (16) are de- 
picted in Fig. 2. It is seen that the assumptions (I), (11), and 
4)” However, for @ < 40° these 


assumptions are not violated and, with appropriate boundary 


(II1) are violated for, say, @ 


layer corrections, the analysis indicates that slightly supersonic 
flows may be turned through small angles. For instance, such a 
channel with entrance and exit sections located at @ = +40 


respectively could be used to turn a flow with w = 1.31a* through 
a total angle of approximately 14 
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One of the solutions of Eq. (9) for g = 1 also yields a flow pat- 
tern with U’ and V + C given by the right hand sides of Eqs 
14), that is, a constant V-velocity component is superimposed 


on the flow field described by Eqs. (14) 
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Simultaneous Measurements of Velocity and 
Temperature Profiles for Adiabatic 
Supersonic Flow of Air in a Tube 


Joseph Kaye,* George A. Brown,t Julius C. Westmoreland, t and 
Edward F. Kurtz, Jr. ** 
Massachusetts Institute of Technology, Cambridge, Mass 


March 1, 1955 


ELOCITY AND TEMPERATURE PROFILES for adiabatic super- 
pian flow of air in a round tube have been obtained simul 
taneously by means of a small combination probe which measures 
both impact pressure and temperature of the air stream from the 
tube wall to the centerline. The need for simultaneous meas- 
urements of velocity and temperature profiles recently became 
evident in a research program whose goal is to obtain reliable 
heat-transfer coefficients for supersonic flow. This program, 
sponsored by the Office of Naval Research, has produced accurate 
and reliable measurements for the rate of heat transfer to a 
supersonic stream of air’? and has also obtained theoretical 
predictions by solution of the equations of momentum, energy, 
and continuity for flow in the entrance region of a tube.* The 
interpretation of the experimental results and the assumptions 
used in the theory have been based, however, on a simple in 
tuitive model used to understand the phenomena occurring in 
supersonic flow in a round tube. This model assumes a laminar 
boundary layer originates at the tube inlet and grows in thickness 
in the direction of flow; coexistent with this laminar boundary 
layer is a central core of fluid which undergoes an isentropic 
process and which decreases in size in the direction of flow 

Previous measurements of velocity profiles' by means of small 
impact pressure probes have demonstrated the existence of this 
laminar boundary layer for a wide range of inlet stagnation condi 
tions, and for both adiabatic and diabatic flow. It was also 
shown‘ that the measured Mach Number profiles of a laminar 
boundary layer for tube flow agree well with similar profiles, 
measured and predicted, for a laminar boundary layer on a flat 
plate with supersonic velocities in the free stream. The present 
measurements of velocity profiles, taken with a new and differ 
ent type of combination probe, confirm quantitatively previous 
velocity measurements. In addition, new experimental evidence 
based on the measured temperature profiles has been obtained 
to confirm independently the existence of the laminar boundary 
layer in the tube flow 

The details of the apparatus used in this work are given in 
previous reports." ? The details of a combination probe are 
shown in Fig. 1 The dimensions of this probe are of the same 
order of magnitude as those used for the simple impact probes.‘ 
lhe thermocouple tip is almost at the inlet plane of the probe 
A new and more accurate traversing rig was constructed for 
positioning of the combination probe, but the method of use is 


essentially the same as described previously The combination 


* Associate Professor of Mechanical Engineering, M.I.1 
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+ Research Engineer, Mechanical Engineering, M.1.T 
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probe was used to measure impact pressure of the air stream from 
tube wall to tube centerline in the usual fashion 

For measurements of temperature, the combination probe 
was calibrated mainly in a subsonic stream of known temperature, 
utilizing a probe recovery factor, defined by 
(T, — T)/(T. — T) (1 


an 


where 7, is the measured probe temperature, 7 and 7) are the 
static and stagnation temperatures of the stream, respectively 
The stagnation conditions were measured upstream of a subsonic 
nozzle, and the static temperature was computed from the 
measured wall pressure at the probe and the condition of isen 
tropic flow for the core of the subsonic stream. The probe was 
also calibrated in supersonic flow in a similar fashion for a Mach 
Number of about 2.3. The results of the calibrations for sub 
sonic and supersonic flow are shown in Fig. 2. A constant value 
of probe recovery factor between 0.83 and 0.85 fits the experi 
mental data, as shown by the solid lines in Fig. 2 

In order to reduce the size of this note, only two representative 
samples of the measured and computed quantities are shown 
in detail in Fig. 3. These two samples were selected to show 
the measured and computed profiles for a low value and for a 
high value of the inlet diameter Reynolds Number, for the 
case of adiabatic supersonic flow in the 0.5 in. tube with a laminar 
boundary layer present. The ratio of the local probe temper 
ature to the local value at the centerline is shown plotted versus 
the radial position, y/a, at the top of Fig. 3. The effect of in 
creasing the inlet diameter Reynolds Number from 35,600 to 
106,000 is to reduce slightly the thickness of the thermal boundary 
layer, without a major change in the shape of the probe-temper 
ature profile in the laminar boundary layer 

The measured static pressure at tube wall and the local impact 
probe pressure were used to compute the local Mach Number, 
employing the Rayleigh pitot-tube equation for supersonic flow 
and the assumption of isentropic deceleration for subsonic flow 
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The local stagnation temperature was computed from the probe 


temperature by 


T,>/To = tp + (1 —1,)T/7 (2 
using a constant value of the probe recovery factor of 0.84 
The local stream velocity was conputed from 
u MV ¢kRT (3 


rhe profiles for velocity ratio and Mach Number ratio are in 
good agreement with those measured previously with a simple 
impact probe These results will not be discussed further 

rhe profiles of static-temperature ratio are also shown in Fig. 38 
rhe effect of increasing the inlet diameter Reynolds Number is 
to decrease but slightly the thickness of the thermal boundary 
layer These profiles indicate the degree of reproducibility ob 
tained by this combination probe as a temperature measuring 
device 

The most interesting profiles shown in Fig. 8 are those for the 
ratio of local stagnation temperature to the centerline value 
In the first attempts to obtain such profiles, the experimental val 
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FIG. 4 - COMPARISON OF MEASURED PROFILES FOR A 
BATIC TUBE FLOW WITH PREDICTED PROFILES FOR 
ADIABATIC PLATE FLOW 


ues were taken at widely separated intervals of the radial positior 
Although some preliminary measurements indicated that ther 
could be a region with an excess of local stagnation temperature 
compared with the centerline stagnation temperature, it was 
decided to eliminate the possibilities of errors in such measure 
ments by repeating them carefully with small intervals of th 
radial position of the probe. The results are shown in Fig 
many more similar results were obtained for other conditions 
The region of this excess of stagnation temperature or enthalpy it 
the laminar boundary layer was studied in detail in order to com 
pare these profiles with theoretical predictions 

Fig. 4 compares the experimental profiles for adiabatic tube 
flow with the theoretical predictions of Van Driest® for super 
sonic flow on a flat plate with a laminar boundary layer. The 
abscissa is the customary one used in the analyses of flat-plate 


flows—namely, 
n (y/L) V Re, ! 


where L is the distance from the leading edge, and the length 
Reynolds Number, Re, is given by 


Re; u..L/v 
where u., and v,, denote the free-stream velocity and kinemati 
viscosity, respectively. For the case of tube flow, however, 


where the centerline velocity and fluid properties vary with th 


distance downstream from tube entrance, a convenient though 
arbitrary definition of 7 is used here—namely, 

n \y BE) * Rez 6 
where L is the distance from entrance plane of the tube and the 


Reynolds Number is arbitrarily chosen for convenience to be 
Ren ul, Li; V ¢ 


in which the fluid properties are to be evaluated at tube entrance 


. or at station 1, and by assum 


(to correspond to flat-plate flow 
ing that an isentropic core exists for the flow at tube entrance 

The experimental profiles for velocity and Mach Number for 
tube flow are in excellent agreement with the predictions for 
plate flow, corroborating similar measurements published pr 


viously. 
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fhe measured stagnation temperature profiles for tube flow 
ell with those predicted for plate flow by Van Driest.é 


ree Wel 


In fact, Fig. 4 shows that the measured excess of stagnation tem 
perature in the laminar boundary layer agrees within a few tenths 
ff one per cent with the excess predicted by Van Driest, based 
iking into account the actual variation of properties of air 
vith temperature 

measurements made with the combination velocity and 


he 
perature probe were used to secure independent checks on the 
ipersonic flow phenomena in a tube, since velocity, temperature, 
ind static pressure are sufficient to define the local behavior of 

flow These checks included the integration of the local flow 

ditions to obtain the mass rate of flow of air, and the integra 
ion of the local energy flux to obtain the total enthalpy flux 
The integrated flow and enthalpy agreed, within experimental! 
or, with the flow and flux based on the calibrations of the 
ipersonic nozzle and on the measured stagnation conditions 
mparison of the defect and excess of stagnation enthalpy flux 
the laminar boundary layer at a given section in the flow is 


so being investigated as an independent check on the reliability 


the combination probe 
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INTRODUCTION 


i: 4 PREVIOUS READERS’ Forum note’ the brachistocronic prol 

has been treated for a constant mass aircraft and for a 
propulsion mechanism operated in such a way that the thrust 
T) depends on the velocity (V) and altitude (/) only 


With the present note the writers generalize their previous co1 


clusions te the case of an aircraft whose mass m = m(ft) is arbi 
trarily time-dependent and whose thrust 7 = 7(h, V, ¢) isa gen 
eric function of velocity, altitude and time (/ A more general 
field of problems is now investigated, insofar as the horizontal 
distance traveled by the ircraft is considered, as it is the case 
with analyses involving the range or with interception problems 
In addition, the variability = o(h) of the acceleration of gravity 


with the iltitude is also accounted for, having in mind possible 
ipplications of the present theory to flight patterns penetrating 
into the upper regions of the atmosphere 

rhe following hypotheses are used: (1) the small angle be 


the vectors thrust and velocity is neglected; (2) only flight 
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paths contained in a vertical plane are considered; (3) the aero 
dynamic lag is disregarded—i.e., lift (1) and drag (D) forces are 
calculated as in steady flight. This last hypothesis is equivalent 
to assuming a drag function of the form D = D(h, V, L), the 
effects of the Mach and Reynolds Numbers on the polar being 
accounted for 

The expressions for the velocity components (#, 4) and the 
equations of the motion on the tangent and on the normal to 
the flight path are written as follows 


] it — Vcosé i) 1 
lh=h—- Vsiag =0 2 

J + gsine + (D—T)/m =0 3 
L — m(ecos@+ V@) =0 | 


where x denotes horizontal distance, 6 is the inclination of the 
flight path with respect to a horizontal plane, and the dot sign 


refers to (total) derivative with respect to the time 


NECESSARY CONDITIONS FOR THE EXISTENCE OF AN EXTREMUM 


The determination of the path of minimum time between a 
be 


set of initial conditions 1 and a set of final conditions 2 te 
reached by the aircraft consists in general in a variational prob 
lem of Mayer® type, insofar as the time to be minimized appear 
in Eqs. (1), (2), (3), and (4 By assuming ¢ as an independent 
variable the brachistocronic problem is formulated as follows 


“to determine the four unknown functions V(t), 6(¢), A(t), a 
hex 


which make dt stationary under the suksidiary conditionst 
e 
(1), (2), and (3 

By introducing the three Lagrangian multipliers A,(‘), Ag( 
\3(¢), the following necessary conditions? are obtained for the 


optimum path 
we 4 Ad Ty — Dy)/m ..D,6 + d sin @ + A, cos 8 0 (5 


d(\3VD,)/dt + X Dy), sin 6 — cos 9 


V(A2 cos 6 \, sin 6 0) 6 
Ako + AT, — Di)/m — d sin 6 D,, cos 0 0) 7 
0 8 
where 
D;, = (OD/dh)y. 1 Dy, = (oD/oV); Dy; OD /oL 
T, = (OT/dh); Ty = (0T/OV), dg/dh 
Eqs. (1), (2 3), (5), (6), (7), and (8) form a differential 


system whose general solution depends on 8 integration con 
stants, to be defined by the boundary conditions of the problem 
It is to be noted that Eqs 5), (6), (7), and (8) are homoge neou 


in A;, Ae, A As a consequence, an arbitrary (provided ~ 0 
value may be imposed to one of the three multipliers—e.g., 
Xx = 1—at a generic point of the flight pattern in view of the fact 
that this circumstance does not affect at all the remaining un 
known functions h(t), V(t), x(t), A(t In addition to the above 
dle condition, seven conditions may be imposed to the flight path, 
for instance, the initial set of values for h, ¢, V, 4, , and the final 
set for h, V, 0, x 

The following general remarks are of importance 

l The above differential system holds also for the extremal] 
of several analogous problems, such as maximum range, max1 
mum increase in altitude, and maximum increase in velocity 
In each case the initial and final set of values for the remaining 
variables must be given 

2) With the help of convenient boundary conditions 


versality conditions) the above equations may also be used to 


+ If the lift L is eliminated from I 3) by using 1. (4 ne obtair 
an expression of the type J é £ 0 
t Eqs. (5) to (8) may be formally deduced as Euler juat r 
— 3 
¥ (5 n. J.) dt 
1 
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solve the interception problem—i.e., to determine the brachisto- 
cronic path for an aircraft A which has to reach another aircraft 
B, performing an arbitrary maneuver x,(¢), A#,(t) in the same 
vertical plane. In addition, also problems of conditioned optima 
may be treated with the above equations, for instance the im 
portant case (h 7 h,) where the aircraft is forbidden to go to alti 
tudes below the initial one, because of the presence of obstacles, 
such as the ground level 

(3) Eq. (1) could be replaced with an equivalent isoperimetric 
*2 


condition by prescribing the value for J = | Vcos@dt. Insuch 
J 1 


would depend on 
seven integration constants V,0 
prescribed at 1 and (h, V, 6) at 2, in addition to the idle condition 
(1), however, has been preferably 


i case the solution of the differential system 
For instance, (h, ¢, may be 
and to the value for J. Eq 
used by the writers, with the main object of treating the actual 
variational problem within the general frame of Mayer problems 


SPECIAL CASES 


The integration of the above differential system is to be per- 
formed by approximated methods in the most general case 
Under particular circumstances, however, the problem may be 
simplified. 

(A) 

If no condition is imposed to the values of x; and x2, Eqs. (1) 
and (8) disappear. The new set of Euler’s equations is obtained 
as a particular case of Eqs. (5), (6), (7), and (8) for \; = 0. The 
number of integration constants reduces to 6, with an idle condi- 
tion at an arbitrary point of the flight path. The values for h, 
V, 6, t at point 1 and for h, V, 6 at point 2 may be prescribed 
Should the Lagrangian multipliers be eliminated, the order of the 
resulting system would be reduced to 5, the idle condition disap- 


Free Horizontal Distance 


pearing. 
(B) Constant Mass, Thrust Independent of the Time 

The hypotheses m = 0, 07 /Oot = O lead to a first integral for 
the differential system constituted by Eqs. (5), (6), (7), and (8): 


A V (9) 


+ D,V@) + V(r. sin 6 + A; cos 8) = C 


By replacing one of the above equations with Eq. (9), a new 
system, homogeneous in )\;, Ax, Az; and C, is obtained. It follows 
that an alternative form of the idle condition could be C = 1. 

If the time is eliminated by using Eq. (2), the number of inte- 
gration constants for the differential system reduces to seven 
After accounting for the idle condition, six other conditions may 
be imposed—e.g., by prescribing the values for h, V, 6, x at 
points 1 and 2 

For the particular case T = D = O the order of differential 
system obtained by eliminating the time variable and the La 
grange multipliers reduces to three and the solution degenerates 
into the well-known brachistocronic of the gravitational field 
i.e., the cycloidal path. 

(C) Free Horizontal Distance, Constant Mass, Thrust Independent 
of the Time 

For \; = 0, m = O, O07 /ot = O, the elimination of the time 
variable and of the Lagrange multipliers leads to a fourth order 
problem, which has been already treated in a previous Readers’ 
Forum Note‘ by the same authors. 

(D) Free Horizontal Distance, Negligible Centrifugal Forces 

For L = 
as the terms \3D,6 and d(\3D,V)/dt disappear. 
moreover, the elimination of the Lagrange multipliers \. and ) 
from Eqs. (5), (6), and (7) leads to the following result: 


mg cos 6 Euler’s equations (5) and (6) modify insofar 
For \; = 0, 


Wi + (Dy — Ty) gu + We using 4 
D, — T), + mM D, g’ cos@é@ = 0 (10) 
where 
u = (Drztan6 — 1)/V (11) 


Eqs. (2), (3) and (10) form a differential system of order three, 
for which the initial set (h, V, ¢) and the final set (h, V) may be 
a parabolic polar with 


prescribed. It is to be noted that for 
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coefficients depending on both the Mach and Reynolds Numbers 


the solution cos 6 = O is found in addition to the solution (10 
(E) Free Horizontal Distance, Required Lift Independent of 


Centripeta! Acceleration and Path Inclination 


aod 
‘ 


For \, = 0, L = mg, Eqs. (5), (6), and (7) lead to a solutio; 
which is composed of arcs of equation: 
cos 8 = UV 9 
and ares of equation: 
T—- D+ (Di — T,)V?/2 + (Ty — Dy)V 4 
V2Dy g'n 
Eq. (18) may also be written as 
0/oV) (TV — DV) —(V (0/oh)(T — D) =0 (14 
The above expression is a relationship of the type f(h, V, /) = 
For the particular case m = 0, O7/Ot = 0, Eq. (14) degenerates 


into an expression in finite terms between 4 and V (see references 
1, 2, 4). 
(F) Free Horizontal Distance, Negligible Induced Drag 

For \, = 0, Di = 0, Eqs. (12), (13), and (14) hold with obvious 
Do, where Do = 
0 Eq. (14) 
finite terms, as shown by the junior 


simplifications (D = zero-lift drag For the 


particular case 07'/Ol = reduces to an equation ir 


iuthor in reference 5 
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Note on a Formula of H. W. Emmons 


J. A. Steketee 

Department of Mathematics, University of Toronto, Ontario, Canada 

March 9, 1955 

I THE PAPER “The Laminar-Turbulent Transition in a Boun 
dary Layer’! Prof. H. W theory of 

starting 


Emmons develops a 


boundary-layer transition from the observation of 


“turbulence spots’’ which grow while they are washed away by 
the flow. 

Fundamental for the applications of the theory is the formula 
(4) which gives an expression for the probability that a point on 
the watertable at a certain moment will be covered with turbu 
flow. If the mean flow is steady one is justified in inter 


preting (4) as the expression for the fraction of the time that a 
point on the watertable will be 


lent 


covered with turbulent flow 


was felt to be rather compl! 
cated through the introduction of the (PP) d Vo functions which 
It is the purpose of this 


The derivation of the formula (4 


were introduced to eliminate ‘‘overlaps.”’ 
which is believed to be 


note to give a derivation of formula (4 
for the W(PP)) dl 


more straightforward as there is no need 
functions. 

We first consider the following case of probability theory 
there are given simple, independent events. The probability 
that the kth event occurs is px, while the probability that the &th 
event does not occur is 1 — p,. The p 
(If all the p,’s are equal we have the case of ” Bernoulli 
ind we 


’s are constant and un- 


equal. 


trials The simple events are not mutually exclusive 


occu 
of th 
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be tl 
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isk for the probability that at least one of the first k events 
occurs There is no loss of generality in asking for at least one 
of the first & events as it will always be possible to renumber the 
events in such a way that the events we wish to consider will 
be the first & in the new numbering system Let this probability 


be denoted as fi, 4 We have then immediately that 


D1, 2 = pi + po — pipe = 1 — (1 — pri) (1 —p 


= pit Pur Pipue =1—-(1- ps 


In particular for k = n 
bo2=1l- Tai-p (2 
1 l 

This answers our question but to relate it to formula (4) of 
Emmons’ paper we proceed as follows: we suppose that a vol 
l i ' 
volume of the 7th cell will be denoted by AV;, and P; wili be a 
Each event of our 


V is given and let this volume be divided in » cells The 





point in this cell and representing this cell 
problem can then be associated with a cell in the volume V 
The probability that the 7th event occurs becomes in this new 
situation the prob ibility that the event in the 7th cell occurs and 
If we put 


we write g( P; )AV; instead of p 
B, = If (1 — pi 


this expression becomes in the new notation 


j a Bi ae 
B, = 1 — g(P;)Al = exp > In [1 — g(PPAI 
( | f 


j of »  [g( PAV 


By =exp 2 — 2 


ind in particular 


n | (PAV | 
Pi n - C&zp >» — >> 3 
t=1 ( m l m { 


We suppose now that the number of cells in V increases indefi 

nitely in such a way that each of the AV;’s becomes smaller and 

smaller. The second and higher terms in the series for the 

logarithm can then be neglected compared to the first term as 

they are of second or higher order of smallness in the infinitesimal 
n 

AV;, while the S becomes an integral over V’ 

i=! 

When we further associate a point ? with V and write f(?) in 


stead of ¢ , » for the probability that at least one event occurs in 
V, we finally obtain the relation 
7 
Py=i—exp}— | g( P. avi } 
Jl 
which is formula (4) in the paper of Emmons. The volume V 
associated with P is the volume of the backward cone of P as 
described in the paper 
It is now easy to construct the analogues of the other formulas 
obtained by Emmons. When we introduce v,, defined by 


Vm = Pi, m+1 — Pry 


it follows immediately from Eq. (1) that 


mm = parr I (1 — (5 
i=1 


Further we rewrite Eq. (1) in the form 
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m 


The v,, corresponds to the ¥( PP») dVo while our Eq. (5) ts easily 
reduced to a form similar to (A13) if we use the expression de 
veloped for B Finally our Eq. (6) corresponds to the integral 
equation (2) in the Emmons paper 
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A Simplified Expression for the Dihedral Effect 
of a Flexible Wing 


W. P. Rodden 

Assistant Engineer, Missiles Division, The RAND Corporation 
Santa Monica, Calif 

March 14, 1955 


SYMBOLS 


a local angle of attack, rad 
B angle of sideslip, rad 
I loca! dihedral angle, rad 


wing spanwise coordinate, ft 


Zz wing deflection, ft 


b wing span, ft 
I wing lift, Ibs 
ui airplane gross weight, Ibs 


airplane normal load factor 
wing tip spring constant, ft./Ib 


f rolling angular velocity, rad. per se« 
| airplane velocity, ft. per sec 

( rolling moment coefficient 

( 0C1/O(pb/27 

Cig oC 1/08 


Subscripts 


rigid 
flexible 


wing tip 


DISCUSSION 


B MAKING TWO REASONABLE APPROXIMATIONS a variation of 
sideslip rolling moment coefficient with load factor can be 
obtained for an airplane The first approximation is to replace 
the wing normal load factor by the airplane normal load factor 
the second is to assume the wing deflection mode to be parabolic 

The local angle of attack of a sideslipping wing with variable 


dihedral is 


where Fr = Tl, + 
9 


= 5, (2y/b)? ys 
and the tip deflection, which is proportional to the wing lift 
can be written 
=kL, =kWn 
if k is defined as a spring constant depending on the wing rigidity 


and spanwise lift distribution 
The variation in dihedral becomes 


Ty = dz/dy = 82 y/b?) = S(RW/b?)n,) 3 
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The deviation of the local angle of attack from the rigid con and 


dition then is 
, ae Cr, = (Cig)o + AC 
Aa = BIy = 8(RW/b?)n By (4) 8 3 1g S 
This is o » same form as the spanwise angle of attack vari viplece . ‘ ; . 
4 a is of the same for ihe ‘ 8 rhis expression shows the dependence of the dihedral effect 
atio where é : 
ation in roll where the gross weight and the normal load factor. It allows a simp}, 
a = py/V (5 correction to wind-tunnel data from rigid models and also offer 


a - a convenient estimate of dihedral effect in the preliminary desig 
‘oO $s Teaso 


stage since spring constants can be readily calculated from estj 
a oc _ oe a A 3, (6 mated rigidities and the Schrenk approximation to the lift di 
O( pb/2V) O|(4kW/b)nB tkWn,; P tribution and variations of C;, with plan form and Mach Number 
a are giv en in the literature The author has obtained good cor 
relation of this theory with some Douglas El Segundo low-speed 
ACi, = (4k1N/b)n.Cip (7 wind-tunnel tests of a flexible model 
. 
Some Static Aeroelastic Design Considerations span. However, where the control surface is actuated at one end 
Relating to Flap-Type Control Surfaces of a torsion bar, the assumption of constant 6 is invalid and could, 
of course, yield completely misleading results. In the following, 
some qualitative results from a series of recent aeroelastic caleu- 
Martin Zlotnick lations are discussed, and the condensed derivation of the aero 
Aerodynamicist, Aerophysics Section, Republic Aviation elastic equation for the independently flexible control surface is 
Corporation ‘ pct : é ; cipgte ; 
March 10, 1955 given rhe latter is quite straightforward and is included o 
account of its previous omission in the literature 
The aeroelastic equation is derived from the relations between 
PF  gasnecitg METHODS FOR CALCULATING the aeroelastic loading (1) the applied aerodynamic load and the resulting streamwise 
of flap-type control surfaces have appeared in the liter change in angular displacement due to flexibility, (2) the stream- 
ature, the author is not aware of any published paper that con wise angular displacement and the resulting aerodynamic load 
siders the flapped surface and the main surface to be independ and (3) the angular displacement at equilibrium and the initial 
ently flexible. The assumption usually made in aeroelastic an displacement 
alysis is that the angle between the flap surface and the main If as and 6, are the streamwise angular displacements on the 
surface (the angle of deflection, 5) remains constant spanwise main surface and flap, respectively, due to the equilibrium aero 
Such an assumption is completely valid where the control surface dynamic load while /,,/qg and /;/q are the equilibrium loadings on 
actuators are not sensitive to load and are distributed along the these surfaces, then relations (1) are 
— pte! 4 ag Pm X, + oS ly al ‘| 
qd qd ( 
i we oy t oS ry fm + gSj) I 5 osm! I 
g qd g q 


where the superscript I denotes a loading due to a and II a loading due to 6, where S denotes the appropriate matrix of structural influence 
coefficients, and where the first subscript indicates the affected surface and the second the affecting surface. The subscript m refers to 
the main surface and / refers to the flap 


Relations (2) are 


q q 
' 


where A denotes the appropriate matrix of aerodynamic influence coefficients and where a and 6 represent the equilibrium angular dis 
placements of the main surface and flap, respectively. Relations (3) are simply 


where ay and 6) are the initial angular displacements 
These relations can be combined and expressed as a pair of simultaneous matrix equations or as a single partitioned matrix equation 


for a and 6 in terms of ag and do as follows: 


[7] Sia Aen Tt Sut Ate SmmlA ms + Smg"Ay fal j a { { 
* LS'Awn + Sy'hrm — 1! Sa Aay + SpA, taf VSo + aw! 
where / denotes the unit matrix 

The structural influence coefficients giving the flap deflection due to a load applied to the main surface (the S;,, matrices) will not always 
be available from the usual stress analysis. A valid assumption for most practical configurations where the main surface is rigid com- 
pared to the flap and the loads are applied near the elastic axis of the main surface is that 5, is due only to aerodynamic torsion loads on 
the flap. This restriction implies that the flap will act as a rigid extension to the main surface when the main surface is under load 
Under these conditions 
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t is, the flap deflection due to loads on the main surface is equal to the main surface deflection 


[7 Smm'A mm + Smg'A 


vhere the primes on the S;,; matrices indicate that they represent 
flap deflections due only to the aerodynamic torsion loads 

It might be expected that if the independent flexibility of the 
flap were neglected, the total load on both surfaces would be un 
fected because of compensating effects—that is, neglecting the 
flap torsion would impose more load on the main surface and 
result in a greater loss in lift on the main surface but none on the 


flap, 


whereas including the flap torsion would give a greater 
loss in lift on the flap but a smaller loss on the main surface 
Certainly this could be the case for some configurations, although 
the correct estimate of the distribution of load would depend on 
the correct assumptions in all cases 

However, the calculations made for practical configurations do 
not always show this trend. It was found that the effectiveness 
factor due to 6, ny; = (Cys) fexible/(Cy, rigid aS Computed on the 
assumption of independent flexibility was on the order of half 
that computed on the basis of constant 6 for some configurations 
The results based on the constant 6 computation may lead to the 
conclusion that the control is satisfactory when in fact it is 
not. It may be noted here, however, that the parameter ny, 
was practically unaffected by the flexibility of the flap 

It would be worth while to make a series of calculations to 
determine if there is a ratio of (GJ) of the flap to (GJ) of the main 
surface when the independently flexible flap would give about 
the same effectiveness as the constant 6 flap for a given range of 
dynamic pressure. This would indicate the circumstances under 
which distributed actuators could be dispensed with in favor of 
the simpler design involving a single actuator and a torsion bar 
\ single actuator with a torsion bar might also be favorable if the 

iin object is to delay aeroelastic reversal at high dynamic pres 
sures where only relatively small effectiveness is required. It is 
possible, however, that the most favorable design where high 
aeroelastic effectiveness is a critical factor will be with actuators 
distributed spanwise 

In considering the effectiveness of an independently flexible 
flap-type control surface, the parameter, ny;, is of primary im 
portance as a “figure of merit’’ only when the flap is actuated by 
position-sensitive servos. If the airplane is manually controlled 
or “hinge-moment sensitive’’ servos are used, the appropriate 
figure of merit’? would be n;; = ny;/ny5 Where ny; is the hinge 
moment ratio corresponding to n,,. In either case, the factor 
ny, is the parameter which corresponds to the normal force 


available for a given hinge moment and is an aeroelastic param 


eter of primary importance for the control system design 


On Approximate Solutions for the Compressible 
Flow Around Oscillating Thin Wings 


Marten Landahl, Holt Ashley, and Erik Mollo-Christensen 
Massachusetts Institute of Technology, Cambridge, Mass 
March 14,1955 


 pegeg A THIN WING with root semichord ) = 1 oscillating 
in a free stream of velocity U = 1 with reduced frequency 
of k = wh/U, w being the circular velocity. The velocity of 
sound is then 1/7 and, assuming linearization, the perturbation 


velocity potential of the flow, (x, y, 2, ¢) = o(x, vy, ze, follows 


FORUM dS 1 


Eq. + may then be rewritten 


the equation 


The boundary condition to be fulfilled on the wing is 
g(x, y, 0) = hil(x, y) + tkh(x, vy) = W(x, 2 


where the position of the wing mean surface at any instant is 
defined by 

= h(x, ye 3 
The pressure difference across the wing is given by 


Ap 
t = —4do, + ike)e! } 


a 
{ 


Exact solutions of Eq. (1) can only be obtained in special cases 


such as the two-dimensional supersonic case. However, by 
neglecting one or several terms in Eq. (1), simple solutions can 
be obtained. Three particular approximations will be discussed 


here, applicable for different ranges of k and .\/ as indicated be 


low: 
— M*¢g,, + gy + ¢: 2kill?g, + k?M¢ 0 5 
for 
VW? > 1 and/or k | 
Puy + O22 — ZkiM*¢ k?A*¢ 0 6 
for 
k>(M—- 1 
Guy + oz + RAMP = UV fork > 1 7 


It is also necessary in all these cases to impose the condition 
(1 — M?) = 0(M 

The validity of Eq. (5) for 17? > 1 is immediately evident by 
comparison with Eq. (1); the alternative condition for large 
values of & and supersonic flow was established by Hjelte! using 


the Laplace transformation. Simple solutions of Eq. (5) are 


possible to obtain by introducing 


giving 
— Mz. + ¥ t ¥:: = UV (9 


Lomax, Heaslet, and Fuller? also used a corresponding approxi 
mation for treating indicial problems 
Eq. (6) has been applied for a small aspect ratio delta wing in 
reference 3. Eq. (7) was treated in reference 4 
An easy way to estimate the errors involved in the different 
approximations is to calculate the two-dimensional case where 
») 


the ¢,, term drops out; it is then readily shown that Eqs. (2), (3 


(8) and (9) lead to the remarkably simple solution 
Ap/q = —[4 wa Wie 10 


The result agrees with that of Lighthill® for a wing of zero thick 
ness and constitutes, in fact, the linearized counterpart of the 
‘piston theory,’’ numerical results of which were given by two 
of the authors.6 Although Table 1 of reference 6 is correct, the 
derivation of ‘‘piston theory” given there is superseded by the 
one in the present paper, which is more satisfactory in every 
respect and involves less restrictions on the sizes of the param- 
eters 

The solution of Eq. (6) can directly be obtained from the f = 1 


solution by means of a simple transformation 
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Eq. (7) gives the one-dimensional acoustic theory. Retaining 
for comparison all terms in the boundary condition, Eq. (2), and 
the pressure coefficient, Eq. (4), the result reads 

Ap — wr + tk ‘ict AY 
7 ikM 
For high reduced frequencies the result thus approaches that of 
Eq. (10). 

To illustrate the ranges of validity of the various theories, the 
coefficients L;’ and L,’ (defined by Garrick and Rubinow’ as the 
real and imaginary parts of the dimensionless lift on an airfoil 
oscillating about its leading edge) have been calculated according 
to the above results and are plotted in Figs. 1 and 2. 

It is seen from the graphs that for flutter, where usually 1/k = 
0(1), the transonic approximation could probably be used in the 
range, M = 0.9 — 1.25. The piston theory seems adequate 
for M > 2.5. For high frequency flutter, any of the theories 
can be used. The one-dimensional acoustic theory, however, 
requires the frequency to be very high. Eq. (7) leads to much 
more accurate solutions in the three-dimensional case when the 
aspect ratio is sufficiently small (cf. reference 4). 

As an amplification of the numerical results presented in refer 
ence 6, it is worth noting that ‘‘piston theory” formulas like Eq 
(10) are proving most useful in dealing with complicated high- 
speed aeroelastic problems. For example, the supersonic bend- 
ing-torsion flutter calculation performed by Garrick and Rubinow’ 
with considerable numerical labor reduces, in the most general 
case, to evaluation of a single radical. In the special circum- 
stances when the bending-torsion frequency ratio is small and 
the wing center of gravity does not lie too close to the mid- 
chord line, the expression for reduced flutter speed reads 
Up (Xqa + a)? + (1/3) 


~V uM,’ 
—T M. a [ °° 9 ° « 
bwe \ (Xa +a) (Xe + a)? + rea — x7q + (1/3) 


(The notation here is that of ref. 7 and is familiar to aeroelasti 
cians.) Calculations by Eq. (12) reproduce curves like those in 
Figs. 17 and 18 of reference 7 to an accuracy of better than 10 per 
cent above M = 2.5. The asymptotic behavior for large values 


of the Mach Number is, of course, exact 


REFERENCES 


1 Hjelte, F., Methods for Calculating Pressure Distributions on Oscillating 
Wings of Delta Type at Supersonic and Transonic Flow, Swedish K.T.H 











aie weg 
8 
mu: & 
4 
~ 
\ 
mu: 12 
5 
M: 2 
0 J 
t@) | 2 3 
1 
k 
Fic. 2 


AERO T._N. 39 (to be published 

2? Lomax, H., Heaslet, M. A., Fuller, F. B., and Sluder, | Tu und Th 
Dimensional Unsteady Lift Problems in High-Speed Flight, NACA Rept 
1077, 1952 

Landahl, M. T., The Flow Around Oscillating Low Aspect Ratio Wing 

Transonic Speeds, Swedish K.T.H. AERO T.N. 40, 1954 

4 Merbt, H., and Landahl, M. T., Aerodynamic Forces on Oscillating L 
Aspect Ratio Wings in Compressible Flow, Swedish K.T.H. AERO T.N. 3( 
1953 

5 Lighthill, M. J., Oscillating Airfoils at High Mach Number, Journal of the 
Aeronautical Sciences, Vol. 20, No. 6, pp. 402-407, June, 1953 

6 Mollo-Christensen, E., and Ashley, H., Applicability of Piston Theor 
the Flow Around Wings in Unsteady Motion, Journal of the Aeronautica 
Sciences, Vol. 27, No. 11, pp. 779-780, November, 1954 

7 Garrick, I. E., and Rubinow, S. J., Flutter and Oscillating Air-l 
Calculations for an Airfoil in Two-Dimens al Supersonic Flow, NACA 
Rept. 846, 1946 


8 Nelson, H. C., and Berman, J. H., Calculations on the Forces and M 
ments for an Oscillating Wing-Aileron Combination in Two-Dimensional 
Potential Flow at Sonic Speed, NACA TN 2590, 1952 

® Anonymous, Tables of Aerodynamic Coe ficients for an Oscillating Wing 


Flap System in a Subsonic Compressible Flow, N.L.L., Amsterdam, Rept. I 
151, 1954 


A Remark on ‘“‘A Simple Laminar Boundary 
Layer with Secondary Flow’’ 


M. Z. v. Krzywoblocki 

Professor, Department of Aeronautical Engineering, University 
of Illinois, Urbana, III 

March 15, 1955 


T° REFERENCE 1 Loos talks about an exact solution of th 
boundary-layer equations. This probably refers to the 
Blasius solution as well, which is cited in the text. Let meé 
mention that with the usually assumed boundary conditions wé 
cannot talk in the present state of available mathematical tools 
about exact solutions of the boundary-layer equations. It 1s 
possible that there exists no solution at all and that the so called 
solutions of the boundary-layer equations represent something 
which cannot be judged by mathematics in its present state, OF 
that there exist more than one solution, perhaps even the W hole 
family of solutions. Moreover, H. Weyl, Oudart and ‘ )strowsk} 
have shown that the Blasius solution breaks down at a certaim 
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point. All these things are too complicated to be discussed in 
: engineering journal, but I may refer the reader to my papers 
» that subject published in recent years in various mathematical 
scientific journals 


REFERENCE 


with Secondary Flou 


H.G 1 Simple Laminar Boundary Layei 


rnal of the Aeronautical Sciences, Vol. 22, No. 1, pp. 35-40, January 


A Short Cut in Computing Product of Inertia 


Kennetn 

Professor of Aeromechanics 
Worcester, Mass 

March 17, 1955 


G. Merriam 
Worcester Polytechnic Institute, 


I’ PUBLISHED PROPERTIES Of structural sections, it is customary 
to include values of J,;, Jy, and Jmin, from which one can com- 
using the well-known relationship 


Reis * Eas * De OR 


pute ra Ww 


From the geometric properties of the inertia circle, it is not 


dificult to show that 


(Izy)? - Il, = I mar! min 


Using the two relationships, the calculation of /;,, from known 


values of J,, Jy, and min, becomes a simple operation 


Interaction Equation for the Buckling of 
Unstiffened Cylinders Under Combined 
Bending, Torsion, and Internal Pressure 


Bertram Klein 
Research Engineer, Los Angeles, Calif 
March 16, 1955 


| pes ENTLY IT IS NECESSARY to determine margins of safety 
for unstiffened cylinders under combined loadings in the 
presence of internal pressure. The cases of compression and 
internal pressure! and torsion and internal pressure? have been 
considered. Also previously treated is the case of bending and 
torsion without pressure.’ It is the purpose of this note to show 
how to treat the case of simultaneous action of bending, torsion, 
ind internal pressure 


The proposed interaction equation for buckling is 


Ry? + Rs? = Rp(1 — Ro?) + 1 
where 
Ry = o/(e¢r); R, = t/( ter Rp = p/P. 
nd 
0 applied extreme fiber compressive stress including 
effect of internal pressure 
%r)p = allowable extreme fiber compressive stress for bend- 
ing in the presence of internal pressure 
T ipplied shear stress 
Ter illowable shear stress without pressure 
pb = internal pressure 
Per = allowable external hydrostatic pressure 
rhis interaction equation reduces to the above mentioned re- 


sults when the proper quantities are set equal to zero. Further- 
more, the shapes of the interaction curves plotted in Fig. 1 are 
what one expects for this type of problem 

In conclusion, it must be remembered to check for possible 


lracture under combined stress at the extreme tensile fibers of 


on 
ne 


FORUM 








°0 a2 04 06 O8pl0 12 14 16 
A 


Fic. 1 Interaction curves for the buckling of unstiffened 
cylinders under combined bending, torsion, and internal pressure 


the cylinder whenever internal pressure is present 
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ure on the Buckling of Thin-walled Circular Cylinders under Torsion, NACA 
ARR NO. L4 E 27, May, 1944 

’ Bruhn, E. F., Tests on Thin-walled Celluloid Cylinders to Determine 
Interaction Curves under Combined Bending, Torsior Compression of 





Tension Loads, NACA TN 951, 1945 


The Effect of Axial Flow in a Burner on 
Combustion Oscillation Frequencies 


R. T. Newton and J. C. Truman* 

Small Aircraft Engine Department, General Electric Company, 
Lynn, Mass., and Aeronautical Engineering Department, USAF 
Institute of Technology, Wright-Patterson Air Force Base, Ohio 
Respectively 

March 21, 1955 


Cn TYPES OF COMBUSTION oscillations are found to have 
frequencies equal to resonant frequencies of the gas in the 
combustion chamber. It was shown in a previous note! that 


the effect of an axial flow through the chamber is to decrease the 


longitudinal resonant frequency by a factor (1 — 17?); that is 
f, = f,(1 — M* (] 
where 
f, = longitudinal resonant frequency of the gas in the 
chamber, cycles per sec 
f,. = longitudinal resonant frequency of the gas in the chamber 


for J = Uv. cy cles per sec 
WV = average Mach Number of the axial flow 


To examine the effect of axial flow on transverse resonant 
frequencies, a result obtained by Woolston and Runyan? may be 
employed In considering a wing oscillating in a subsonic com- 
pressible flow of Mach Number 


section, they have shown that the transverse resonant frequencies 


Vin a tunnel of circular test 


ire given by 


fr = (Camn/#D) V1 VW 2 
where 
1 transverse resonant frequency of the gas in the tunnel 
or chamber ), cevcles per see 
( = velocity of sound 
D = tunnel (or chamber ) diameter 
Qmn = Characteristic numbers for Bessel functions 
* Now, Small Aircraft Engine Department, General Electric Company 


Lynn, Mass 








584 JOURNAL OF THE 


Eq. (2) may be written in a form analogous to that of Eq. (1 


as 


fy =fnmv 1— M (3) 


where, again, the subscript zero denotes the condition JJ = 0 
Eq. (3) also follows if one considers the axial and transverse veloc 


AERONAUTICAL 


ity components of a pressure wave transmitted by a signal source 


located on the axis. This approach is noted by Runyan and 


Watkins 


These results are based on the assumption that the Mach 


Number of the axial flow is constant. In most combustion 
systems, Mach Number varies significantly throughout the 
burner, both along the axis and in a transverse plane. The 


accuracy of these results, applied to a given burner, may there- 
fore be questionable, although, qualitatively, one would expect 


the predicted reductions in frequencies to occur 
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Some Considerations on the 
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Note on the Compressible Laminar Boundary 
Layer with Heat Transfer and Pressure 
Gradient 


Eli Reshotko and Clarence B. Cohen 
Lewis Flight Propulsion Laboratory, NACA, Cleveland, Ohio 
March 16, 1955 


SYMBOLS 


ratio of velocity within the boundary layer to local external 


velocity 
(0 wall temperature gradient parameter 
V/ Mach Number 
Vu Nusselt Number, Vu v(OT/ OY) w/ Tau 1 
Rey Reynolds Number based on wall fluid properties 
Ta adiabatic wall temperature 
Te stagnation temperature 
Ts, free-stream stagnation temperature 
1 wall temperature 
3 pressure gradient parameter 
n boundary layer similarity variable 
o Prandtl Number 


DISCUSSION 


I EVY!' POINTS OUT ON THE BASIS of approximate calculations 
that for a Prandtl Number (oc) different than one, there is 
—e'(0 


Mach Number effect on the heat-transfer parameter 
tie/tin) = 4 


even 


However, Levy’s approximate calculations 


seem to be qualitatively inaccurate. Reproduced here 
(Fig. 1 
a number of exact solutions to Eqs. (14) and (18) of reference 
1 obtained by the writers using an IBM Card-Programmed Cal 


The fact that the quantity g’(0)/[(7./T.,) — 


is Fig. 18 of Levy’s paper on which are also indicated 


culating machine 
1] indicated by Levy for 1] ~ o, o = 0.7 is finite through the 
value T,,/T;, = 
solutions obtained by the present writers for 7,/7., = 1 show 
that ¢/(0) ¥ O for this case, thus ¢/(0)/|(7,./7,,) — 1] 
However, for JJ ~ <, ¢’(0 
M = V5, g'(O) is zero when 7,,/7,, = 0.90, corresponding in 


1 requires a recovery factor of one. The exact 


is infinite 


is zero when 7),./7,, = 0.80 and for 


both cases to recovery factors of 0.80. For cases where Ty, + 


T,,(0 # 1, 7 # O) the heat transfer quantity related to Nu- 
V Rew» would be [—2'(0)]/[(T/Ts,) — (Taw/Ts,)|. The present 
solutions for \J = and WM = 7/5 as well as solutions for 


SCIENCES AUGUST, 1955 
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Fic. 1 Effect of Mach Number on shear and heat-transfer 
functions for 8 = 1, Prandtl Number = ().7 


M = Oare plotted in Fig. 2 in terms of this suggested parameter 


It is seen from Fig. 2 that the additional term in Eq. (18) of 


reference | contributes only slightly to the heat-transfer paran 


eter when referred to the adiabatic wall temperature. Further 
more, the values for JJ = 0, o = 0.7 are easily obtained from 
Prandtl Number one solutions (arbitrary Mach Number 


discussed in reference 2. Thus the correction for Prandtl Nun 
ber which is used for low-speed flows appears to be adequate for 
high-speed flows 

In reference 1, Levy has briefly commented on the uniqueness 
of the solutions. Some further details may be of general interest 
Eqs. (14) and (17) of reference 1 have been solved by digital 
calculation for o = 1 by the writers and the results have been 
tabulated.’ For adverse pressure gradients a lack of uniqueness 
has been observed Previously, Hartree‘* also observed for co 
stant fluid properties that the boundary conditions: 
f’'(O = U, ri = | are not sufficient to determine a un 
solution in the region 6 0. That is, there are an infinite nu 
ber of values of the shear function at the wall, , for a give 
pressure gradient parameter, 8. Some of the properties of th 


solutions can be related to the coefficients in the following « 


pression for velocity ratio valid for large n with P? l 
f’ = 1+ lal(n—k 23+ 1 + (a,/2 y= 4 
exp } —i(m — x)?/2 { T- a(n K)~ 


Here f’ is the ratio of velocity within the boundary layer to th 
local external velocity, ai, a2 and « are arbitrary integration cou- 
stants and a; is an integration constant related to the distributio! 


of stagnation enthalpy in the asymptotic region through the re 


lation 


9 


i T;, = 1 + as(n — x) exp | —[(n — «)*?/2]} - 


For 8 > 0 (favorable pressure gradient) a2 is necessarily zero m 
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order to satisfy the boundary condition f’(=) = 1. A unique 
solution results under these circumstances. For 8 < 0, a is 
chosen as zero to avoid infinite values of displacement thickness 
related to the integral f, (1 — f’)dn) 
the equations for adverse pressure gradients with a: = 0, two solu 
tions with finite displacement thickness still remain. These 


However, in solving 


correspond to the maximum and minimum values of wall shear 
obtainable without incurring velocities within the boundary 
layer in excess of the free-stream velocity. For Hartree’s 
problem (7/7, = 1), 7s/Ts, is identically unity so that a; is 
identically zero. Hartree‘ presented the maximum wall shear 
solutions. The minimum wall shear solutions for this problem 
have a negative shear at the wall (further details are given in 
reference 3). 

The phenomenon of velocity ratios in excess of unity in the 
outer portion of the boundary layer has been termed ‘‘velocity 
overshoot”’ in reference 3. Two distinct types of velocity over- 
shoot are possible. The first of these occurs when az ¥ 0 (and 
38 < 0) and therefore must be excluded. The second occurs in 
favorable pressure gradient flow over heated surfaces. The exist- 
ence of this type of velocity overshoot can be demonstrated analy- 
tically. For favorable pressure gradients, (8 > 0) with ae set 
equal to zero as previously mentioned, the a; term in Eq. 1 is dom- 
inant for large ». Thus (f’ — 1) and a; are of the same sign 
Hence for a heated wall (a; positive, Eq. 2), the velocity ratio 
must approach unity from above. This excess velocity is as 
sociated with the acceleration of a layer of fluid in the outer por 
tion of the boundary layer with density less than the external den 
sity. Since this layer is subject to the external pressure field and 
is restrained only slightly by the viscous forces acting on it, it is 
accelerated more than the external flow. This velocity over- 
shoot may in certain circumstances be large enough to produce 
negative momentum thicknesses. This is completely unrelated 
to the phenomenon of negative displacement thickness discussed 
in reference 2. The displacement thickness can become negative 
for large rates of cooling with favorable pressure gradient because 
of the high density within the lower portions of the cooled 
boundary layer. 

With regard to the integral technique for the calculation of 
laminar boundary layers proposed in reference 1, it is surprising 
that it does not agree with Rott and Crabtree5 The integral 
technique devised by the present writers for compressible laminar 
boundary layer with heat transfer and pressure gradient? does 
reduce to that of Rott and Crabtree for the insulated surface 
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Engineering Relations for Friction and Heat 
Transfer to Surfaces in High Velocity Flow 


E.R. G. Eckert * 
Professor of Mechanical Engineering, University of Minnesota 
March 23, 1955 


SYMBOLS 


friction factor 
specific heat at constant pressure 


h = heat-transfer coefficient 
i = enthalpy 
k = thermal conductivity 
q = heat flow per unit area and time 
r = temperature recovery factor 
ri = enthalpy recovery factor 
x = distance from leading edge 
7 absolute temperature 
V = velocity 
Y = ratio of specific heats at constant pressure and volume 
“ = viscosity 
p density 
r = shearing stress 
Re = Reynolds Number 
Ya = Mach Number 
St = Stanton Number 
Pr = Prandtl Number 
Subscripts 
i based on enthalpy 
m maximum 
r recovery 
§ in stream outside boundary layer 
u at wall surface 
* reference condition 


DISCUSSION 


IT THE DESIGN OF MISSILES for large supersonic Mach Numbers, 
special consideration must be given to the temperature which 
the skin of the missile assumes. This requires calculations de- 
termining the heat transfer and skin temperature for a consider- 
able number of flight conditions. Such calculations are very 
time-consuming unless simple relations are available which 
describe with sufficient accuracy the friction and heat transfer 
to surfaces in high velocity flow. It has been pointed out pre 
viously that values for those parameters calculated with rela- 
tionships which hold for constant property fluids describe actual 
conditions in a supersonic boundary layer in many cases satis 
factorily when the properties of the fluid are introduced at a proper 
reference temperature (see references 1, 2). This is in itself sur- 
prising because the properties vary for large supersonic velocities 
extensively throughout the boundary layer. The validity of this 
procedure, however, has been verified again recently in a survey 
(see reference 3*) on heat transfer at high velocities in which a 
relationship for the reference temperature was developed. The 
introduction of the properties at this reference temperature 
describes friction factors and heat-transfer coefficients on a flat 
plate with constant wall temperature in such a way that it agrees 


* Permission of the Wright Air Development Center, Air Research and 


Development Command, USAF, for publication is gratefully acknowledged. 
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it high temperature. Prandtl Numbers as given in the Gas 
Tables and as recommended by the National Bureau of Standards 
jiffer, for instance, at a temperature of 2000°R by 12 per cent 
and accordingly the recovery factors calculated with those 
Prandtl Numbers differ by 6 per cent and heat-transfer coefficients 
by 8 per cent 

TURBULENT BOUNDARY LAYER 


Good agreement with available experimental information is 
iso obtained when the calculation procedure outlined above is 
used and when appropriate relationships describing friction 
factor, recovery factor, and heat-transfer coefficient in incom 
pressible turbulent flows are introduced. A relation which 
is in common use to describe the turbulent friction factor of a 
flat plate is the equation by Karman-Schoenherr.* Very good 
wreement with experimental results in the whole Reynolds 
Number range is also obtained by a relation for the local friction 
factor which has been published by Schultz Grunow. 


cr = 0.370/(logiRe*)?-584 (15 


The temperature or enthalpy recovery factor can be de 
scribed with good accuracy by the equation 


r= \/Pr* 16 
ind the Stanton Number by 
St = (cy/2) (Pr*) 2/8 (17 


The last relation has been proposed by A. P. Colburn and seems 
to be the one which also satisfactorily describes newer experi 
ments for air and generally in a Prandtl Number range not too 
far away from the value one. A comparison between friction 
factors determined by the proposed method and values obtained 
from experiments is shown in Fig. 1. In this figure the ratio of 
the actual friction factor cy in a supersonic boundary layer of air 
to the friction factor cy; for an incompressble fluid is plotted over 
the Mach Number. The friction factors in the group cy/cys; are 
not based on properties at the reference temperature but on the 
properties at the static temperature in the stream outside the 
boundary layer, a procedure which is frequently used in the 
presentation of such data. The points indicate the results of ex 
perimental investigations (for references see reference 3) and the 
lashed lines give the relations obtained from the above described 
alculation procedure. It can be recognized that very good 
igreement is obtained not only for the case where the wall is 
unheated (the wall temperature equal to the recovery temper 
iture) but also for a wall temperature equal to the stream tem 
perature. The two experimental points at a Mach Number 7.2 
ind for the condition 7, = 1.8 7; are somewhat lower than the 
alculated curve. However, it has been pointed out in the cor- 
responding publication that the accuracy of these tests is not 
is high as the one for the measurements at Mach Number 3.8 and 
Ty = T;. Information on heat transfer is more scarce than on 
friction. There are, however, indications that the recommended 
procedure describes also the conditions for heat transfer satis 


factorily. No experimental information is available at large 


Mach Numbers, at large temperature differences, and at high 
temperatures, and the above recommendations have to be con 
sidered as tentative for these ranges 

It should be pointed out once more that the relations presented 
here are restricted to constant pressure along the surface, to two 
dimensional flow, and to constant wall temperature. They can 
be used to obtain information on rotationally symmetric flow 
when an appropriate transformation is made and also as a basis 
to calculate friction and heat transfer for a wall temperature 
which varies along the surface. The relations are also restricted 
to the condition in which the boundary-layer equations describe 
friction and heat transfer satisfactorily, especially for conditions 
where no slip or boundary-layer shock interaction occurs and 
where the air within the boundary layer is not dissociated as a 
Consequence of high temperatures. The slip condition restricts 
the relation to flight below approximately 30 miles. The maxi 
mum temperature to be expected within the boundary layer can 
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be determined from Fig. 2. From this figure it can also be esti 
mated whether dissociation in the boundary layer will occur 
In this figure the maximum temperature 7,, within the boundary 
layer is plotted over the Mach Number for flow along a flat plate 
When the surface is not cooled or heated then the recovery tem- 
perature is the largest temperature within the boundary layer 
Usually the surface temperature of missiles will be kept at lower 
values either by heat storage in the skin or by radiation. For 
two temperature conditions the maximum temperature within 
the boundary layer is also indicated in Fig. 2. It may be ob 
served that for a ratio of the absolute wall temperature 7, to the 
static temperature in the stream 7, equal to 4 dissociation must 
not be expected at an altitude of around 50 miles below a Mach 
Number around 14 and for zero miles altitude for Mach Num 
bers up to 18. Recent calculations’ * have also indicated that 
i certain amount of dissociation will not effect heat transfer and 
friction when the dissociation does not occur in the layers im 
mediately adjacent to the surface. This would indicate that 
the recommended relations can be used even up to higher Mach 


Numbers 
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| apa ED INTEREST has been shown recently in these columns 
in the problem of calculating the flow in the transonic part 
of two-dimensional supersonic wind-tunnel nozzles.’ * It ap 
pears, however, that an earlier paper by Cherry® giving the solu 
tion of the exact hodograph equation for the nozzle problem in a 
practical form has escaped the knowledge of aerodynamicists 
It may therefore be desirable to point out briefly the usefulness 
of Cherry's method of nozzle design. 

Cherry’s solution is in the form of a series of eigenfunctions of 
the hodograph equation, the first of which possesses the par 
ticular sonic singularity of nozzle flows with nonzero pressure 
gradient at the throat. The solution for any prescribed analytic 
distribution of pressure along the axis of symmetry can be ob 
tained by a suitable choice of the coefficients of the series 

The pressure distribution on the axis of symmetry need not, 
however, be prescribed in detail. Each eigenfunction is an exact 
solution of the hodograph equation, and a finite number of them 
suffice to meet all physical requirements of nozzle design. A 
nozzle for f = 2.41 has been designed for the Sydney University 
Supersonic Tunnel from computations kindly supplied by Pro 
fessor Cherry, and it was found that only two eigenfunctions 
were required! The Mach Number distributions on axis and 


liner are shown in Fig. 1 
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The low-speed contraction was completed, for 1/ < 0.34, by 
methods of incompressible, three-dimensional duct design. A 
the other end of the nozzle, the hodograph solution was joined 
to the simple wave region by a narrow transition strip computed 
by the conventional method of characteristics, in order to avoid 
any discontinuity of pressure gradient and liner curvature; this 
transition is indicated by the dashed segments of the curves 
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